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Flutter Analvsis of Complex Airplanes by 
EK xperimental Methods 


E. BERKELEY KINNAMAN* 
Boeing Airplane Company 


ABSTRACT 


The analysis of the flutter potentialities of a highly elastic air- 
plane having a wing of high aspect ratio with large flexibly mcun- 
ted concentrated masses distributed along the span has proved to 
be amenable to solution by the use of elastically and dynamically 
scaled wind-tunnel models. Experimental flutter testing at Boe- 
ing has progressed from a simple half-wing cantilever model of 
1946 to today’s almost completely scaled flutter models tested in 
reasonable simulation of free flight. 

Continued development of experimental flutter analysis tech 
niques and correlation with full-scale flight testing and theory are 
needed to permit more accurate predictions of the flutter charac 
final design and flight 


teristics of an airplane prior to the 


stages. 


INTRODUCTION 


T \HE PREVENTION OF THE DYNAMIC INSTABILITY of an 
“flutter,” is approaching ‘“‘big 
industry 


airplane, called 
business’ in the 
agencies. 


aircraft and research 


Primarily, this field of engineering is preventive medi- 
cine for the high-performance aircraft now in progress of 
design, since flutter engineers are required to ensure 
that flutter is not designed into the airplane. 

Being certain that a complex airplane will be free 
from flutter is becoming increasingly more difficult if 
not impossible. The growth of the facilities available 
to the flutter engineer have not in fact been so rapid 
as the increase in design problems. At Boeing, we 
have been faced with the problem of wing-supported 
concentrated masses from the initiation of the airplane 
design. We have been faced with problems such as low 
frequencies, highly swept wings, variable density, and 
interaction of components. 


Presented at the Aeroelasticity Session, Twentieth Annual 
Meeting, I.A.S., New York, January 28-February 1, 1952 
‘Vibration Unit Chief 


PRELIMINARY FLUTTER ANALYSIS 


The approach to the analysis of the flutter character- 
istics of an airplane design is similar to most other engi- 
neering problems. Ideally, we need to answer rapidly 
the preliminary design engineers’ questions, first, on 
the probability of flutter of the airplane configuration 
as a whole and, second, as to the effects, flutterwise, of 
all the variations that show promise of increasing the 
efficiency of the flight article. 

Obviously, long and involved calculations are im- 
practical in this early stage. In most instances, the 
flutter analyst can do no more than call on intuition 
based on experience and on the available literature to 
recommend preferable assemblies of components. As 
the design begins to take shape and coordination be- 
tween the aerodynamics, power plant, structures, and 
equipment has resolved the airplane into one general 
pattern with several variations, the usefulness of a more 
complete flutter investigation of the airplane and com- 
ponents becomes apparent. It is at this point that the 
need for a facility to calculate the flutter characteristics 
in an expeditious manner is a necessity if the airplane 
design is to benefit from these calculations. The desk 
calculator, IBM, differential analyzer, and analog all 
have features that limit their usefulness. We have no 
experience with the high-speed digital computer; how- 
ever, we are hoping that it will markedly reduce the 
number of limitations. 

Recognizing the fact that engineering calculations 
start with significant simplifying assumptions, it then 
becomes a matter of choice of assumptions necessary to 
accomplish the task and selection of the best analytical 
procedure. 

A flutter engineer’s definition of an airplane is con- 
siderably different from that of other aeronautical engi- 


neers—naimely, an infinite number of masses connected 
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by springs and dampers in a geometrical pattern, each 
mass being subjected to aerodynamic, elastic, acceler- 
ating, and damping forces. 

The necessary parameters for inclusions in the flutter 
analysis must be resolved from the drawings into num- 
bers. 


Immediately the simplifying assumption of a 
finite number of masses (whether 5, 10, or 1,000) is 
This is done by arbitrarily selecting seg- 
ments of the wing, fuselage, and empennage surfaces as 


necessary. 
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the mass and aerodynamic area parameters. Since 
these segments are structural compenents of the air. 
plane, they also comprise the elastic structure. 
Determination of the mass characteristics of the seg. 
ments requires skilled weight engineers to calculate the 
detail ingredients to the required accuracy. The elastic 
properties of the complex structure built up from thin 
sheet, stiffeners, multiple spars, etc., present a more 
difficult analysis problem. On the assumption that 
primary modes of the components are the important 
ones, it is possible to consider that the simple E/ and 
JG stiffnesses Available information 
indicates that this may be oversimplification due to 
effects of shear lag, discontinuous structure, and redun- 


are sufficient. 


dancies. Rigorous analysis and consideration of higher 
order modes would require the inclusion of all of these. 
However, since these details are not available until 
well into the detail design stage of the airplane, it is 
necessary to proceed with the simpler stiffness data. 
The next decision required is that regarding the aero- 
dynamic terms. Cross-coupling of air forces between 
segments, tip effects, compressibility, and downwash on 
the horizontal tail surfaces all require consideration 
prior to determining both the method of making the 


flutter calculation and the facility to use. 


CHOICE OF COMPUTING FACILITY 


After reviewing the assumptions required to conduct 
either mathematical or experimental flutter analyses, it 
becomes apparent that there are several advantages in 
exploiting the mechanical analog or wind-tunnel flutter 
model. The model can 
elastic terms for all components as would be used in the 
mathematical calculations and, within the limitations 
of the test air velocity, all aerodynamic terms. The 
model makes flutter calculations so rapidly that numer- 
ous changes in configuration (six to 18 per test day) can 
This makes it possible to determine 


include the same mass and 


easily be made. 
rapidly the flutter sensitivity of many of the influential 
parameters. 

At Boeing, we began our work on experimental low- 
speed wind-tunnel models several years ago, with the 
expectation of obtaining a check of the calculated 
flutter characteristics of an airplane. It appeared ad- 
visable to make the model include the same assumptions 
used in the flutter calculations, which included the rigid 
body freedoms; wing vertical bending and wing torsion 
stiffnesses; and infinite stiffnesses of the fuselage, em- 
The 
calculations indicated that flutter would occur at 120 


pennage, and engine nacelle attaching struts. 


per cent design speed in an antisymmetrical mode. 
The model agreed that the flutter speed was approxi- 
mately as predicted analytically; however, it fluttered 
in a symmetrical mode. This was one of the first indi- 
cations that the analytical approach might have major 
There began the struggle to reduce the 
number of assumptions, refine the analyses, and include 


deficiencies. 
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additional aerodynamic elastic and mass terms as 
additional degrees of freedom. 

Our models have become increasingly complex as 
we have gained more experience with the various items 
influencing flutter modes. We are finding that the 
number of test configurations required to satisfy even 
the minimum of possible couplings of parameters is 
continually increasing. 

We initially tested the easiest changes, such as span- 
wise and chordwise location of wing-mounted concen- 
trated weights and only one airplane weight that ex- 
perience and intuition had said was the most critical. 
We are now testing 14 or more gross weight configura- 
tions involving various distributions of useful weight 
and variation in stiffness and stiffness ratio of the wings 
and nacelle struts, thus multiplying the number of tests. 

The present models scale the elastic, mass, and geo- 
metric dimensions for horizontal and vertical tail sur- 
faces, fuselage, wings, and nacelles and struts. Weare 
not including movable control surfaces on our present 
models because of mechanical difficulties, as well as the 
questionable character of the low Reynolds Number of 
the small chord surfaces. Simulation of the jet stream 
of the engines has not been attempted; however, the 
gyroscopic effect of the engines has been included for 
some cases. Simulation of different altitudes has been 
made by incorporating the true scale of model-to-air- 
weight ratio and to obtain the proper load factor wing 
deflection scale. 


PARAMETER SENSITIVITY 


Some flutter modes have régimes of variation of 
parameters where extreme sensitivity is shown, and 
small changes result in large changes in flutter speed. 
This has been shown many times by the models. For 
example, a 20 per cent increase in wing torsional stiff- 
ness resulted in more than 100 per cent increase in 
flutter speed; 5 per cent reduction in outboard wing 
fuel raised flutter speed more than 100 per cent; 3 per 
cent reduction in fuel in external fuel tank reduced 
flutter speed by 20 per cent with a major change in 
mode. Figs. 1 and 2 show the abrupt change in flutter 
mode and speed for two variables. One mode was 
sensitive to relatively minor changes in model assembly, 
such as slight friction contact between the balsa wing 
sections and lead weights simulating fuel; small angles 
of yaw, 1° or less, of the external fuel tank; change from 
a small tubular weight to a streamlined tank; or the 
addition of a small amount of liquid fuel in the tank. 
An investigation of the influence of rotational fre- 
quency of an all movable horizontal tail surface re- 
sulted in wide changes of mode shape. Fig. 3 shows 
this variation. Investigations of the nacelle strut 
stiffness have shown this to be a significant variable. 

Significant changes in wing flutter also have occurred 
with variations in fuel weight distribution in the fuse- 
lage. 
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Fic. 4. 


Typical wing spar design. 


One interesting and significant result to which we had 
not previously given sufficient importance was experi- 
enced during a series of tests of identical configurations 
on models scaled for simulation of two different alti- 
tudes. Flutter curves at one altitude were shifted a 
considerable amount from the curves obtained at the 
other altitude for one series of parameter changes. 
Where one flutter mode had been the critical one, an- 
other mode became critical at a lower speed for the 
second altitude. In addition to the significance of the 
air density change, the problem of being sure that all of 
the necessary parameter variations had been investiga- 
ted is highlighted by these results. 


MopeEL WING DESIGN 


The above phenomena have been experienced over a 
period of several years and the testing of a number of 
models, and they appreciably influenced the details 
of our model construction. A good example of this was 
the development of our wing design. The first model 
wing was conceived as the simplest possible structure. 
The choice of scale parameters, geometry to be as 
large as possible for ease of manufacture and speed 
scale to permit simulation of the maximum airplane 
speed within the speed range of available tunnels, re- 
sulted in a mass and stiffness scale that permitted the 
simple wing design. Tests showed that a wing of 
solid balsa could, with suitable spanwise and chordwise 
slots, be made to scale the wing vertical bending and 
torsional stiffnesses. It was assumed that scaling the 
fore-and-aft stiffness was not necessary. We were 
fortunate that, by trial-and-error slotting, the first 
set of wings easily achieved the proper stiffness. Sub- 
sequent wings were more difficult. After it appeared 
that we would be expending more effort in experimental 
model testing, we decided to develop the structural 
design so that the proper stiffnesses would be predicted 
rather than rely on trial-and-error procedure. A spar 
of rectangular cross section, tapered spanwise as neces- 
sary and with flanges to control stiffness in that plane, 
was selected and resulted in excellent accuracy in 








580 JOURNAL OF THE 





Side view of initial soft spring support 


Fic. 5. 


stiffness and stiffness distribution, as seen in Fig. 4. 
The flutter modes of the first model with this wing 
design involved large fore-and-aft displacements of the 
wing tip. Review of the stiffness data for the model 
and airplane uncoverd an oversight in computation of 
the fore-and-aft wing stiffness which had resulted in 
extra flexibility for the model. Correction of this 
deficiency raised the flutter speed approximately 20 per 
cent and reduced the fore-and-aft deflection to a small 
magnitude. 

Consideration of the significance of the fore-and-aft 
stiffness leads logically to the recognition that the 
static vertical bending deflection is also important 
flutterwise because of the coupling of fore-and-aft 
wing bending with wing torsion. The additional 
aerodynamic coefficients of the fore-and-aft movement 
can also be effective. To obtain the scaled wing de- 
flection, the static aerodynamic span load distribution 
must be scaled unless preformed spars are used. It also 
follows that, if the aerodynamic coefficients are non- 
linear, the wing must be operating at the same Cy, as 
the full-scale airplane. 

We have been flying the models at unit load factor 
and in most cases designing for true scale deflection. 
If it is desired to investigate other than unit load factor 
wing deflection, then additional wing spar stiffnesses 
are required for simulating stiffness changes due to 
skin wrinkling or to obtain the desired deflection curve 
unless some loading provision that will not interfere 
with the flutter mode is used. 


MOopDEL SUPPORT 


The support of the dynamically and elastically scaled 
model is a major problem in a successful wind-tunnel 
flutter model program. Our suspension systems have 
varied considerably in attempting to achieve the de- 
sired boundary conditions. Provision for the rigid 
body motions for symmetrical or antisymmetrical 
simulations, ensurance that the flutter modes are 
influenced as little as possible by the supporting means, 
prevention of low-frequency unstable oscillations of the 
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model, and safety provisions to prevent destruction 0 
the model in case of flutter are some of the design re 
quirements. 

Our first support installation was chosen after review 
ing the available information of flutter model testing 
by others and was patterned on the spring suspension 
used by Junkers in Germany. Fig. 5 shows this in. 
stallation. The model was supported on music wir 
attached to soft cantilever springs at the floor and ceil 
ing. Multiple wires attached to the wing and fuse. 
lage leading to a manual pull-down, which pulled the 
wires taut to restrain the model firmly. A single drag 
wire was attached to the nose of the fuselage. Two 
antidrag wires diagonally aft to each side and attached 
to the tunnel wall by soft springs positioned the mode 
laterally. This resulted in elastic restraint in all rigid 
body rotations and translations, except fore and aft 
which was rigid. The fish net stretched across the tun- 
nel was to catch the broken pieces of the model as it dis. 
integrated during flutter. The first test run on this in- 
stallation introduced us to one of the most troublesome 
problems of complete model testing—that of instability 
in pitch and vertical translation. Several trial-and-er- 
ror changes were required before it became apparent 
that, even though the vertical restraining springs were 
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Fic. 6. Front view of initial soft spring support 
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Fic. 7. Plan view of initial soft spring support 
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extremely soft, their fore-and-aft location with respect 
to the airplane center of gravity was important. Hav- 
ing the elastic center behind the c.g. is beneficial. 
Figs. 6 and 7 are sketches of the most successful ar- 
rangement. Rigidity in side translation provided by 
the cross drag bracing was used to help control Dutch- 
roll instability. 

A second approach was to design independent sup- 
ports for symmetrical and antisymmetrical testing. 
Fig. 8 is a side view of the symmetrical installation with 
the model attached to the end of a rigid strut which 
was free to travel vertically on four roller bearings. 
The model was free to pitch on a bearing at the strut 
attachment located at the model c.g. Provision for 
large vertical displacement of the model approximately 
40 in. was incorporated on the premise that the model 
instability difficulties were primarily due to the limited 
displacement and the nonlinear galloping of the model 
off the In addition, a hydraulic shock 
absorber was designed for limiting vertical translation 
and bounce-off. In spite of these precautions, consider- 
able difficulty was experienced with instability at high 
air speeds for many of the test configurations. 

The extra weight of the strut was a disadvantage in 
that it was impossible for the lighter gross weights to be 
Fig. 9 is a photo 


snubbers. 


scaled for true weight distributions. 
of the strut and upper bearing support. 
the attachment of the strut to the fuselage spar and 


Fig. 10 shows 


pitch bearing. 

Fig. 11 is a side view of the antisymmetrical support 
providing fixity in drag and vertical translation, free- 
dom in pitch and yaw, and elastic restraint in roll and 
side translation. Instability in roll and side transla- 
tion made it necessary to restrain elastically the model 
in side translation to obtain a favorable frequency ratio. 
Positioning of the modei in roll was accomplished by a 
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Fic. 8. Side view of strut support for symmetrical modes. 
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Fic. 9. Detail of strut and bearing support 





strut attachment. 


Fic. 10. Detail of model 
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modes 


Fic. 11 


remotely operated electric motor rotating a torque 
tube attached to the trunnion and fuselage spar. Fig. 
12 gives additional details of this installation. 

One test result not anticipated was the sensitivity of 
the antisymmetrical flutter mode to the fore-and-aft 
position of the trunnion attachment on the body spar; 
moving from 10 to 18 per cent M.A.C. raised the flutter 
speed by 110 per cent. 

A major fault in the use of two independent installa- 
tions for symmetrical and antisymmetrical investiga- 
tions considerably devaluates their usefulness. The 
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Details of side translation mount. 
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Fic. 13. Side view of vertical rod installation 
artificial locking of the rigid body motions change the 
flutter characteristics of the model for the condition 
not being studied. For example, many cases of sym- 
metrical flutter occurred with the model mounted on 
the antisymmetrical mount at speeds considerably 
lower than the required maximum speed. 

The flutter model suspension system now being used 
for our testing has improved features. Fig. 13 is a 
schematic side view of the vertical rod installation. 

A !/,-in. diameter rod, ground to uniform diameter 
and hard-chrome plated, is installed under controlled 
tension. A short trunnion sliding on the rod and 
attached to the fuselage spar through a pitch bearing 
located at the model center of gravity provides all of 
the rigid body degrees of freedom; zero frequencies in 
vertical translation, pitch and yaw; elastic restrain 
provided by the rod, in roll and lateral and fore-and- 
aft translation at approximately one-half the frequency 
of the lowest antisymmetrical mode. The shock 
absorbers at the limits of travel on the rod are secondary 
safety stops on this installation. The flexible cables 
attached to the trunnion above and below the model 
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lead outside the tunnel to a servomechanism to suppress 
the vertical translation-pitch instability. A low-ire- 
quency pass filter reduces the effect of the cable re- 
straint, although the cables are slack except when needed 
as controlled by the intelligence of the mechanism. 
Our automatic damper is a skilled, dexterous, and hard- 
working engineer holding the ends of the cables. Oscil- 
lograph records of his force-time-displacement show 
good 90° phase. The design of an electromechanical 
servomechanism is on our list of work to be done, since 
our present models sometimes demand the maximum 
of frequency and force response from our human servo- 
mechanism. 

Locks for stopping the flutter to prevent damage to 
the model are not used. It has been found that the 
models are sufficiently rugged to withstand the majority 
of the flutter cases. Numerous minor failures have 
occurred during the developmental period; however, 
only one major failure has resulted from flutter. A 
high-frequency flutter mode of the horizontal tail 
surface failed an underdesigned spar fitting, causing 
The resultant divergence of 

The time for failure was so 


loss of the whole surface. 
the model failed the wings. 
short that the operation of any safety system would 


have been improbable. 








Fic. 14. Vertical rod installation. 
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FLUTTER 
Fig. 14 shows the tunnel installation. The faired 
island below the model surrounds fixed instrumentation 
in the tunnel not associated with the flutter testing. 
An attempt was made to remove the side translation 
elastic restraint and to reduce the roll spring rate by 
the installation shown in Fig. 15. This was a combina- 
tion of the vertical rod support and a pinended link. 
Dutch roll occurred at extremely low speeds and could 
not be eliminated by reasonable measures. This system 
was quickly abandoned. 


MopEL CONSTRUCTION 


The details of the model construction have developed 
in the same manner as the model supports. The solid 
balsa structure of Fig. 16 has changed to hollow boxes 
covered with balsa and silkspan tissue paper. Fig. 17 
is one design of a wing segment. The single point 
attachment at the spar prevents the stiff box from af- 
The slots between segments 


fecting the spar stiffness. 





Fic. 16. Solid balsa wing and fuselage 
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Fic. 18. Model details 
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Fic. 19. Model details of plastic liquid-containing tanks 





Empennage assembly. 


Fic 20 


are covered with thin rubber strips cemented on the 
final assembly. Control of weight, center of gravity, 
and pitch moment of inertia per segment is achieved by 
actual weighing, balancing, and swinging on a bifilar 
pendulum, then adding the necessary ballast weights 

In Fig. 18, part of the fuselage shell has been re- 
moved, revealing instrumentation leads, fuselage spar, 
and trunnion details. Fig. 19 shows liquid-containing 
tanks for fuel simulation; Fig. 20, the empennage. 
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The small electric motor, worm and gear, and screw 
attached to the leading edge of the stabilizer and fuse- 
lages spar raises or lowers the stabilizer angle for flying 
the model in the tunnel. Aileron and rudder control 
has not been necessary. 

Since structural damping can be an important param- 
eter, the model construction must be closely controlled 
to prevent rubbing friction or the use of materials with 
high hysteresis. Our design is easily controlled, and the 
resultant model damping is conservatively less than 
that of the airplane. 


FUTURE DEVELOPMENT 


experimental flutter 
The use of a 


Continued development of 
analysis is required along several lines. 
variable density wind tunnel to investigate further the 
influence of altitude on flutter modes is needed. Flut- 
ter testing of a complete scaled model at a one-to-one 
speed scale further to reduce aerodynamic assumptions 
would be highly beneficial. It is necessary to continue 
investigation of the influence of tunnel walls on oscil- 
lating air forces. The simulation of the elastic prop- 
erties of the full-scale wings, especially the lower as- 
pect ratios, is a good research problem. 
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A major gap exists in the knowledge of how well the 
models or mathematical analyses simulate the full-scale 


airplane flutter characteristics. Full-scale flight-flut 
ter testing of a number of configurations, which have 
been previously investigated in the wind tunnel, would 
be necessary to provide the minimum information 
It should be emphasized that negative flight results 
would not serve this purpose. Testing techniques and 
safety measures would have to be developed to _ per- 
mit obtaining the required information with a near 
approach to flutter. 


CONCLUSIONS 


Experience has shown that flutter analysis by use of 
the mechanical analog or wind-tunnel flutter model 
is an expeditious means for investigating the flutter 
characteristics of a complex airplane. 

Continued effort to develop this technique and to 
reduce the number of significant assumptions of elastic 
simulation, aerodynamic compressibility, and tunnel- 
wall effects will be well worth while. Correlation of 
flutter analyses and full-scale flight-flutter investiga- 
tions is necessary to assist in improving the accuracy 
of flutter predictions. 
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Analysis of Oval Rings by Fourier Series’ 


ROLAND J. BENJAMIN? 
l’niversity af Illinots 


SUMMARY 


The distinguishing feature of this method of analysis of thin 
oval rings is the representation of the shape of the ring by the 
use of Fourier series. By means of standard computational 
schemes, the coefficients of Fourier series of sine and/or cosine 
terms may be so selected as to cause the graph of the function 
represented by the series to approximate closely the shape of 
the ring under consideration. 

Once the ring has thus been represented by an analytical 
expression, any of the common strain energy methods may be 
employed in the solution of the problem. Because of its sim 
plicity, the method of elastic weights is used herein, and a 
sample problem of a general oval ring subjected to concentrated 
loads is worked out. 

The possibilities of application of the method to fuselage ring 
design are pointed out by the solution of the problem of an 
elliptical ring subjected to a concentrated load and held in 
equilibrium by distributed skin forces. 

The accuracy of the method, established by analyzing true 
To 
reduce computation, some of the higher harmonics of the Fourier 
series may be neglected with but negligible effect upon the result 


elliptical rings whose solutions are already known, is high. 


The method is general and may be applied to rings having 


two, one, or no axes of symmetry 


INTRODUCTION 


_ PROBLEM OF ANALYSIS of a thin closed ring 
occurs frequently in the design of aeronautical 
structures. For rings of circular or elliptical shapes, 
numerous solutions are available in the literature.!~° 
The more general case of rings with but one axis of 
symmetry has been treated by several writers®~* in 
broad terms, but in the applications of the proposed 
methods of solution to specific rings (other than cir- 
cular) difficulties are encountered in attempting to 
For such rings some 
usually 


evaluate the integrals required. 


method of graphical integration must be 
employed. 

The method of solution proposed in this paper avoids 
such procedures and permits solution of the problem 
It is applicable to rings having 
two, one, or no axes of symmetry. However, the case 
of complete lack of symmetry is of little importance, 


since the rings encountered in aircraft design, particu- 


on an analytical basis. 


larly fuselage rings, normally have at least one axis of 
symmetry. 
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The basic step in this method is the analytical 
expression of the coordinates of a point on the ring by 
means of Fourier series. 


DETERMINATION OF FOURIER SERIES 

To 
locate points on the ring, the vertical axis of symmetry 
is chosen as the y-axis. The x-axis may be taken 
through any convenient point; in this case it is midway 


Let us now consider the ring shown in Fig. |. 


between the extreme ends of the ring. 

It is desired to express the coordinates (x, y) of any 
point on the ring in terms of the parameter s, where 5s 
is the distance along the curve measured clockwise 
from the top of the ring. In order to obtain simple 
periodic functions, it is convenient to replace s by the 
parameter @ defined as 


6 = ms/c (1) 


where c is the semicircumference of the ring. Thus, as 
the ring is traversed from the uppermost point in the 
clockwise direction, s increases from zero to 2c and @ 
increases from zero to 27. It should be pointed out 
that @ is merely the quantity defined by Eq. (1) and 
is not the measure of any geometric angle. 

We wish to find expressions for x and y as functions 


of in the form of Fourier series. In the computational 








Oval ring for Fourier analysis 


Fic. 1. 
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TABLE 1 
Coordinates of Points on Ring of Fig. 1 

Point S x y 
0 0 0) 10.00 
l e/iZ 2.20 9.52 
2 c/6 4.01 Ss. ir 
3 c/4 5.25 6.27 
4 c/3 6.20 4 20 
5 5c/12 6.83 1.99 
6 c/2 7.16 —().28 
4 7c/12 7.14 —2.56 
8 2c/3 6.80 —4- 93 
9 3c/4 5.83 —6.90 
10 5c /6 4.25 —8.54 
1] l1c/12 2.24 —9.61 
12 ( 0 —10.00 
13 13c/12 —2.24 —9 61 
14 7c/6 —4.25 —8. 54 
15 ne /4 —5.83 —6.90 
16 40/3 —6.80 —4.83 
Li 17¢/12 —7.14 —2.56 
18 3c/2 —7.16 —0.2 
19 19¢/12 —6.83 1.99 
20 5c/3 —6.20 4.20 
21 7c/4 —5.25 6.27 
22 llc/6 —4.01 8.17 
23 23c¢/12 —2.20 9.52 
24 2c 0 10.00 


scheme to be used, the coordinates of equally spaced 
points around the ring are required; it is particularly 
convenient to use 24 such points. These values, as 
well as the value of c, may easily be determined from a 
drawing of the ring with the aid of dividers. The 
values of x, y, and s for the various points are tabulated 
in Table 1. The value of c was found to be 27.46. 
Evidently x is an odd function of 6 and y is an even 
Hence, x may be represented by a sine 
Following the method 


function of @. 
series and y by a cosine series. 
of Fourier analysis described by Den Hartog’® and 
using the values in Table 1, we get 


x = (1/12)(9.06 sin @ — 2.97 sin 26 + 4.40 sin 30 + 
1.09 sin 40+ 0.22 sin 506+ 0.46 sin 60 + 
0.10 sin 76 — 0.15 sin 86+ 0.04 sin 94 — 
0.05 sin 108 — 0.14 sin 116) (2) 
y = (1/12)(—2.57 + 115.91 cos 6 + 1.65 cos 26 + 


3.94 cos 36 + 1.04 cos 40 + 0.17 cos 56 + 
0.04 cos 66 + 0.01 cos 78 — 0.16 cos 86 — 
0.06 cos 99 — 0.01 cos 108 — 0.03 cos 116 + 
0.01 cos 126) (3 


APPLICATION TO OVAL RING 


Since we now have expressions representing the shape 
of an arbitrarily chosen oval ring, we may proceed to 
find the bending moment at any section of the ring due 
As an example, let the 
The vertical forces 


to the action of applied loads. 
ring be loaded as shown by Fig. 2a. 


P/2 located at the points s = 2c/3 and s = 4c/3 are 
resisted by the single force P acting at the top of the 
ring. 


Any convenient strain energy method of analysis 
The method of elastic weights® will 
“dummy- 
In either 


may be followed. 
be used in the illustrative examples; the 
load” method has been used with equal ease. 
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case it is assumed that Hooke’s Law is valid and that 
the strain energy due to axial and shear forces is 
negligible as compared with that due to bending. 

Since both the ring and its loading are symmetrical 
about the vertical centerline, the ring may be assumed 
to be cut at the top and bottom and only one-half of it 
need be considered, as in Fig. 2b. 

In applying the method of elastic weights the ring 
is assumed rigidly fixed at the top section. Instead 
of applying the indeterminate moment and forces 
directly to the cut section at a, they are assumed to be 
acting at the end of a rigid arm L rigidly attached to 
the ring. 

An elastic weight dw is defined by 


dw = ds/EI (4 


where ds is a differential length of the ring. E and / 
have their usual meaning and are here considered to 
be constant throughout the ring. 

The length of the fictitious arm L is such that its 
free end will fall on the origin of coordinates. This 
origin must be located at the “center of gravity of the 
elastic weights’ of the complete ring such that 


Sxdw= fxydw= fydw=0 
or, if EJ is constant, 
Sxds= fxyds = fyds =0 (5 


Since the ring is symmetrical about the y-axis, the 
first two conditions of Eq. (5) are automatically satis- 
fied. To satisfy the third condition, the vertical loca- 
tion of the origin, which was arbitrarily chosen to 
carry out the Fourier analysis, must be changed. Let 
us denote by y’ the ordinate of a point referred to the 
present axis and by y the ordinate of the same point 
referred to the axis through the c.g. of the elastic 
weights. 

Considering the present axis, 


So’ y' ds = cF (6 


where c is the semicircumference of the ring and 7 is 
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Fic. 2. Probiem of oval ring. 
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the ordinate of the c.g. of the elastic weights. If 7 is 
now found and the origin translated to the point (0, 9), 
the condition prescribed for the use of the elastic 
weight method will have been fulfilled 

The expression for the quantity now denoted by y 
(3) in terms of 6. Recalling that 


, 


is given by Eq. 
§ = ms/c, we have 
mds c 
dé = or ds = -d0 
c T 
Then the integral of Eq. (6) becomes 


. 


0 us 0 


This may easily be evaluated by term-by-term inte- 
gration of Eq. (3), giving 


Thus, 


and the origin must be transferred to the point 
0, —2.57/12). The effect of this transfer is to elimt- 
nate the constant in the expression for y’, and there 


results 


12)(115.91 cos 6 + 1.65 cos 26 + 3.94 cos 30 + 
1.04 cos 46 + 0.17 cos 56 + 0.04 cos 64 + 
0.01 cos 76 — 0.16 cos 86 — 0.06 cos 948 — 
0.01 cos 106 — 0.03 cos 118 + 0.01 cos 126) ( 


y= (1 


d) 

We see that the location of the axis through the c.g. 
of the elastic weights is particularly simple when the 
Fourier series method is used. First, the series ex- 
is found referred to any 
Then, by merely dropping the initial 
the axis 


pression for y’ convenient 
horizontal axis. 
constant, the expression 
through the c.g. of the elastic weights is obtained. 

We have now satisfied the conditions required for 
the use of the elastic weight method. Returning to 
Fig. 2b, let X,,, X,, and X, denote the indeterminate 
respec- 


for y referred to 


moment, horizontal force, and vertical force, 
tively, assumed to be acting at the end of the rigid 
arm L. These quantities are taken as positive in the 
directions indicated. 

At any point P(x, y) of the ring, the bending moment 
M is given by 


M = Myot+Xnt+Xxt+ Xay (8) 


where /, is the moment that would exist because of 
the external loads on the ring if XY,,, X,, and X, were 
It may be shown that 


= S My ds —_ S Mox ds 
ore: Sods ’ oe S x ds 


= S Moy ds Sy? ds (9) 


all zero. 


RIN 
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These are the basic equations of the elastic weight 
When MM, is not symmetrical about the 
vertical axis, the integrations must be 
completely around the ring. If, as in the 
case, there is vertical symmetry of MJ, the integrations 
need be carried out over but one-half of the ring, and 
X, will always be found to be zero. 


method. 
carried out 
present 


SOLUTION OF SAMPLE PROBLEM 


For the problem of Fig. 2b, we define 7%)’ as the 
M, at any section between a and } and JM,” 


M, at any section between 6 and c. 


value of 
as the value of 
Thus, 


M,’ = 0, M,” = (P/2 


Ax — €) 


where e is the x-coordinate of the point of application 
of the load and is found to be 6.80. 
Now 


*0 ot) 70 
Aun = / Mp ds / ds = / My ds/c 


where 
*O °*2c/3 *0 
/ M, ds = i M,’ ds + I M,”" ds 
, ; 20/3 
Pc >) Pe 0 
= x dé — ds 
2m J2x/3 2 2¢/3 


Integrating each of these terms separately 


Pc *0 
y x dé = 
2 J2x/3 


by term-by-term integration of Eq. | 


Pe v 2Pec 2P(6.80)c 
— ds = — = : = 2:27Pc 
2 SJ2c/3 6 6 


Then, 


—1.75Pc 


2), and 


Xm = (—1.75 + 2.27)Pc/c = 0.52P 


Because of symmetry, X, = 0. 


DETERMINATION OF 


--f{ Muy ds fr 


0 3 
fu Moy ds = Lf” My'y ds + c My"y ds 
| Pc 1d Pec o 
o= f xy dé — / y dé 
2 J2x/3 Qe J22/3 


The product xy, in which x and y are represented by 


series expressions, requires investigation. 
y = a,cos@ + a2cos 206+... 
b, sin @ + bd: sin 26+... 


Since 


+ diz cos 126{ 


10) 
+ by sin 1104 ‘= 


x 


obviously the product will be of the form 


12 11 


> Dd ancos mé-b, sin nd 


m=l1n=1 


xy = 
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We now recall the trigonometric identity 
cos mé-sin n6 = (1/2){sin(z + m)é@ + sin (7—m)é 


Using this identity and neglecting harmonics higher 
than those originally present, we obtain 
xy = (1/2)(A, sin @ + Asin 26+...+4+ 

Ky sin 126) (11) 
The values of the constants A, in the above expression 
are functions of the coefficients of Eq. (10) and are 
given by 

12 


K,= > a(bi-; + 5:4 (i = 1,2,3,..., 12) 


l 


subject to the definitions b_, = —b,, 6b) = 0, and 


b, = Oforn> 12. For example, 


Ky = a,\(O + be) + ao( —b; + b3) + a3(—b; + bs) + 
a4 — b, + bsg) +...¢ Qj» — by, + 0) 


Ky = ay(b, + bs) + a(O + bs) + a3(—b, + 55) + 
Gai — b., + bg) : Pe — Qy2\ — bin T 0) 


ete. Egy. (11) for xy may easily be integrated, giving 


*0 
/ xy dé 


_- (Ki + Ao + Ay + Ay + Az + Kg + Kyo + Aun) 


where Ay, Ay, Ay, and Ay have dropped out. Thus 
these constants need not be computed. Using the 


Fourier coeflicients previously determined for this ring, 
the above expression may be evaluated. 


ala be re. 
xy dd = - (78.63) = —9.39P« 


The term 


*O 
/ vy dé 


presents no difficulties; term-by-term integration of 


Eq. (7) may be used, giving 


Pe “0 
= vy dé = 9.00P« 
2a J 2s : 
“0 
/ Muy ds = —0.39Pc 
The denominator of the expression for X), is 


0 c 70 
/ yds / y? dé 
B © Js 
( 0 ( “0 Ae 12 
vy? dé dm Cos mé-a, cos né } dé 
wn WSJ m lp I 


Thus, 
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For m unequal to n the value of this integral is zero. 


For m equal to n, 


“| . 12 


*0 
Cc c ( 
/ y? dé = > a,? cos? n6 dé = 7. és 
Z , a > & 
wT Jax us ~ n 1 on l 


Substituting values of a, from Eq. (7) 
70 
/ yids = —47.06: 


X, = — (—0.39Pc/—47.06c) = —P/121 


Therefore, 


The final expressions for the bending moment are 


from a to b, 
M = P(0.52 — (y/121)] 


from 6 toc, 


The normal and shear forces at any section may be 
found in terms of P and X, by the usual methods, 
taking into account the angle the section makes with 
the horizontal. 

By substituting appropriate values for the Fourier 
coefficients a, and 0,, this solution may be applied to 
any ring with one axis of symmetry loaded as indicated 
in Fig. 2a, regardless of size or shape. The values to 
be used would, of course, be found from a Fourier 


analysis of the particular ring. 


ACCURACY OF THE METHOD 


As a check on its accuracy this method was applied 
to a true elliptical ring subjected to various conditions 
of loading. The solution of such rings has previously 
been worked out by Burke! through the use of elliptic 
integrals. The comparison of the results of the two 
methods of analysis for two specific loadings chosen 
at random is given by Table 2. It will be noted that 
the differences between the results of the two methods 
are small. There is no convenient way to check the 
accuracy of the Fourier series method for nonelliptical 
rings, but the close agreement with Burke’s method 
for elliptical rings indicates that the accuracy should 
be high. 

Also in Table 2 is illustrated the effect of using a 
reduced number of terms in the series. The extremely 
small magnitudes of the coefficients of the higher har- 
monics suggested the idea of dropping them entirely. 
For the true ellipse, the complete series contained only 
six terms because of the presence of both horizontal 
and vertical symmetry. When the last four terms of 
each series were dropped off, the effect upon the results 
was negligible. However, the computational effort 
required was greatly reduced, particularly when one 
series was to be multiplied by the other. 





ti 
be 
ne 
m 


al 


al 


th 
lo} 
ax 
th 
of 


ac 
ce 
th 
dit 


me 


Zero, 


are: 


y be 
10ds, 
with 


urier 
d to 
ated 
2s to 
urier 


lied 
ions 
usly 
iptic 
two 
osen 
that 
10ds 

the 
‘ical 
hod 
yuld 





ANALYSIS OF OVAL 


TABLE 2 


RINGS 


BY FOURIER SERIES 5&9 


9 


Comparison of Results of Different Analyses of Elliptical Rings (Major Axis Assumed Vertical 


Type Vertical forces 
of P Upward at s iD 
Loading Vertical Forces: P/2 Downward at s = 3c/4 
P Upward at s = 0 P/2 Downward at 5e/4 
Method P Downward at s = ¢ Horizontal Forces: 
of (Ends of Major Axis ( to the right at c/3 
Analysis No Horizontal Forces Q to the left at s = 5¢/3 
Method of Burke (elliptic integrals X, = XY, =0 \ 0) 
Xm = —1.76P X, 0.046P + 0.760 
Xm = —(5.16Q + 1.97P 
Fourier series, six terms X, = YX, =0 yy, = 0 
Xm = —1.77P Y, = 0.046P + 0.760 
£ —(5.180 + 1.98P 
Fourier series, two terms X; X, = 0 xX, = 0 
Xs —76P X, = 0.047? + 0.76Q 
Xe —(5.18Q + 1.98P 


APPLICATION TO FUSELAGE RING 


In airplane structures the loads applied to a fuselage 
ring are not resisted by concentrated forces as in the 
preceding example problem. Instead, the ring is held 
in equilibrium under load by the combined action of a 
series of concentrated forces supplied by the stringers 
plus a distributed reaction due to shearing forces set 
up in the skin. To simplify the analysis, the assump- 
tion is often made that the effect of the stringers may 
be accounted for by assuming an effective skin thick- 
ness greater than the actual thickness. This approxi- 
mation will be used here. 


EXAMPLE 


The problem is illustrated by Fig. 3a. For simplicity, 
an elliptical ring is used rather than a general oval ring, 
although the procedure is identical in both cases. 

The concentrated load P acting on the elliptical ring 
is opposed by shearing forces in the skin of effective 
thickness ¢. These shearing forces are a result of the 
longitudinal bending of the fuselage about a horizontal 
axis and are distributed around the circumference of 
the ring in some as yet unknown manner. The moment 
of inertia of the ring cross section is assumed constant. 
The effect of assuming that the skin shearing forces 
act at the neutral axis of the ring instead of at the 
centerline of the skin is regarded as negligible, since 
the depth of the ring is small as compared to its other 
dimensions. 

The are length s along the neutral axis of the ring is 
measured clockwise from the top, and the semicircum- 
ference is denoted by c. 

Following the method previously described, Fourier 
analyses may be carried out for the figure formed by 
the neutral axis of the ring. Since the elliptical ring 
has two axes of symmetry, all even terms drop out 


leaving 
y = a, cos 6 + a3 cos 30 + a; cos 50 + ... + ay, cos 116] 


= b, sin @ + b; sin 36 + bs sin 50 +... + d, sin 110) 
(12) 


= 


where 6 = ms/c and a, and 6, are the coefficients de 
termined from the Fourier analyses. 

The elastic weights method will be used with notation 
as indicated by Fig. 3b. The initial equation of this 


method is 


M=M+Xn+Xx 


+ Xny 
The 


first step in solving the problem is to determine an 


in which the quantities are as previously defined. 


expression for J/o, the bending moment that would 
exist at any section of the ring if Y,,, X,, and X, were 
all zero. 

In dealing with thin sections it is convenient to work 
with the shear flow g, defined as the shearing stress 
times the section thickness. For the case of longi- 
tudinal bending of the fuselage about a horizontal axis, 
the expression for g-may be found from any textbook 
of aeronautical structures to be (this shear flow de- 
velops only if the cross section of the fuselage at the 
location of the ring does not distort either in the plane 
of the cross section or perpendicularly to it) 


V : 
g = j / yvydA 


in which g is the shear flow in a fuselage cross section 
for which the vertical shear is |’ and at a point whose 


The 


distance from the bending axis is denoted by e. 











Fic. 3. 


Problem of elliptical ring with shear forces 
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distance to the outer extremity of the section from the 
bending axis is f, and J is the moment of inertia of the 
entire fuselage cross section about the bending axis. 
For this problem V = F and dA = tds = (tc/x) dé. 


| © te 
g= y— do 
I 0 sy 


where ¢ denotes the particular value of 6 at the point 
Now, 


Thus, 


in question. 


™ 11 


P ; tc 
[= / ydA = ff yas = — =, a," 
0 2 n =1,3,5 


as was found in the previous example, and 


7 


? , a3 . 
| y dé = =asing + = sin 3¢ + 
0 vo 
(a5/5) sin 5@ +... + (€y,/11) sin 11d 


Therefore, 
oP 


7“ T ) a," 


- a. 
(a sin @ + ~ sin d@ + ...+ 
o 
M1. , 
_- Ll¢ (13) 


Now consider Fig. 3c. Point A (x4, va, S = S4) iS any 
point on the ring at which we wish to determine the 
moment Wp due to the shear forces in the upper portion 
of the ring. If we can find an expression for the ele- 
mentary moment produced at A by the skin shear 
force at any particular point B (xz, yp, s = sz), then 
an expression for \/) at A may be obtained by letting 
the location of the point B vary from A to the free 
outer end of the half ring and summing the elementary 
moments by integration. 

If we define 8 = ms4/c and @ = mSz/c, then x4 and 
ya are functions of 8 and xg and yz are functions of ¢. 


COMPONENTS OF SKIN SHEAR FORCE 


To find the moment of the elementary shear force 
q ds at B about A, it is convenient to resolve this force 
into horizontal and vertical components, g,ds and 
g, ds, respectively. From Fig. 3c, with the horizontal 
component positive toward the left and the vertical 
component positive upward, 


grds = qdscosx, g,ds = qdssinX 


f F Py, E 23 » B 
) =— -OS 2 
E Grva AD P x re 2 cos 


The integration of the third term is similar, giving 
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The expressions for sine \ and cosine \ may be ob- 


tained from Eq. (12) by differentiation—that is, 


sind = —dy/ds = 
(r/c)(a, sin @ + 3a; sin 86 +... + 1Llay sin 11¢ 


cos \ = dx/ds = 
(1r/c)(b; cos @ + 363 cos 36 +... + 11d cos 11¢ 


Since g is also represented in series form, the multipli- 
cation may be carried out and simplified by trigono- 


metric identities. Disregarding higher harmonics, 


r . ‘ 
q: = (Be sin 26 + By sin 4¢ + 
Cc Zz, an” 


. + By sin 129) 


l4 
P 
= — (A.cos 2 + Aycos 4¢ + 
ly c>oa,? ? 
.-+Apcos 12¢) 
where 
a a 
A, = z: = [4 + Ja; — (1 — ja; 
1,3,5 a 
11 a 
B,; = > [(@ — fb_5 — G+ is 
) i a ae J 
(¢ = 2, 4, 6,..., 12) subject to the definitions a_, = a,, 
b_, = —),, a, = bd, = Otorn>11. It is to be noted 


that A; and B; are constants dependent only upon the 
shape of the given ring. 


DETERMINATION OF JJ, 


Having now expressions for the horizontal and ver- 
tical components of the shear force, the value of VW) at 
the point A is given by 

Cc "8B 
M, = r [gr(va — YB) + Gy(X4 — Xp) ldo 
T 0 


or 


¢ of) 3 
M, = / G2Va — I G2¥B + 
Tv 0 0 
on: } *8 
/ QyXa — I ave] do (15 
0 0 


For any particular point A the value of x, and yy, 
are known fixed quantities. Performing the above 


indicated integration gives, for the first term, 


F By B. . By By 
, eos 4b +... +7, cos 128 — + +...+ 


2 } 12 


f , Px, (2 in 26 + A, . 18 + “ Ay . e+ 6D :) 
= = F sin 2 sin ¢ ser S yA ‘ n- 
GyXa Ab i r 5 4 12 sin a 


The second and fourth terms are somewhat more complicated, since xg and yg are expressed by Eqs. (14) with 


6 replaced by ¢, the specific value of @ at B. 


However, by multiplying out and simplifying, these two terms 
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may be integrated. We find 


~B P [ ” ‘ C; : Cu os C; . Cu 
: Gr¥p dd = — 2 a,” pm et 3 cos36+...4+ T cos 118 —{Q+ 3 +... . 


> 


- I ee Du D; Du\ 
/ qyxp dpe = — - -| D, cos 8B + — cos38+...+ cos 118+({Di+— +...+ + 
0 2c Do an 3 11 3 11 


oe 


bs b; b bs b, b 
2 > a,” (0, cos 8 + ~ cos 38 + ~cos58+...+ 7 cos 118) —2 > a,” (0, = = F268 4 ") 
» >) e 


” 


C, and D; are constants given by 


11 11 


C; = p> a,;(B;_; + Bi +;), D, = + b(A;—; — Aiz;) 
j 1,3,5 ) 1,3,5 os 
where we define B_, = —B,, A_, = A,, and A, = B, = O for x = 0 and 2 > 11. 
Combining the various terms, the expression for WV, for the section A at which @ = 8 becomes 
a ' C; + D; Ci +D 
M = ,| (Gi + Dd(e0s 8 = 1) + (oe 36 — 1) +...4+ — “(cos 118 — 1) | + 
2 >a,” 3 11 
P bs by 
b, (cos 8 — 1) + — (cos8B — 1) +... 4+ (cos 118 —1)] — 
7 3 1] 
Py B, - B, By _ 
Ya] 2 (cos 2B — 1) + (cos48 —1)+...4+ 19 (cos 128 — 1) 


a 
Px { A 2 A 4 A 12 
- sin 28 sin +8 awe “sin 128+ 6S-a,? (16) 
Fan \e" + 4 + + 12 sl 2, Gn ) 


This represents the bending moment at A in the determinate structure due to the distributed skin reactions. 
It should be recognized that x4, y4, and 6 are constant for any given section A but change in value as the loca- 
The reason for leaving x4 and y, in this form rather than using their series ex- 


tion of the section is changed. 
Ordinarily, the coordinates of the point at which the moment 


pressions is to simplify the computation of Wp. 
is desired will be known and may be substituted directly in Eq. (16). 
Now that an expression for /) has been determined, the quantities X,,, X,, and X, may be found. It is im- 


mediately obvious that X, must be zero due to symmetry. 


EVALUATION OF X,, AND X,, 


According to the elastic weights method, 


¢ ‘ rT r l x 
Xn = — (f Moy ds f ds) =— (f M, dp f as) =— | M, dg 
0 0 0 0 T 0 


Carrying out the integration, 
P 


~ 27 > Ss. 


ae ( f Moy ds | f y? is) =- ( f Moy dB | | y? is) 
0 / 0 0 / 0 


Upon multiplying out and integrating, this becomes 


: a ; C: D: ; Cy + D ; 
X) odin 9 co ia EG + D, — Ei) + a; ( 4 _ Fs) +...t+@n ( == —_ Fn) + dF; + b3F3 + 
2r( >> a,2)? 3 11 


+ bus] 4 | ett... todd eeo es +o 2(S4S4 + Ge) 
oe = 4,0 303 sees 1 eee =a eee 
uf u x >. a," 1,01 a; Qy91u 9 12 9 4 12 


C3 + D; Cu ss Dy - 
+ + (17) 


Ke (c + Di + er 
3 1] 


Now, 


= 4 2 
(18) 
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The constants are evaluated as follows: 


11 
i = > 6} ({Ai4,/(i+7)] — [A 
1.3.5 
1 
G, = 2, a,(b;~-; + 5445) 
) 1,3, 
subject to B_, = —B,, A_, = An, b-n = —by, and B,, 


The bending moment at any section of the ring may 


now be found from 
M = M) + S & 4 XnVo 


where MM, for the particular section is found from 
Eq. (16). X,, and X, are the values found from Eqs. 
(17) and (18), and yo is the ordinate of the point in 
question. 

The expressions for the normal and shear forces are 


N = (X, + gr ds) cos X + [ g, ds sin X 
0 J/0 


V = (X,+ / gz ds) sin X + / gy, ds cos X 
0 0 


in which sy is the value of the are length at the section 
under consideration and the other quantities are as 
previously defined. 

The several constants a;, 6;, A;, Bi, . 
tions only of the shape of the ring and are calculated 


. . G, are func- 


without reference to any particular section. Once they 
have been computed, the bending moment at a number 
of sections may be readily found. 

Because of the representation of the shape of the 
ring by Fourier series, the foregoing analysis is valid 
not only for a true elliptical ring but also for any ring 
symmetrical about two axes. The effect of different 
shapes of rings is merely to change the values of the 
constants used in the general expressions derived above. 


EFFECT OF HIGHER TERMS 


As before, the higher terms of the series for x and y 
may be neglected with but small effect upon the final 
result. In a particular case, when all six terms of the 
series were used, the results were found to be 


X, = —0.142P, Xm = 1.39P 


while, when only the first three terms of the series were 
used, it was found that 


X, = —0.145P, ae 1.41P 


The error is about 2 per cent, while the amount of 
computation saved is significant. It is particularly 
desirable to neglect many of the higher terms in the 


case of rings having but one axis of symmetry when 


aj} (Bi_;/(i — j)] + [B 


(¢+7)}} (¢ = 1, 3, 3, re ht 

( ws 1) ei “a 

1 DI (¢ = 1, 3, 5, me 
(2 = 2, 4, 6, , 12) 


=A. = 6b, = Ofer sn = 0 and «> 11. 


the series for x and y each contain both odd and even 
terms. 


CONCLUSIONS 


The method here presented is simple, in itself, and 
it requires no extensive mathematical background. It 
is applicable to any ring with at least one axis of sym- 
metry, whether or not the loading is symmetrical. For 
nonsymmetrical loadings, the entire ring, rather than 
only one-half, must be considered and the quantity 
X, must be evaluated. 

It is recognized that the use of Fourier series pre- 
sents certain disadvantages. The most important of 
these from the practical point of view is that the 
process of deriving original solutions is somewhat 
cumbersome, particularly when distributed skin re- 
actions are considered. However, once a solution has 
been arrived at for a given oval ring and condition of 
loading, the bending moment at any number of sections 
may be quickly determined. Also, the solution may be 
applied to any other ring with similar loading merely 
by varying the constants. As was pointed out, the 
amount of computation required may be considerably 
reduced by using fewer terms in the original series with 
but little sacrifice of accuracy. 

This method of analysis was originally intended to 
apply to aircraft fuselage rings. It may, however, be 
advantageously applied to any ring having one axis 
of symmetry. Since commonly only concentrated 
loads will be involved in other than aeronautical 
applications, the solutions will be simple. 
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Two-Dimensional Airfoils at Moderate 
Hypersonic Velocities 


WILLIAM H. 


Consolidated Vultee 


(1) SumMMarRyY 


Expressions for lift, drag, and pitching moment coefficients 
are derived for a large number of airfoil sections in two-dimen- 
sional hypersonic flow. These expressions, differing from those 
obtained using linearized, second-order, and third-order theories, 
can be used throughout a moderate range of hypersonic Mach 
Numbers. The fundamental parameter of these formulas is the 
two-dimensional hypersonic flow similarity parameter, K = 
Mé. To obtain these expressions, the shock-wave and expan- 
sion-wave equations for pressure ratios across such waves are 
expanded in powers of A. Following this, it is shown that the 
first three terms so obtained for both wave types are practically 
identical for values of A less than about 1.0. Further, it is 
shown that they closely approach the exact pressure-ratio ex- 
pressions. The application of the expressions thus derived to 
yarious airfoil sections is then demonstrated, and a table of air- 
foil coefficients is presented. It is shown that the airfoil formu- 
las can be applied to finite wings with reasonable accuracy for 


design estimations 


(II) SymMBo_s AND NOMENCLATURE 


A = plan-form area 
A.R = aspect ratio 
Cp = drag/g/ = drag coefficient per unit span 
C, = pressure coefficient 
c. = F./4l = chordwise force coefficient per unit span 
C, = F,/cl = normal force coefficient per unit span 
Co = Mo/<l/? = pitching moment coefficient per unit 
span (positive in stall) 
f = f(t, 8) = thickness distribution function [see Eqs. 
(3) and (4)] 
F,, F, = chordwise and normal force components 
H = stagnation pressure 
| = indices [see Eq. (49) |] 
K = M6, hypersonic similarity parameter 
l = chord length 
M = Mach Number 
Mo = pitching moment, positive in stall 
p = static pressure 
q = dynamic pressure 
= area ratio (see Fig. 8) 
S = span 
t = thickness 
é., = chordwise location of pitching moment center 
9 = body coordinates (see Fig. 1) 
Greek Sym bols 
a = angle of attack 
B = a/b 
Y = ratio of specific heats 
6 = ?¢/l, thickness ratio 
n = 2y/t, dimensionless coordinate (see Fig. 1) 
\ = shock-wave angle 
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— = 2x//, dimensionless coordinate (see Fig. 1) 

w = flow deflection angle, positive ‘‘into’’ a corner, negative 
“around” a corner 

Subscripts 

Q = free stream when used with Mach umber; origin 
when used with pitching moment 

1 = local 

/ = lower surface 

u% = upper surface 


The symbol O|[ | denotes “‘order of’ term within brackets 


(III) Discussion 


y HAS BEEN DEMONSTRATED that long-range missiles 
require hypersonict Mach Numbers for. efficient 
operation.’ ? To date, considerable work has been de- 
voted to the aerothermodynamics of hypersonic flight 
in the upper atmosphere. In this region the molecular 
mean free paths are comparable to or exceed the sig- 
nificant dimensions of a missile and its external flow field, 
and the concepts of the kinetic theory of gases appar- 
ently apply.*~* This paper, however, is concerned 
with the aerodynamics of hypersonic flight at alti- 
tudes wherein the laws of continuum flow theory apply. 

Ivey, Klunker, and Bowen" proposed a method of 
treating the continuum flow about simple bodies at ex- 
treme hypersonic Mach Numbers. This method, 
based on Newtonian flow concepts, might be expected 
to apply at Mach Numbers of the order of 10 to 15 
(say) or higher if the effects of viscosity and molecular 
vibration and dissociation in the flow field do not 
modify the flow behavior appreciably. Since linear- 
ized theory is of dubious value above Mach Numbers 
of about 3, the missile designer is left with a void in the 
continuum flow region between Mach Numbers of 
about 3 and 10. Here he must either apply lengthy 
numerical methods, such as the method of character- 
istics, or make use of extrapolations to obtain pressure 
distributions and force and moment coefficients. It 
is the purpose of this paper to present a method of ob- 
taining such aerodynamic coefficients for two-dimen- 
sional airfoils in this Mach Number region. The 
method appears to be applicable for airfoils of sensible 
thicknesses and angles of attack at Mach Numbers 
between about 3.0 and 12. The coefficients obtained 
using this method compare well with those obtained 


+ For the purposes of identification it will be assumed that 
hypersonic Mach Numbers refer to Mach Numbers of 3 and 


above. This definition is admittedly arbitrary. 
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within dimensionless coordinate system (£, 7). 
using the shock-wave and expansion-wave equations in 
the customary piecewise manner. 

In 1946 Tsien, using the equation of motion for two- 
dimensional compressible flow about slender and thin 
bodies, found a similarity parameter relating hyper- 
sonic flows about families of similar bodies.'! This 
similarity parameter, when properly used, relates the 
pressure distribution and force and moment coefficients 
due to pressure stresses on a given body at a given hyper- 
sonic Mach Number to the pressure distribution and 
force and moment coefficients for bodies of the same 
thickness distribution at different hypersonic Mach 
Numbers related to the given Mach Number by the 
similarity rule. Indeed, if 6 is defined as being the 
thickness ratio of the airfoil, if 


| 

| 
lA 
lA 


-= 
s 


lA 


n = 2y/t, —] <1 (2) 
are defined as being dimensionless coordinates of the 


airfoil (Fig. 1), and if 


fi(& B) = (dn/dé), 8 (3) 
filé, B) = — (dn/dé), + B (4) 
where 
B= a/é (5) 
then, as will be demonstrated subsequently, 
C/O? = PCR, B;-£) (6) 
therefore, 
C,/6 = Y(K, B) (7) 
C./8 = ACK, Bp) (8) 


C,/6 = 0 (K, B) (9) 
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where 


is the hypersonic similarity parameter. 


Linnell has utilized this similarity parameter to de- 
rive and present expressions for pressure, force, and 
moment coefficients for two-dimensional airfoils in 
hypersonic flows.'? His expressions for pressure co- 
efficient, however, were different for flow fields behind 
compression waves and expansion waves and, as such, 
precluded the obtaining of force and moment coefficients 
for airfoil sections in closed form. Nevertheless, 
Linnell’s method represented some saving in labor 
over the application of shock-wave and expansion-wave 
formulas to two-dimensional airfoils when A was of the 
order of 1 or less. Smelt showed that, for hypersonic 
flows of low values of A, one expression for the pressure 
ratios across shock waves and expansion waves could 
be derived which approximated the full Rankine- 
Hugoniot and Prandtl-Meyer pressure ratios accu- 
rately.'* In obtaining these expressions, Smelt used a 
slightly different similarity parameter not intrinsically 
different from Tsien’s. In essence, this removed the 
dependency of the pressure ratio for any point on the 
surface of a two-dimensional airfoil on the number and 
type of wave fronts between the point and the undis- 
turbed free stream. Airfoil coefficients in closed form 
could then be obtained. 


The present paper derives expansions in powers of 
K from the Rankine-Hugoniot and Prandtl-Meyer 
pressure coefficient ratios. Because Tsien’s A is the 
parameter of the expansion, the expressions obtained 
differ from Smelt’s expressions. The two expressions 
derived in this manner are shown to be almost identical 
and closely approximated by one three-term equation 
for K < 1.0. Using the three-term approximation, 
general and particular expressions for the aerody- 
namic coefficients of two-dimensional airfoils are ob- 
tained in closed form. 


(IV) PRESENTATION OF THE THEORY 


(a) The Expansion-Wave Equation for Pressure Coefficient 


The Prandtl-Meyer equation relating Mach Numbers 
across a supersonic isentropic expansion through an 
angle w can be written [see reference 14, making use of 


Eq. (13.10b) | as 
Y= ‘) 
(M,? — 1)") 
(; + 1 : 


(? + +) ‘ 
—w = tan 
y— |] 


y¥-—1 2 “ 1/9 
tan~! aS (Mo? — 1) + 


{tan—! [M,? — 1] — tan? [M,2 — 1]”*} (11) 
The symbols in Eq. (11) are illustrated in Fig. 2a. For 
Mo2[(y — 1)/(y + 1)] > 1, Mi > Mo, Eq. (11) can be 
approximated by 
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TWO-DIMENSIONAL AIRFOILS 


— 240-37)" 
Oe a= th a 


of 4 (= pe =) it (12) 
3(¥ —_ 1)? M,'M,? f 


Now, for isentropic flow, we can write [Eq. (45) refer- 


ence 15 | 


pi _ {Mo? 1 + [2/(v — 1) Mo?}\”' 


by M21 + [2/(y¥ — DMEIS (13) 


1 


and, for My > [2/(y — 1)]”’, 


2y/y-1 y¥/y-1 
pr _ (**) fi 4 of 4 } (14) 
pb \M ( (y — 12MM dS 


or, neglecting terms of the second order in the brackets 
and combining the first term of Eq. (12) with Eq. (14) 
yields 


; Po yM,? 
‘ 2y/¥ l 
2 $14 (7—'\un ~i> «5 
7M," ' 2 f 


Further expansion yields, for |[(y — 1)/2].Mow| < 1, 
I . 


2w l + | 
en A a eB ee 


Mo 2 6 
( 1) (38 — y) Myo 
rr MiB al Ln (16) 
1S 
Now, from Fig. 1, 
w = f(t, B)b 
therefore, from Eq. (16), 
we ) 
Cp _ eee mes Ley 
& K- 2 
“= . g 2 = "24 
(y+ DK sg (y+ 1) (3 y)K?f (17) 
6 41S 


where, from the limit of convergence for Eq. (16), if 


imi <5 (18) 


The expression corresponding to Eq. (17) for pressure 
coefficient for surfaces behind oblique shock waves will 


now be derived. 


(b) The Shock-Wave Equation for Pressure Coefficient 


The Rankine-Hugoniot equation for pressure ratio 
across an oblique shock wave can be written {see refer- 
ence 14, Eq. (4.29)] as 

b 2 5 alia) — | 
Pr Y_ M, sin? x» — 2 (19) 
fo ¥ 3 * Tt 1 


Further [from reference 14 Eq. (4.27) ], 


: / 1 sin A sin w ] 
sin? \ — Yt = ee = () (20) 
2 cos (A — w) M,? 
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2b SHOCK WAVE 
Fic. 2. Definition of angles used in shock- and expansion-wave 


equations. 


The symbols appearing in Eqs. (19) and (20) are illus- 
trated in Fig. 2b. Combination of Eqs. (19) and (20 


yields 
= 4 j } 
>= am? A. — (21 
y¥+ 1 M.?(y + 1) 
Now, 
sinw fo 4 %Q _ wy +0[ * } 
cos(A—w) | 6 2 120 ]/ 


(29 


Combining Eqs. (20) and (22) and solving for roots of 


sin A yields 


: yr i 4 wo (A — w)* FE ) 
sin A = w<1l— + O ~+ 
4 ( 6° 2 120}4 


lyjfy+i1 /- wv]? 40 = 
w}| |1—O 
TC lel tar | 


where the positive sign is used since sin \ > 0. For 
w<ar/2andd < 1, 


6 is 
sind = 7+ 14h 4/1 + © | t (a4 
4 { (y + 1)2Mow?} 


bo 
wo 


If Eqs. (21) and (24) are combined, there results 


: y¥+ 1 y¥+1\? 2 it" a 
CG, = = ot + ot [( r y'+(52)] (25 


This equation agrees with expressions presented by 
Linnell'? and Ivey and Cline."* Again, if, from Fig. 1. 
w = f(&, B)éand if K = Md, then for 
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pansions using exact, three-term approximation, second order, 
and linear theories 


(y+ DAF (vy + I) Mw 


| | 
Eq. (25) can be expanded into 
y+ 1 (y + 1)7K 


== f+ f2 + ——f* + 
K P 16 


a 4 “3/5 
fee (26) 
1,024 


where, from the limit of convergence for Eq. (26), if 


font 
wis * 
to 


in?) 


y = 1-4, 
SK) < 1.66 (27) 
and, as can be shown for w ~ 0 and (A — w)? < 0.10, 


using Eqs. (20) and (22), 


Eq. (26) completes the derivation of the basic expan- 
sion and shock-wave equations for pressure coefficient. 


(c) The Equation for Isentropic Low Pressure Coefficients 


It is apparent that Eqs. (17) and (26) are identical 
up to the third term. In fact, if it is assumed that 
Eq. (17) holds for both shock- and expansion-wave proc- 
esses, the error introduced by the assumption is 10 per 
cent or less of the third term of Eq. (26) in the interval 
of convergence of Eq. (26) for y = 1.4. Essentially, 
this implies that compressions at wedge angles satisfy- 
ing inequalities (27) and (28) are isentropic for all 
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practical purposes. Laitone” and others have pointed 
out that oblique shocks with small flow deflection angles 
and normal shocks at low supersonic Mach Numbers 
are essentially isentropic as is illustrated by the small 
loss in total pressure across such waves. The total 
pressure loss across an oblique shock wave of incli- 
nation angle \ is identical to the total pressure loss 
through a normal shock with free-stream Mach Num- 
ber \/) sin X. Therefore, the maximum total pressure 
loss for the shock waves concerned within this paper 
can be calculated using Eq. (24) to obtain maximum 
Mo sin \ by substituting Af = \/jw= 1. For this case 
My sin X = 1.766, and, from Table III of reference 15, 
H,/Hy = 0.829. 


entropy when compared to the entropy increase pos- 


This is a relatively small increase in 


sible with larger flow deflections even at the lower 
range of Mach Numbers in the range of inequality (28), 

Because both expansions and compressions are essen- 
tially isentropic within the interval 


fK| < 1.0, My > 3.19 


the three-term approximation derived from Eq. (17) 
will be used to represent both compression and ex- 


pansion-wave pressure coefficients. That is, 


ie 2,. 87! .& + TK, ad 
82 - K 7; 2 kil oa 6 f (29) 
Fig. 3 presents the pressure coefficients for the hyper- 
sonic flow over the upper and lower surfaces of a flat 
plate. For the case of the flat plate 6 = 0 and from 
Egs. (3) and (4), f6 = *w, where the negative sign ap- 
plies to the upper surface. As can be seen from Fig. 
3, the three-term approximation agrees very well with 


shock-expansion theory.'* Linear theory and second- 


order theory expressions are also presented in Fig. 3. 
The three-term approximation exhibits better agree- 
ment with the shock-expansion theory than does either 
of the other two theories. Third-order theory, while 
not plotted in Fig. 3, would agree substantially with 
the three-term approximation. 

Fig. 4 illustrates the region of validity of the three- 
term approximation in the /), w plane. The region 
of validity of the linear theory has been superimposed 
on this plot for comparison purposes. This somewhat 
arbitrary region was obtained from the linear theory 
bounding inequality 

Cle 2 re l 
 |{Mo? — 1]"3| ~~ Mo? 
or, being liberal with this inequality and letting A 
Mow, 
K| < 0.125 [1 — (1/M?)]o' 


As can be seen in Fig. 3, Eq. (29) has a relatively wide 
region of applicability when compared with the region 
of applicability of the linear theory shown on the same 
plot. It is interesting to note that the lower Mach 
Number limit of the three-term approximation is in the 
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TWO-DIMENSIONAL AIRFOILS 
vicinity of the practical upper Mach Number limit of 
the linear theory. 


(V) Two-DIMENSIONAL AIRFOIL COEFFICIENTS 


a) Derivation of Airfoil Coefficients 


Eq. (29) presents an expression for pressure coeffi- 
cient which is an explicit function of thickness ratio, 
§: chordwise dimensionless coordinate, £; normalized 
angle of attack, 8; hypersonic flow similarity param- 
eter, A; and the ratio of specific heats, y. Because 
it has been shown that this equation will adequately 
represent the pressure coefficient at any point (&, 7) on 
the surface of an airfoil without reference to the num- 
ber and type of wave fronts between the point and the 
free stream, it is possible to obtain expressions for the 
common airfoil coefficients in closed form for any section 
shape for which the locus of the surface can be ex- 
pressed in algebraic form. The resulting expressions, 
of course, can only be expected to be valid for values of 
My greater than 3.19 and for | fK| < 1.0. 
for chordwise force, drag, normal force, and pitching 


Equations 


moment coefficients of a large number of airfoil sections 
are presented in this section following the derivation of 
the general expressions suitable for any airfoil section. 

If an airfoil is situated in an x, y coordinate system 
such that the chord lies along the x axis and the mid- 
chord point is the origin, it can be shown that 


| “l/2 dy dy 
.. — oe ( < ) 2 ( - ) | , 
l / 1/2 | 7 dx !,, Co dx/, - 
, * . 9) 
Cc, = / (c,, = €,,) a (31) 
f 1/2 


Cn = C.+ al, (32) 
“1/2 
Cs, = P f (Cry — Cy)x dx (33) 
Co, = Cay + (X0/1)C, (34) 
for thin airfoils and small angles of attack. Referring 


to Fig. | and using the dimensionless coordinates &, 7 
defined by Eqs. (1) and (2), Eqs. (30)—(34) become 


? am fi 
C,/8? = (1/2) / (Padu + P fi) dé — (a if (35) 
oe 1 
_1 
C,/6 = (1/2) f (P, — P,) dé (36) 
—1 
7] 
Cp/5? = (1/2) / (Pufu + Pifi) dé (37) 
_] 
: 1 
C,,/6? = (1/4) f (P, — P,t/dé (38) 
—1 
Cy,/8? = (Cy/8*) + (&a/26*)C, (39) 


If we define 
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REGION OF SUBSONIC FLOW 





—_ - a .4 m 
Ia| 
Region of validity of three-term approximation com 
pared with region of validity of linear theory 


Fic. 4. 


Da 
C; = K (40 


then Eq. (29) becomes, 
P(t, B) = Af + Of + Gf? (41 


Further, defining a recurring integral /;; as 


ee) Ms ji =0,1 Beats 
ly = (5) J eras Vj=1234f (42 


Eqs. (36)-(38) can be written, using Eqs. (41) and 
(42), 


C,/ = » [Loj): — Uoy)ulC, (43) 
4 
(C,/6*) + (aC,/6*) = uo [(Los)u + (Los): J Cj = Cp/8? 
it (44) 
3 £, 


C,/2® = DY (Kyu — Ty; = Cy/8 — 


j=1 ” 26" 
(45) 


Once the functions /,(é, 8) and f,(&, 8) are expressed in 
algebraic form for any airfoil, the integrals (43) can be 
evaluated and the airfoil coefficients obtained from Eqs. 
(43), (44), and (45). The limit fK < 1.0 can be ex- 
pressed, using Eqs. (3) and (4), and solving for the 
limiting value of angle of attack, 


< (1.0/M) — 


Qmax.| > 
where Jf, > 3.19. These inequalities define the range 
of validity of the airfoil coefficients obtained using Eqs. 
(43), (44), and (45). 

Table | presents airfoil coefficients for some common 


(dy/dX) max (46) 


supersonic sections using these equations. Once values 
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TWO-DIMENSIONAL 


of the C; are determined by the thickness ratio, Mach 
Number, and ratio of specific heats, the coefficients are 
obtained as functions of 8 from Table 1. 


(b) Numerical Examples and Comparison with Theory 


Fig. 5 represents normal force coefficient versus angle 
of attack for a flat plate in the region of validity of 
inequalities (46). Also Newtonian 


flow’ !*"" result applicable for 1J) ~ © obtained from 


shown is the 


Eqs. (15), (25), and (31 


C, = (y + la? 


Since the expression for normal force coefficient ob- 
tained using the three-term approximation is valid up 
to Mach Numbers of about 10 or 12, Fig. 5 suggests 
that shock-expansion theory be used for Mach Numbers 
between about 10 and those Mach Numbers where the 
theory of reference 10 is adequate. 

Reference 19 presents experimentally determined 
values of normal force, pitching moment, and drag 
) per cent thick double-wedge rec- 
Figs. 


coefficient on a 
tangular plan-form wing at Mach Number 6.9. 
jand 7 present a comparison between the experiment- 
ally and theoretically determined values of normal force 
coeflicient, center of pressure location, and drag coeffi- 
cient, respectively. Inspection of Figs. 6 and 7 reveals 
that the present theory appears to agree with the experi- 
mental evidence even though the wing tested had an 
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Fic. 6. Lift coefficient and center of pressure location versus 
angle of attack for symmetrical double-wedge airfoil at 1/7 = 6.9, 
AR. = 1.0, and 6 = 0.05. , Three-term approximation; 
, A experiment (reference 19). 


AIRFOILS AT 


HYPERSONIC VELOCITIES 599 





.05 7 
| 





| | 











DRAG COEFFICIENT~Cp 




















10 





ANGLE OF ATTACK~a° 
Drag coefficient versus angle of attack for symmetrical 
1.0, and 6 = 0.05 
, experiment (reference 
19). The addition of a skin-friction drag increment equal to 
0.6028 brings theory and experiment into better agreement. 


Fic. 7. 
double-wedge airfoil at JJ = 6.9, A.R. = 
, Three-term approximation; 


aspect ratio of 1 and the present theory applies to in- 
finite aspect ratio wings. 

The agreement between the three-term approxima- 
tion for infinite aspect ratio wings with measurements 
on a wing of finite aspect ratio is not surprising consider- 
ing that the region of three-dimensional flow is con- 
fined to the region shaded on the plan form of Fig. 8. 
For a flat-plate rectangular plan-form airfoil, the ratio 
of the area of three-dimensional flow to the total plan- 


form area is, approximately, 


R = [M2 — 1])7~7“(A.R.) 


This ratio is plotted versus A.R.[1/,? — 1]! ? in Fig. 8 
for Mo greater than 3.19. As can be seen from Fig. 8, 
this ratio is less than 0.33 for plan forms of aspect ratios 
equal to or greater than 1.0. Since the lift in the three- 
dimensional flow is somewhere near one-half that indi- 
cated by two-dimensional theories such as that pre- 
sented in this paper, the error introduced by using two- 
dimensional theory is somewhat less than 33 per cent 
for aspect ratios less than or equal to 1. This consider- 
ation indicates that two-dimensional theory can be used 
for three-dimensional thin wings for engineering esti- 


mates with small error. 


(VI) CONCLUSIONS 


(1) An expression for pressure coefficient as a func- 
tion of the hypersonic similarity parameter and flow 
deflection angle has been derived suitable for use in 
the region My > 3.19 and | fA) < 1 with reasonable ac- 
curacy. 

(2) General and particular equations for lift, nor- 
mal force, drag, chord force, and pitching moment co- 
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region for rectangular flat plate, MW) > 3.19. 
efficients of two-dimensional airfoil sections are derived 
and presented in closed form for use when J/,) > 3.19 
and | fK| < 1.0: 

(3) It is shown that the airfoil section formulas can 
be used for finite wings in the region of applicability of 
the theory for design estimates. 

(4) The comparatively simple expression for the air- 
foil section coefficients can be used to determine pro- 
files optimum with respect to aerodynamic character- 
istics, bending strength, torsional strength, etc. 


REFERENCES 


1 Ivey, H. R., N., Jr., and Oborny, L. F., /ntro 
duction to the Problem of Rocket-Powered Aircraft Performance, 
N.A.C.A. T.N. No. 1401, December, 1947 


Bowen, E. 


SCIENCES-SEPTEMBER, 1952 


* Ley, W., The Upper Atmosphere 
ploitation, Aeronautical Engineering Review, Vol. 10, 


Its Exploration and EB, 
No. 2, p 
20, February, 1951 
> Zahm, A. F 
Institute, Vol. 217, pp. 153-166, 1934 
The Gas Kinetics of Very High 
translation of ZWB Forschungskericht M. 972, May, 
N.A.C.A. T.M. No. 1270, May, 1950 
‘Tsien, H. S., 
Gases, Journal of the Aeronautical Sciences, Vol. 13, No. 12, pp 
653-664, December, 1946. 


Superaerodynamics, 


Flight Speed 
1938 


+ Sanger, E., 


Superaerodynamics, Mechanics of 


5 Snow, R. M., Aerodynamics in a Highly Attenuated Atmo 
phere, U.S. Navy Publication RLA-4 

7 Ashley, H., Applications of the Theory of Free Molecule Floy 
to Aeronautics, Institute of the Aeronautical Sciences, Fairchild 
Fund Paper No. 164 

8 Stalder, J. R., and Jukoff, D., Heat Transfer to Bodies Trav 
ing at High Speed in the Upper Atmosphere, N.A.C.A. T.N. No 
1682, 1948. 

® Heineman, M., Theory of Drag in Highly Rarefied Gases, 
Communications on Pure and Applied Mathematics, Vol. 1, No 
3, pp. 259-273, September, 1948 

10 Ivey, H. R., Klunker, E. B., and Bowen, E. N., A Method 
for Determining the Aerodynamic Characteristics of Two and Thre 
Dimensional Shapes at Hypersonic Speeds, N.A.C.A. T.N. No 
1613, July, 1948 

1! Tsien, H. S., 
of Mathematics and Physics, Vol. XXV, No. 3, p 
1946. 

2 Linnell, R. D., 
Flows, Journal of the Aeronautical Sciences, Vol 
January, 1949. 

13 Smelt, R., Problems of Missiles at 
Ordnance Laboratory Report No. 1131, 


Similarity Laws of Hypersonic Flows, 


247, 


Journal 
October, 
Airfoils in Hypersoni 
16, No. 1, p. 22, 


Two-Dimensional 


Extreme Naval 
Symposium on Ordnance 


Speeds, 


Aeroballistics, June, 1949 

14 Liepman, H. W., and Puckett, A. E., /ntroduction to the Aer 
dynamics of a Compressible Fluid; John Wiley & Sons, Inc., New 
York, 1947 

16 Anon., Votes and Tables for Use in the Analysis of Superson 
Flow, N.A.C.A. T.N. No. 1428, December, 1947 

16 Ivey, H. R., and Cline, C. W., Effect of Heat-Capacity Lag 
on the Flow Through Oblique Shock Waves, N.A.C.A. T.N. No 
2196, October, 1950. 
7 Laitone, E. V., Exact and Approximate i 
Dimensional Oblique Shock Flow, Journal of the Aeronautical Sci- 
ences, Vol. 14, No. 1, p 

8 Ivey, H. R., Stickle, G. W., and Schuettler, A., Charts for 
Determining the Characteristics of Sharp Nosed Airfoils in Tu 
Dimensional Flow at Supersonic Speeds, N.A.C.A. T.N. No. 1148, 
January, 1947 

9 McLellan, C. H., Exploratory Wind-Tunnel Investigation 
Wings and Bodies at MM = 6.9, Journal of the Aeronautical Sci- 
18, No. 10, p. 641, October, 1951 


1 Solutions of Twe 


25, January, 1947 


ences, Vol 





Journal of the Frankliy 


Ra rened 





aoe. | 


In 
these 
stabil 
nves' 
Pres 
Meetir 
* Ae 
T Ac 
T Ac 


und Ex 


io. 2, p 
‘ranklin 


Speed 


15 


1938 


Rare fied 
12, pp 


| Atmos 


ule Flow 


‘airchild 
s Trav 
N. No 
1 Gases, 


1, No 


Method 
id Thre 
N. No 


Journal 
Yctober, 


personu 


lp. 2, 


Naval 


rdnance 


he Aero- 
c., New 


pe rsonu 


ily Lag 
N No 
of Two- 
ical Sci- 


tarts for 
in Two- 
o. 1143, 


ation of 
ical Sci- 





The Effects on Dynamic Lateral Stability and 
Controllability of Large Artificial Variations 
in the Rotary Stability Derivatives 


ROBERT O. SCHADE,* ORDWAY B. GATES, JRr.,f ann JAMES L. HASSELL, Jr. 
Langley Aeronautical Laboratory, N.A.C.A. 


SUMMARY 


This paper presents the results of a theroetical and experi- 
mental investigation of the effects of large artificial variations of 
several stability derivatives on dynamic lateral stability and con- 
trollability. The derivatives considered were the yawing 
moment-due-to-yawing velocity C,,, the rolling-moment-due- 
to-rolling velocity C),, the yawing-moment-due-to-rolling veloc- 
ity Cup the rolling-moment-due-to-yawing velocity C,, and 
the yawing-moment-due-to-rolling acceleration C,;. Although 
the results do not apply directly to airplanes or flight conditions 
ther than the two herein investigated, it is believed that the 
trends obtained might give an indication of the general effects 
of artificial stabilizing devices that cause changes in these deriva- 
tives. It was found that increases in oscillatory stability could 
obtained by changes in any of the derivatives, but in some 
cases these increases in oscillatory stability were accompanied by 
jetrimental effects on controllability and on the stability of the 
periodic modes. It appears from the results that the simul- 
taneous change in two of the derivatives might be more satis- 
factory than a change in a single derivative for improving the 
famping of the oscillation without impairing the other flight 


characteristics 


INTRODUCTION 


gel MUCH INTEREST has been shown in arti- 
ficial stabilization devices as a means for im- 
proving the damping of the lateral oscillation of air- 
planes. Such devices have been used on a number of 
airplanes in the form of yaw dampers. The yaw 
damper is essentially a gyroscope sensitive to the yaw- 
ing velocity with its output used to control a servo- 
mechanism that operates the rudder to damp the yaw- 
ing motion of the airplane. In an idealized system, 
such a device produces a variation in the stability de- 
ivative C,,,, the yawing moment due to yawing; similar 
levices can be used to vary the damping-in-roll deriv- 
ative Cy, the yawing-moment-due-to-rolling derivative 
(,, or the rolling-moment-due-to-yawing derivative 
C,,. 

In order to study the effects of large variations of 
these four rotary stability derivatives on the dynamic 
stability and control characteristics of airplanes, an 
investigation is being conducted at the Langley Aero- 
Presented at the Aerodynamics Session, Twentieth Annual 
Meeting, I.A.S., New York, January 28-31, February 1, 1952. 
* Aeronautical Engineer. 

t Aeronautical Research Scientist. 

t Aeronautical Research Scientist. 
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nautical Laboratory of the N.A.C.A. The effect of 
another parameter, C,,;, the yawing-moment-due-to- 
rolling acceleration, has also been studied, since it 
affects the stability of an airplane in a manner similar 
to the product-of-inertia parameter that is known to 
have a powerful effect on the lateral oscillation. 


Increasing the value of either of the damping deriv- 
atives C,,, or C;, increases the total damping of the air- 
plane, while changing the value of either of the cross 
derivatives C,, or C;, or the acceleration derivative 
C,,; primarily causes a redistribution of the natural 
damping of the system among the various modes of 
motion. Because of the cross coupling between the 
rolling and yawing motions in a lateral oscillation, it is 
possible to increase the damping of both of these mo- 
tions by the use of any one of the five derivatives under 
consideration. This point is emphasized because some 
groups working with automatic stabilization systems 
have the erroneous conception that an airplane must 
be stabilized independently about the rolling and yaw- 
ing axes. 


One phase of the investigation involved the correla- 
tion of dynamic lateral stability calculations for a 45° 
sweptback wing airplane (hereafter referred to as air- 
plane A) with experimental results obtained in the 
Langley free-flight tunnel from a '!/g-scale model of 
this airplane through a wide range of values of C,,, Cip, 
Cy», and C;,. These calculations and tests were made 
for a low-speed condition. The second phase of the 
investigation involved a theoretical analysis of the 
effect of several artificial stabilization devices on the 
lateral stability and controllability of a straight-wing 
airplane (hereafter referred to as airplane B) for a high- 
speed flight condition. This phase of the investigation 
included studies of the effects of wide variations of 
Can Cn Ce” 


np 

** In order to simplify the analysis, most of the calculations 
were made assuming idealized artificial-stabilization systems, but 
a few check calculations were made taking into account the small 
constant time lag of the stabilization device used in the tests. 
Extensive calculations were not made by this method, however, 
because of the extremely laborious process involved and because 
a systematic determination of the effect of time lag on stability 
throughout the variation of the four derivatives was considered 
beyond the scope of the present investigation. 
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Although the results do not apply directly to air- 
planes or flight conditions other than those investigated, 
it is believed that the trends presented might give a 
qualitative indication of the general effects of large 
variations of these stability derivatives. 


SYMBOLS AND COEFFICIENTS 


The forces and moments are referred to the stability axes, which 
are defined as an orthogonal system of axes intersecting at the 
airplane center of gravity in which the Z-axis is in the plane of 
symmetry and perpendicular to the relative wind, the Y-axis is 
in the plane of symmetry and perpendicular to the Z-axis, and 
the }-axis is perpendicular to the plane of symmetry. The 


symbols and coefficients aré defined as follows: 


foe = angle of roll, rad 

y = angle of yaw, rad 

3 = angle of sideslip, rad 

, = yawing angular velocity, rad. per sec. or deg. per sec. 
(dy /dt) 

p = rolling angular velocity, rad. per sec. or deg. per sec 
(do/dt) 

p = rolling angular acceleration, rad. per sec. per sec 
(d*¢/dt?) 

V = air speed, ft. per sec. 

p = mass density of air, slugs per cu.ft. 

q = dynamic pressure, lbs. per sq.ft. [(1/2)pV?} 

b = wingspan, ft 

S = wing area, sq.ft 

A = aspect ratio (b?/S) 

t = time, sec 

y = nondimensional time parameter based on span ( Vt/b) 

N = yawing moment, ft.lbs 

D = differential operator (d/ds,) 

1/T:/, = damping factor, reciprocal of the time to damp to 
half amplitude (positive values indicate stability 
and negative values indicate instability or time to 
double amplitude), 1/sec. 

W = weight of airplane, lbs 

W’/S = wing loading, lbs. per sq.ft 

m = mass of airplane, slugs (W/g) 

g = acceleration due to gravity, ft. per sec. per sec 

Mb = relative-density factor (m/pSb) 

n = inclination of principal longitudinal axis of airplane 
with respect to flight path, positive when principal 
axis is above flight path at the nose, deg 

y = angle of flight path to horizontal axis, positive in a 
climb, deg. 

ky, = radius of gyration in roll about principal longitudinal 
axis, ft. 

kz. = radius of gyration in yaw about principal vertical axis, 
ft. 

Kx, = nondimensional radius of gyration in roll about prin- 
cipal longitudinal axis (kx,/b) 

Kz, = nondimensional radius of gyration in yaw about prin- 
cipal vertical axis (kz,/b) 

Ky = nondimensional radius of gyration in roll about longi- 
tudinal stability axis (VV Ky,? cos*n + Kz,’ sin?) 

Kz = nondimensional radius of gyration in yaw about ver- 
tical stability axis (‘VV Kz,? cos? + Ay,? sin? 7) 

kyz = product-of-inertia factor [kz,? — kx,” sin 7 cos 7] 

Kyz = nondimensional product-of-inertia parameter 

Kz, — Kyx,? sin 7 cos 7] 

C, = trim lift coefficient (W cos y/q5S) 

C = rolling-moment coefficient (rolling moment /g.Sd) 

c. = yawing-moment coefficient (yawing moment/q.Sb) 


Cy = lateral-force coefficient (lateral force /q.S) 
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Cig = 0C;/08 

Gc 8 = 0C,/08 

Cy, = OCy/O0B 

Cc = 0C,/(rb/2V 

G, = OC,/0(pb/2V 

Ci = 0C,/d(pb/2V 

Ci = 0C,/O(rb/2V 

Cup = OC, /O( pb?/4V2 

Cy, = dCy/d(pb/2V 

Cy = OC; /d0(rb/2V 

Cup = OC,/06 

6 = rudder deflection, rad 

6,/r, = gearing ratio, rate of change of rudder deflection wit} 
yawing angular velocity 

6,/p = gearing ratio, rate of change of rudder deflection wit 
rolling angular velocity 

6 p = gearing ratio, rate of change of rudder deflection wit 


rolling angular acceleration 


SCOPE OF THE INVESTIGATION 


The basic aerodynamic and mass characteristics oj 
the two airplanes for the conditions chosen for this 


investigation are presented in Table 1, and, unless 
otherwise noted, all calculations were based on thes 
values. The dynamic lateral stability characteristics 
of each of these airplanes were calculated by the con- 
ventional methods of reference 1. This investigation 
was made to study the trends of the various rotan 
stability derivatives which may be useful in artificial 
stabilization systems, and, in order to simplify the 
study, idealized artificial stabilization devices were 
assumed. In practical application, of course, the 
effect of the actual characteristics of the artificial sta- 
bilization system should also be considered. 

Large variations of each of the rotary derivatives 
Crs Cipy Cup,» and C), were accomplished in the tests of 
the free-flight tunnel model of airplane A by having the 
response signal from a small rate gyro located in the 
operate the servomechanisms that moved the 
rudder or With the particular artificial 


stabilizing devices used in these tests, it was possible 


model 
ailerons. 


to obtain essentially pure derivative changes because 
of the relatively small phase lag compared with. the 
natural frequencies of the oscillations, because the gyr 
was mounted so that it would respond to either a pure 
yawing velocity or a pure rolling velocity about the 
stability axes, and because the rudder and ailerons wert 
interconnected to produce either pure yawing moments 
or pure rolling moments. * 

The geometric and mass characteristics of airplane 
A were adjusted to produce a neutrally stable oscilla 


tion for the basic condition. The flight tests of the 


* Deflection of the all-movable vertical tail to produce th 
rolling and yawing moment derivatives also produced changes ! 
the lateral force derivatives, Cy,, the lateral force due to rolling 
and Cy,, the lateral force due to yawing. In the calculations 
however, these changes in Cy, and Cy, were neglected because 
preliminary calculations indicated that even the largest changes 
in these derivatives did not appreciably affect the calculated 


results. 
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EFFECTS OF VARIATIONS 

scale model were made for exactly the same condi- 
tion for which calculations were made. In order to 
provide a direct comparison with the calculated os- 
cillatory stability, the period and damping values ob- 
tained from these experiments have been scaled up to 
represent those of the full-scale airplane. 

Since damping of the lateral oscillation is the primary 
function of any artificial stabilization system, only 
derivative variations in the directions that produce 
improvement in oscillatory stability will be discussed. 
More detailed studies of the effects of artificial sta- 
bilization on each of the airplanes herein discussed are 


to be found in references 2-4. 
RESULTS AND DISCUSSION 


Effect of Cn, 


Airplane A.—The effect of large variations of C,, on 
the dynamic lateral stability of airplane A is presented 
in Fig. 1. These data‘are in the form of period and 
damping of the lateral oscillation and damping of the 
aperiodic or nonoscillatory modes of motion. The 
damping of both the oscillatory motion and the aperi- 
odic motion is expressed in terms of the damping factor 
| T.,,, the reciprocal of the time to damp to half ampli- 
tude. Positive values of this damping factor indicate 
stability, and negative values indicate instability (time 
to double amplitude). Experimental values of the 
period and damping of the short-period lateral oscil- 
lation determined from flights of the model with various 
values of C,, are indicated by the symbols on Fig. 1. 
These values of period and damping were scaled up 
from values obtained from time histories of flights of 
the model such as the typical time histories presented 
in Fig. 2. 

A comparison of the calculated and experimental 
values of period and damping of the lateral oscillation 
indicates good agreement for the various values of C,, 
tested. In the basic condition with a C,, value of —0.2 
indicated by the solid symbol on Fig. 1) the airplane 
had neutral oscillatory stability. Moderate increases 
in the value of C,, in the negative direction caused a 
marked improvement in damping of the short-period 
lateral oscillation. Maximum damping of the oscilla- 
tion was obtained with a C,,. value of about —2.0. 
Further increases in the value of C,,, caused a reduction 
indamping of the lateral oscillation. In fact, oscillatory 
instability was obtained for values of C,,, above —6.0. 
ltis therefore apparent from both the flight-test results 
and the calculations that there exists some optimum 
value of C,,, for best damping of the lateral oscillation. 
For the higher values of C,,, (—3.0 to —7.0) the flight 
tests indicated that the lateral motion of the model 
progressively changed from the normal ‘‘dutch-roll’’ 
iscillation to a pendulum type of oscillation which was 
predominantly roll and sidewise displacement. The 


tme histories of Fig. 2 show that, at a C,, value of 
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Calculated effect of Cy, on stability and comparison 
with experimental data for airplane A 


Fic. 1 


—7.2, the ratio of yaw to roll is approximately one- 
This 
decrease in the ratio of yaw to roll is attributed to the 
fact that increasing the damping in yaw causes partial 
restraint of the yawing motions. 

In order to determine the significance of the increase 


half the value obtained in the basic condition. 


in the damping of the lateral oscillation obtainable by 
varying C,,, the results were compared with the Air 
Force and Navy requirements for damping of the lateral 
oscillation. This comparison shows that airplane A 


would meet the requirements with values of C,, be- 


tween —0.7 and —4.6. The flight behavior in this 
range of C,, was very good, and there was no indica- 
tion of the decreased maneuverability or ‘‘stiffness’’ 


in making turns, which has been experienced with some 
airplanes equipped with yaw dampers. Steady turn- 
ing maneuvers to determine the effect of C,, on the 
maneuverability could not be made in the free-flight 
tunnel because of the restrictions imposed by the size 
of the test section. 

The calculated stability in the negative C,, range 
indicates that the aperiodic modes of motion merge and 
form a second oscillation for C,,, values between —2.0 
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Io increased negatively despite the lightly damped or yp. 
ee * PERIOD, SEC stable oscillations at the higher values of C,,,. 
ti. x Airplane B.—Calculations of the effect of large var; 
aa tae oer ations of C,, on the dynamic lateral stability for gir. 
ba ssenonc} i plane B are presented in Fig. 3. Increasing C,, from 
4 4 the basic airplane value of —0.4 caused definite im- 
yy, . DAMPING FACTOR, op! provement in the damping of the lateral oscillation up 
“A / to a value of C,, of approximately —5.4. Further j; 
ff Oe eMODE / je creases in C,, above this value caused a reduction : 
Pa S tie a damping of the lateral oscillation. 
SS — nen a” Above a C,, value of —6.7, the oscillatory mode 
10 -9 -8 -7 -6 -5 -4 -3 -2 -I 0 


YAWING MOMENT DUE TO YAWING, Cnr 


Fic. 3. Calculated effect of C,, on stability for airplane B. 
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Fic. 4. Calculated effect of Ci, on stability and comparison with 


experimental data for airplane A. 
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Fic. 5. Calculated effect of C,, on stability and comparison with 
experimental data for airplane A. 
and —5.2. This oscillation was so heavily damped 


Al- 
though it was not possible to measure the damping of 


that it was never observed in the model flights. 


the aperiodic modes of motion from the flight tests, the 
pilot was aware of the increasingly good spiral stability 
of the model as the value of C,, was increased nega 
tively, since the model had an increasingly strong tend 
ency to fly uncontrolled along the centerline of the tun- 
nel test section despite the natural gustiness of the air- 
flow. This increase in spiral stability is shown by the 
flight records of Fig. 2 which indicate that it was possible 
to obtain longer uncontrolled flight records as C,,, was 





broke down to form two more aperiodic modes. One 
of these aperiodic modes merged with the mode origin- 
ally designated as the spiral mode at a C,, value of 
—9.0 to form a second oscillation. This second oscil 
lation became increasingly less stable as C,, was in 
creased beyond — 9.0. 

that B met the 


damping requirements from a C,, value of —0.7 up to 


It may be noted airplane 
the highest values investigated. 
Generally speaking, the trends produced by varying 
C,, with airplane B were comparable to those for air- 
plane A. That is, as C,, was varied negatively for 
both airplanes, the damping of the oscillatory modes 
increased up to some optimum value and then de- 


creased, while the aperiodic modes showed increases | 


in damping. 


Effect of C,, Variation 


The effect of large variations of the damping-in- 
roll derivative C;,, on the dynamic lateral stability 
characteristics of airplane A are presented in Fig. 4 
The experimental period and damping values indi- 
cated by the symbols on the figure were obtained from 
flight records similar to those previously presented for 
ae 
mental data with the calculated oscillatory stability for 
lp The 


damping of the oscillation with high negative values oi 


A comparison of the limited amount of expert- 


small values of C,, shows fairly good agreement. 
C,, could not be determined experimentally because ol 
the abnormal flight behavior of the model in that range 
The rolling mode of motion was so heavily damped 
that the model was essentially restrained from rolling 
This tendency toward restraint imposed by high nega 
tive C,, values, although undesirable from the stand- 
point of maneuverability, resulted in flights requiring 
little corrective control when the model was flown ina 
The restraint in roll exper! 
enced in the flight tests is indicated in the calculated 


steady wing-level attitude. 


results that show that one of the aperiodic modes 0! 
motion (the rolling mode) became increasingly stable 
as C;, was increased negatively. 


As C, 
the basic value of —0.3 (indicated by the solid symbol 


was increased in the negative direction from 


on Fig. 4), the damping of the short-period oscillation 


rapidly improved for values of C,, up to about —0.! 


Further negative increases of C,, (to —7.3) produce 
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ZEPrPPECTS OF VARIATIONS tN 
no further damping of the lateral oscillation. The air- 
would just barely meet the damping requirements 


plane ; 
It is appar- 


jor values of C;, between —0.7 and —7.3. 
ent from these results that almost all of the damping 
added to the system by increasing the damping-in-roll 
derivative C;, is absorbed by the aperiodic rolling mode 
of motion, while only a small part of the added damp- 
ing goes into the damping of the oscillation. 

Effect of Cn, 

Airplane A.—The dynamic lateral stability charac- 
teristics of airplane A as affected by large variations of 
the yawing-moment-due-to-rolling derivative C,, are 
presented in Fig. 5. Model flights were possible only 
for a limited range of C,,, because, as will be discussed 
later, excessive amounts of aperiodic instability were 
encountered. The experimentally determined values 
of the period and damping of the oscillation are in fairly 
good agreement with the calculated characteristics. 
~ These results indicate that, as C,,, was varied in the 
positive direction from the basic condition, there was 
a rapid improvement in the damping of the lateral os- 
cillation and that the damping requirements were satis- 
fied with all C,, values greater than 0.1. This im- 
provement in damping, however, was not obtained 
without a loss in stability elsewhere in the system. The 
calculations for values of C,,, between 0 and 0.4 reveal 
that this improvement in oscillatory stability was ob- 
tained at the expense of the normally well-damped 
rolling mode. Even though this mode of motion re- 
mained stable up to the point of its merger with the 
other aperiodic mode at a Ow value of about 0.4, it was 
considerably less damped than in any other flight con- 
dition experienced in these tests. 

The model of airplane A became touchy to fly as C,, 
was increased up to 0.4 and became extremely difficult 
to control for values of C,,, greater than 0.4. Sustained 
flight was impossible for C,,, values greater than 1.0. 
As the value of C,,, was varied through this range (—0.1 
to 1.0), the pilot complained of an increasingly strong 
tendency for the model to go into a tight turn in re- 
sponse to normal aileron control. To the pilot, this 
tendency appeared to be a severe case of spiral insta- 
bility. The results of calculations, however, show that 
the stability of the spiral mode of motion remained 
unchanged up to a C,, value of 0.6, while there was a 
rapid decrease in stability of the rolling mode of motion. 
It appears, therefore, that a decrease in stability of the 
rolling mode of motion might sometimes be mistaken 
for spiral instability. 

The two aperiodic modes merged to form a long- 
period oscillation that became unstable at a C,,, value 
of 0.6. This oscillation was not observed in the flight 
tests because of its extremely long period of over 40 
sec. Immediately after it became unstable, the long- 
period oscillation broke up to form two new aperiodic 
modes of motion, one of which became increasingly 


unstable as C,,. was increased further. 


np 
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Calculated effect of C,, on stability for airplane B with 


two values of Cip. 
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The results of these tests and calculations indicate 
that the derivative C,,, might possibly be useful for re- 
distributing the natural damping of an airplane so as 
to obtain a well-damped lateral oscillation and, at the 
the same time, obtain more maneuverability in roll. It 
is obvious, however, that the use of Cy», is limited to 
values less than those that would cause the undesir- 
able aperiodic motions experienced in these tests. 

Airplane B.—The calculated dynamic lateral sta- 
bility characteristics of airplane B, as affected by large 
variations of C,,, for two values of C;, are presented in 
Fig. 6. The results of the calculations with the basic 
value of C;, (—0.40) are shown in Fig. 6a, and the re- 
sults with C,, increased to —0.80 are shown in Fig 
6b. 

As in the case of airplane A, the results of the calcu- 
lations with the basic value of C,, show that, as C,, is 
increased positively, the damping of the lateral oscilla- 
tion improves rapidly and the rolling mode becomes 
less stable. The damping requirements are satisfied 
for all values of C,,, greater than 0.1. 

It should be pointed out that for this airplane an in- 
crement of C,,, of 0.62 resulted in the same damping of 
the oscillation as did an increment to C,, of —1.6. 
Physically, this indicates that the gearing ratio 6r/r 
required to obtain this damping with a C,, device is 
approximately 2.5 times as great as the gearing ratio 
ér/p required with a C,, device. In many 
however, the rudder motion required by the C,, 
will be larger than that required by the C,, device be 
2.5 times 


cases, 


device 


cause the rate of roll will often be more than 
as great as the rate of yaw. 

The results previously presented for airplane A indi- 
cated that varying the stability derivative C;, had little 
effect on the damping of the lateral oscillation but 
caused large increases in the damping of the rolling 
mode. From these results it appeared probable that, 
if C,, were increased simultaneously with C,,, greater 
damping of the short-period oscillation might be ob 
tained before the rolling mode had become too lightly 
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Fic. 7. Calculated effect of Ci, on stability and comparison with 
experimental data for airplane A. 
damped. Therefore, a case in which both C,, and 


C,, were increased was investigated. The effect of 
C,, on the stability of the lateral modes of motion when 
C,, was increased to —0.80 is shown in Fig. 6b. The 
results of these calculations show that, for a C,,, value 
of 1.1, the damping of the short-period oscillation was 
appreciably increased when the damping of the rolling 
mode was the same as for the basic condition. Forma- 
tion of the long-period oscillation was delayed to a higher 
value of C,,. These results indicate that, by increasing 
C,, aud C,, simultaneously, it might be possible to in- 
crease the damping of the lateral oscillation without 
adversely affecting the lateral motions in other re- 


spects. 


Effect of Ci, 


Airplane A.—The effect of large variations of C;, on 
the dynamic lateral stability characteristics of air- 
plane A is shown in Fig. 7. The calculations indicate 
that increasing the value of C;, in a positive direction 
improved the damping of the lateral oscillation. The 
requirements for damping of the lateral oscillation were 
satisfied with values of C;, greater than 0.8, but, for 
these cases, the spiral mode of motion was so unstable 
that the: overall flight behavior was unsatisfactory. 
Attempts to measure the damping of the oscillation in 
the positive C;, region were not successful, because the 
spiral instability, which resulted from the increase in 
C,,, required that the pilot apply corrective control 
immediately after a disturbance to prevent the model 
from crashing. 

Qualitatively, the theory is in agreement with ex- 
perimental results in that it correctly predicts the 
existence of stability or instability, but quantitatively 
the agreement is poor. The reason for this poor quan- 
titative agreement has not been established. * 


* The calculated values of the damping in the negative C;, 
range are not in agreement with the experimental results. Cal- 
culations were made including a constant time lag of 0.05 sec. 
in an effort to explain this discrepancy. These calculations were 
in much better agreement with the experimental results, indicat- 
ing that the discrepancy between the measured and calculated 
values of damping may therefore be attributed, at least partly, 
to the effect of time lag in the stabilization device 


The flight tests showed that there was a slight jm. 
provement in the damping of the lateral oscillation 
when C;, was increased to 0.3, but the model became 
more difficult to fly despite this increase in damping, 
Both the calculated results and the flight-test results 
indicate that spiral instability occurred at a C,, value of 
0.3 and became more severe as C;, was further increased. 
The aperiodic instability observed in the flight tests 
with positive C;, appeared to the pilot to be similar to 
the difficulty experienced with small positive values of 
Crp (up to 0.6). 
calculations, however, that the spiral mode became 


It is interesting to note from the 
unstable when C), was increased, whereas it was the 
rolling mode that became less stable when C,, was 
increased. ) 

It appears from these results that variation of the 
derivative C;, offers little hope for improving the flight 
characteristics of an airplane. This derivative, how- 
ever, may be used to redistribute the damping between 
the oscillatory and spiral modes if surplus damping of 
the spiral mode is initially present. Preliminary cal 
culations for this airplane have indicated that the 
damping obtained with C,, alone could be improved ap- 
preciably by utilizing C;, to redistribute part of the 
excess damping of the spiral mode to the oscillatory 


mode of motion. 


Effect of C.., 


The possibility of increasing the damping of the 
lateral oscillation of airplane B through use of an arti- 
ficial stabilization system sensitive to rolling acceler 
ation and geared to the rudder was considered. This 
system introduces into the yawing moment equation 
the term 


—(O.N 06,) [06,/0(d?o/dts) | (d*o/dt?) 


which may be considered as an increment to the term 
mkxz (d*o/dt*) in that equation. This system was 
investigated as a consequence of the results presented 
in references 5 and 6, which indicated that the product- 
of-inertia parameter kyz has a stabilizing effect on 
damping of the lateral oscillation if the principal longi- 
tudinal axis is above the flight path at the nose of the 


airplane. The results of reference 5 also indicated that 
the term involving kyz in the yawing equation was 


primarily responsible for the stabilizing effect. There 
fore, it was believed that, if the rudder were deflected 
in proportion to the rolling acceleration, an improve- 
ment in damping of the lateral oscillation would result, 
since the aerodynamic forces produced by the movemeut 
of the rudder would have the same effect as increasing 
the value of kyz in the yawing-moment equation. 
The parameter added to yawing-moment equation may 
be expressed in nondimensional stability derivative 
form as —(1/4)C,,,D,*@. 

The effect on the stability of the lateral modes of 
is increased from 0) to 

The damping of the 


motion as the parameter C,,, 
—0.08 is presented in Fig. 8. 


lateral oscillation increases rapidly as C,,, is increased, 
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EFFECTS OF VARIATIONS 


The 


spiral mode is essentially insensitive to changes in 


put the period becomes increasingly shorter. 


(,,, whereas the rolling mode becomes considerably 
less damped as this parameter is increased. 


of —0.009, the lateral oscillation has ap- 


For a 
value of Cy 
proximately the same damping as for the cases con- 
red previously for airplane B with a C,, value of 
A point of interest, how- 


side 
~2,0 or a C,, value of 0.6. 
ever, is that the gearing ratio 6,/p necessary to ob- 
tain this value of C,,, is of the order of 1 20 the gearing 
ratio required for the C,, and C,,, cases. 

The use of C,,, produced good damping of the lateral 
oscillation, but, as will be shown later, this derivative, 
like some of the other derivatives, also produced unde 


sirable effects on controllability. 


Effect of C.,, Cy, and C, , on Controllability 


Each of the artificial stabilization devices discussed 
in conjunction with airplane B appears capable of im- 
proving considerably the damping of the lateral oscilla- 
tion. The acceptability of each of the devices investi 
gated, however, would depend to a large extent on the 
The 


lateral motions subsequent to application of a rolling 


response of the artificially stabilized airplanes. 
moment coefficient of 0.01 were therefore calculated 
for the case of the basic airplane and for the cases of the 
airplane equipped with three different artificial sta 
bilization devices. Each of these stabilizing devices 
increased the value of one of the following stability 
1.20, Cr, = 0.60, 


The motions were obtained by the con 


derivatives to the values: C,, 
Ow 0.005. 
ventional methods described in reference 1 and are 
presented in the form of time histories of the rolling 
and yawing velocities in Fig. 9. 

A comparison of the motions obtained with the basic 
airplane and the motions with the C,, device indicates 
that there is practically no change in the response 
characteristics. These results may appear to conflict 
with the previously mentioned results obtained with 
some full-scale airplanes equipped with yaw dampers 
in that no appreciable decrease in yawing velocity 
with increased C,, is predicted by the calculations. 
The results shown, however, are for an aileron rolling 
maneuver, and, consequently, the stiffness in yaw should 
not be expected to be as evident as in a turning ma 
netver. 

With the C, device, there is a large increase in the 
rolling velocity which is attributed partly to the lower 
damping of the rolling mode and partly to the large 
favorable yawing motion produced by this derivative. 

With the C,,, device, there is a large adverse yawing 
motion that greatly reduces the rolling velocity during 
the initial portion of the maneuver and which appar- 
ently causes the final rolling velocity to be about the 
same as that for the basic condition, even though the 
damping of the rolling mode is less than half of that for 


the basic condition. It is interesting to note that C,, 
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Fic. 9. Calculated effect of Cn,-, Crz, and Cy, on the response of 
airplane B to the applied rolling moment 
and C,,, the yawing moments produced by rolling 


velocity and acceleration, produce opposite effects on 
controllability, even though both increase the damping 


of the lateral oscillation. 


CONCLUSIONS 


Although the following results do not apply di 
rectly to airplanes or flight conditions other than those 
investigated, it is believed that they provide a quali 
tative indication of the general effects of large vari 
ations of the stability derivatives: 

(1) 
large increase in damping of the lateral oscillation with 


The only single derivative that provided a 


out adversely affecting other flight characteristics was 
C,,. Flight tests of some full-scale airplanes, how 
ever, have indicated that even this derivative might 
result in unsatisfactory stiffness in turning maneu- 
vers. 

(2) 
tendencies before 
damping to the lateral oscillation. 


or C,, caused undesirable spiral 
of 


Increases in —C, 
> 


Increases in C,, 


adding a substantial amount 


caused objectionable stiffness in roll, and increases in 
—C,,; caused objectionable adverse yawing with a con- 
sequent reduction in rolling velocity in aileron rolls. 
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Investigation into the 


Operating Cycle of a 


Two-Dimensional Supersonic Wind Tunnel’ 


G. V. BULLt 
Institute of Aerophysics, Unwersity of Toronto 


SUMMARY 


The operating cycle of an intermittent supersonic wind tunnel, 
while probably recognized, has never been defined rigorously. 
In this paper, the wind-tunnel operating cycle is divided into five 
distinct phases based on the character of the flow along the wind 
tunnel axis. These phases are, in cyclic order, the nonstationary 
starting phase, the quasi-stationary starting phase, the stationary 
phase, the quasi-st itionary closing phase, and the nonstationary 
closing phase. The first of the starting phases is treated both 
theoretically and experimentally for the case of the opening valve 
located downstream from the working section and experimentally 
for the case of the valve located upstream from the nozzle throat 
These theoretical studies are based on an assumed wave pattern 
(subsequent to the breaking of the diaphragm) determined from 
results of shock-tube investigations and by application of an ap 
proximate theory for the solution of two-dimensional nonstation 
ary flow problems due to Guderley. The quasi-stationary start 
ing phase is studied theoretically by extending the stationary 
supersonic-subsonic jet theory in which the jet exit static pres 
sure is less than the pressure in the region into which the jet is 
issuing. The theoretically postulated state of affairs during these 
phases is verified by obtaining experimentally a measured series 
of shadowgraphs and, thereby, determining the time history of 
the wind-tunnel starting process. The effect of condensation of 
water vapor on these phases is studied experimentally. 


INTRODUCTION 


— STANDARD METHOD OF TREATING the theory 
of intermittent supersonic wind tunnels of the 
type shown in Fig. la is based on application of isen- 
tropic steady-state equations to channel flow. In this 
manner the state conditions at the end of the Laval 
nozzle (that is, in the working section) can be related 
to the inlet stagnation state conditions, and the Mach 
Number can be determined as a function of the throat- 
to-working-section area ratio. For any specified Mach 
Number, the ratio of the inlet stagnation pressure to 
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reservoir stagnation pressure determines the horse- 
power requirements of the tunnel and/or the duration 
of supersonic flow in the working section. With no 
recompression downstream of the working section, the 
pressure ratios inlet stagnation pressure to working- 
section static pressure and inlet stagnation pressure to 
reservoir stagnation pressure would be equal. How- 
ever, in practice, pressure recovery between the work- 
ing section and vacuum reservoir is accomplished by 
means of shock phenomena, and the critical reservoir 
stagnation pressure (that is, the reservoir stagnation 
pressure that just maintains supersonic flow in the 
working section) is less than the inlet stagnation pres- 
sure by the amount of pressure loss (or entropy increase) 
through the shock system. By assuming recompres- 
sion through a normal shock wave (i.e., assuming no 
supersonic diffuser) the required pressure ratio for any 
working-section Mach Number can readily be found.! 


In this paper it is desired to consider the complete 
operating cycle of this type of wind tunnel and to con- 


sider in detail the starting processes. 


WIND-TUNNEL OPERATING CYCLE 


It is customary to divide the wind-tunnel operating 
cycle into three parts —-viz., the starting processes, the 
steady-state processes, and the closing processes. For 
purposes of theoretical analysis it is desirable to sub- 
divide the wind-tunnel operating cycle into phases 
based on the character of the flow along the tunnel axis 
With this division the operating cycle 
(1) the nonstationary starting 


(see Fig. Ic). 
consists of five phases: 
phase; (2) the quasi-stationary starting phase; (3) 
the stationary phase; (4) the quasi-stationary closing 
phase; and (5) the nonstationary closing phase. The 
nonstationary starting phase (during which flow condi- 
tions are varying with time along the whole tunnel axis) 
commences with the opening of the valve and termi- 
nates with choking of the Laval nozzle throat. The 
quasi-stationary starting phase (during which flow 
conditions are constant over part of the tunnel exis) 
commences with choking of the nozzle throat and is 
considered to terminate when flow conditions have 
become steady along the tunnel axis up to the entry 
throat to the vacuum With the termi 
nation of the quasi-stationary starting phase, the tunnel 
operation enters the stationary phase, which is usually 
many orders of magnitude greater in duration than 


reservoir. 
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any of the previous phases. This stationary phase is 
terminated when the stagnation pressure in the vacuum 
reservoir has risen to such a value that the flow is sonic 
at the entry to, and subsonic inside, the reservoir. In 
this case a normal shock wave appears at the entry 
throat, and with further rises in reservoir stagnation 
pressure the shock moves upstream into the working 
section and is finally dissipated in the nozzle throat. 
The quasi-stationary closing phase is considered to 
commence with the termination of the stationary phase 
and to terminate with the dissipation of the shock wave 
at the nozzle throat. The nonstationary closing phase 
is the final phase of tunnel operation during which the 
flow all along the tunnel axis is purely subsonic and not 
steady with time. This phase terminates when the 
stagnation conditions of the reservoir and inlet become 
The stationary and closing phases have been 
The remainder of this 


equalized. 
dealt with in the literature. 
paper will consider only the nonstationary starting 
phase and the quasi-stationary starting phase. 


THE NONSTATIONARY STARTING PHASE 


The nature of the phenomena associated with this 
phase depends on the location of the opening valve 
with respect to the throat. Hence, we must consider 


two possible cases. 
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(1) Valve Located Downstream from the Throat 


If the valve is located as shown in Fig. la and can }y 
assumed to open instantaneously (in practice achieved 
by using a cellophane diaphragm), then it is possihk 
to treat this phase by using the shock-tube wave mode] 
shown in Fig. 1b. This wave model postulates that. 
subsequent to rupturing of the diaphragm, a centered 
rarefaction wave proceeds upstream through the work 
ing section and Laval nozzle and a shock wave and con 
tact surface proceed downstream into the vacuum res 
ervoir. It is evident that these two phenomena wil] 
be dissipated in the vacuum reservoir and require no 
further consideration. Consequently, the nature of 
the axial distribution of flow parameters can be deter 
mined by constructing the nonstationary characterist 
diagram. This construction (Fig. Ib) can be carried 
out using a method due to Guderley,* provided the flow 
can be assumed isentropic and that across any cross 
section the flow parameters are constant. Under thes 
assumptions the equations of the characteristics may 


be writtenas. Ct, 


| dp du din F dy 
} au / along uta 


pa dt dt ax di 
and C 
ldp du din F dx 
au( ) along u a 
pa dt dt dx al 


With the tunnel geometry shown in Fig. la, the char 


acteristic field bounded by the C” characteristic OA 
and the C 
one-dimensional so that within this region Eqs. (1) and 


Mu 
| | 
a 


characteristic through A may be treated as 


(2) become: C 


CLA 
const. along — u+a (6 
; 
( 
and C 
2 a i 
im l 
Y 1 \a ad 
dx 
const. along — u a t 
at 


Along the C~ characteristic OA, 


Along the C+ characteristic Az, 


Z a u 
—1})+— =90 
y — l\do ao 
Thus, for any point 7 (on OA) using the relationship for 


the C~ characteristic, 


(dx/dt =U—a 
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OPERATING 


we can write 


2 ¥ l/dx 
a a ) 
y+ 1 y+ 1\dt 
(dx/dt) + a 
Thus, by choosing any arbitrary slope (dx/dt)—; for 


the C~ characteristic, we can find the values of the flow 
parameters at 71, the point of intersection of the C 
characteristic through A, and the selected C~ charac 
teristic. By choosing the appropriate C 
can evaluate any point on the C+ characteristic through 


A. Eqs 


slope, we 


|) and (2) may be written in the form 


dr/dt (au/aoF) (dF /dx) 
du/dt (au/aoF) (dF /dx) 


ind then can be put in difference form 


au dF re 

Ne Na 2 (le ts \) 
aok dx 
au dF 

bh Ma . (te tp) (6 
' aoF dx 


Aisa known point on a C* characteristic, B is a known 
point on a C~ characteristic, C is the point of intersec 
tion of the C* characteristic through A with the C 
characteristic through B.) 

Using the initial conditions and applying Eqs. (5 
and (6), the remainder of the characteristic diagram 
can be constructed by extrapolation along the char 
acteristics; the characteristic diagram for the case in 
vestigated is presented in Fig. lb. ,. From this charac 
teristic diagram it is seen that the effect of converging 
walls is to decrease the absolute value of the C char 
acteristic slope (dx/dt) — as the characteristic is extended 
into the convergent duct; the slope of any such charac 
teristic is a minimum at the throat and decreases in 
the divergent section. From the characteristic rela 
tionships of Eqs. (2), it is seen that u (the induced par 
ticle velocity) must increase as the characteristic pro 
ceeds into the convergent section. From energy con 
siderations if « increases, the velocity of sound de 
creases and, since the flow has been assumed isentropic, 

bp = po(a/a)”'7~” 
and therefore the pressure must decrease. We can 
state then that each characteristic traveling upstream 
causes a negative pressure gradient upstream and post- 
tive pressure gradient downstream of the throat, in 
both cases the gradient being considered in the direc 
tion of particle travel. Consequently, separation will 
occur at the throat, and the induced mass flow will issue 
irom the nozzle throat as a jet bounded by free stream 
lines. The gas flow contained between the tunnel 
walls and the free stream lines will be characterized by 
The induced mass flow is 
This unsteady phe- 


rotationality or vorticity. 
continually increasing with time. 
nomenon is termed the ‘‘nonstationary jet effect’’ and in 
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Fig. 2a a shadowgraph taken about halfway through 
the nonstationary phase shows this effect in the throat 
of a Laval nozzle designed to produce Mach Number 2.9 
flow in a 5 These effects 


are cumulative; after a small portion of the rarefaction 


by 5-cm. working section 


wave has passed the velocity of sound is reached in the 


throat and choking occurs. Once the nozzle has 
choked, the pressure and particle velocity remain con 
stant in and upstream from the throat—that is, the 
steady state has been reached along part of the tunnel 
axis so that the overall axial character of the flow may 


be considered quasi-stationary. 


2) Valve Located Upstream from the Throat 


When the opening valve (or diaphragm) is located 
upstream from the nozzle throat, the phenomenon 
associated with the nonstationary phase changes from 


an isentropic rarefaction to a discontinuous shock phe 
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the nonstationary and quasi-stationary starting phase with the 
diaphragm located upstream from the throat. 


nomenon. If the valve is located in a portion of the 
duct with parallel walls and opens instantaneously 
(again achieved by the use of a cellophane diaphragm), 
then subsequent to opening a normal shock wave is 
formed‘ which travels downstream at a velocity in ex- 
cess of the velocity of sound and is followed by a con- 
tact surface traveling downstream at the induced par- 
ticle velocity. 

When the normal shock enters the converging sub- 
sonic portion of the Laval nozzle, it will suffer 
of regular or Mach reflections depending on the initial 
If the conditions are such that Mach 


a series 


pressure ratio. 
reflection occurs, then the following qualitative dis- 
cussion (determined experimentally) applies. The 
initial plane shock on entering the convergent channel 
suffers Mach reflection from each wall giving rise to a 
reflected shock and contact surface as depicted in Fig. 
fa. The induced particle velocity direction is indicated 
neglecting any curvature of the legs of the Mach shock. 
It must be remembered that the contact surface repre- 
sents a surface across which there is a discontinugty in 
the tangential velocity component, and, henee, the’ par- 
ticle velocity vector is not necessarily parallel to this 
fa the region interior to ABC is as- 
to the 
the region exterior to the shock ACD is also 
but in this 


surface. In Fig. 
sumed constant with a flow velocity parallel 
boundary; 
assumed to be a region of constant flow, 
case the flow velocity is parallel to the axis; the region 
interior to AECD is one in which the velocity vector 
lies somewhere between being parallel to the axis and 
As the initial shock travels 
point of reflection 


being parallel to the wall. 
through the converging duct, the 
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(the so-called triple point) moves away f.om the te. 
flecting wall and the contact surface gets larger—ij¢. 
the regions ABCE and AECD get larger (Fig. 4b) 
The triple points thus approach each other and even. 
tually coincide (Fig. 4c). During coincidence, since 
perfect symmetry is assumed, a contact surface of zero 
strength will be generated. This condition of sym. 
metry (coincidence of triple points) will persist for a 
short period of time during which the shock front wil] 
have traveled downstream. Clearly, then, a contact 
surface of zero strength will be generated parallel to the 
axis, the length of which will be determined by the 
magnitude of the excess forward velocity of the wave 
front over the relative cross velocity of the triple points, 
(If the wave-front velocity was equal in magnitude. to 
the relative cross velocity of the triple points, then this 
zero strength contact surface would be a point.) These 
conditions of symmetry will be destroyed when the 
triple points pass through each other and continue to 
travel in their original direction—that is, toward the 
wall opposite to the one from which they originally 
started. It should 
be noted in this case that, once the reflected shock ADC 


This condition is shown in Fig. 4d. 
crosses over the axis of symmetry, it cancels out the 
induced particle velocity due to the wave originating 
at the other wall so that the flow velocity vector be- 
When the 
fe), inverted 


hind it is again parallel to the tunnel axis. 
triple points collide with the wall (Fig. 
Mach reflection occurs and the triple points again move 
toward the axis of symmetry (Fig. 4f) and the process 
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OPERATING CYCLE OF 
is repeated (Fig. 4g), the number of reflections and 
inversions being determined by the strength of the 
initial shock and the length and slope of the con- 
verging walls. After the shock front passes the throat, 
reflection no longer occurs, but rather (as recently 
shown by shock tube interferograms taken at Prince- 
ton) an expansion wave causes the shock front to be 
curved. The reflected shock and contact phenomena 
break down because of dissipation beyond the throat 
(Fig. 4h). After the passage downstream of the initial 
wave front, the mass flow continues to increase (due 
to the rarefaction wave upstream of the throat) until 
choking occurs. The quasi-stationary phase is then 
entered. 

The case of regular reflection occurring was not stud- 
ied in detail; however, the results achieved indicated 
that this phenomenon was similar to the Mach reflec- 
tion case. Some shadowgraphs for this case are pre- 
sented in Fig. 5. 


THE QUASI-STATIONARY PHASE 


1) Valve Downstream from the Working Section 


During this phase, the nonstationary flow down- 
stream from the throat will develop into stationary flow. 
When the induced particle velocity in the throat has 
attained the sonic value, stationary flow will be estab- 
lished here and just downstream a slight expansion will 
make the flow supersonic. Since there is no stable 
shock-free transonic deceleration, recompression must 
be accomplished by some shock process. From the 
standard theory of supersonic-subsonic jets, in which 
the exit static pressure of the jet is less than the region 
into which the jet is impinging,’ it is seen that this shock 
process can consist of either a single normal shock wave 
or a pair of oblique shocks reflecting from a centerpoint 
of symmetry (cross shocks). However, reflection from 
the centerpoint of symmetry can occur only if the shock 
waves, recompressing from the exit jet pressure to the 
impinged region pressure and causing a flow deflection 
toward the axis, produce behind them a region of gas 
flow ata Mach Number sufficiently high to permit the 
existence of another shock capable of deflecting the flow 
back parallel to the axis. Since the Mach Number be- 
hind a shock (.V/.) decreases as the shock strength in- 
creases and is always less than the Mach Number 
ahead of the shock (./;), then, as J; approaches unity, 
the strength of the shock must approach unity if 1/2 is 
to be greater than i. But during the nonstationary 
phase already discussed, a finite pressure difference 
has already been established prior to choking. Con- 
sequently, immediately subsequent to choking, /, 
will differ only infinitesimally from unity (because of 
expansion adjacent to and on the downstream side 
of the sonic line), while a finite pressure difference will 
exist. Thus, reflection cannot occur from a central 
point of symmetry, and some form of curved normal 
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shock (assuming the sonic line to be curved) must be 
formed slightly downstream from the throat just subse- 
quent to choking. That this state of affairs exists 
experimentally is shown in Fig. 2b, which is a shadow- 
graph taken about 700 microsec. after Fig. 2a. It is 
evident that, after choking, the pressure downstream 
from the nozzle throat will start to fall because of the 
inequality of the fixed mass flow through the nozzle 
throat and the increasing mass flow into the vacuum 
(The mass flow into the reservoir increases 
throat into the vacuum 


reservoir. 
until the entry reservoir 
chokes. ) 

Thus, as the back pressure falls and approaches 
the pressure ahead of the shock, reflection from 
the central point of symmetry will occur—that is, 
the initial single normal shock will develop into a pair 
of crossed oblique shocks. The developed cross-shock 
pattern is shown experimentally in Fig. 2c (taken 1,100 
microsec. after Fig. 2b) and Fig. 2d (taken 200 microsec. 
after Fig. 2c). 

In Fig. 2d, the oblique shocks have rotated toward 
the Mach angle from the position they occupied at 
Fig. 2c. This is, of course, what is to be expected as 
the back pressure continues to fall; when the pressures 
on either side of the shock wave become equal, the 
shock becomes a Mach wave across which expansion 
will occur because of the divergence of the nozzle walls. 
Since this expansion can be considered to occur simul- 
taneously with the equalization of pressures on either 
side of the shocks, the pressure on the downstream side 
of the expansion wave will drop below the back pres 
sure, and, hence, at this point equilibrium is destroyed 
and a pair of shocks must be formed—i.e., the shock 
system can be considered to have moved down- 
stream. 

It is evident that this system, once established, will 
never depart far from equilibrium— that is, the oblique 
shock waves will always lie near the Mach wave. This 
is shown in Figs. 2d, 2e, and 2f, in which the develop- 
ment of the cross-shock pattern downstream can be 
traced. Another point, not shown experimentally but 
which emerges directly from this analysis, is that the 
rate at which the cross-shock system moves down- 
stream depends not only on the rate of fall of back pres 
sure but also on the rate of fall of the pressure ahead 
of the system due to the expanding walls. Only when 
the walls are parallel will the rate of development down- 
stream depend solely on the rate of fall of back pressure. 
Figs. 2g and 2h show the development of the jet through 
the working section and the establishment of steady 
flow in this region. 


2) Valve Located Upstream from the Throat 


No attempt has been made to treat this case in a 
theoretical fashion; however, experimentally it was 
observed (Fig. 4j) that, subsequent to choking, a cross- 
shock jet existed and developed downstream in an iden- 
tical manner to that previously discussed. 
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THE EFFECT OF CONDENSATION PHENOMENA 


Condensation of water vapor due to the falling pres- 


sure and temperature does not alter the nonstationary 


phase discussion because the critical conditions are not 
reached until after choking occurs. The quasi-station 
ary phase is altered insomuch as, in addition to the jet 
shocks, there exist condensation shock waves through 
which some recompression is accomplished. The con 
densation shock phenomena occur almost simultane 
ously with choking, as shown in Figs. 3a and 3b, and the 
resultant standard X-shock pattern moves downstream 
into the expansion part of the nozzle as the jet develops 


downstream (Fig. 3c). Hermann® has shown the 
condensation phenomena in a Laval nozzie for the case 


of stationary flow. 


CONCLUDING DISCUSSION 


The work presented here provides an insight into 
the unsteady part of intermittent supersonic wind 
tunnel operation. The overall picture has been estab 
lished and the starting phases investigated. The divi- 
sion of the operating cycle into five distinct phases 
permits systematic treatment of the problem. It is 
interesting to observe that it is possible to calculate 
the required pressure ratio across the diaphragm prior 
to rupture in order to permit operation in any one of 
these phases. The most interesting cases are the 
initial diaphragm pressure ratio required to choke the 
nozzle throat and the initial pressure ratio required to 
establish supersonic flow in the working section. The 
latter has been dealt with previously in the literature. 





The initial diaphragm pressure ratio, in order to choke. 
can be calculated by measuring the initial slope of the 
choking characteristic from Fig. lb and calculating 
(using this measured value) the velocity of sound and 
induced particle velocity on this characteristic (by the 
method outlined previously). Knowing these values. 
the shock-tube equation’ can be solved to find the re 
quired initial diaphragm pressure ratio. For this par 
ticular case the required initial pressure ratio was cal. 
culated® to be 1.06. 
value was not achieved because of the impossibility of 


Experimental verification of this 


obtaining a diaphragm rupture at such a low pressure 
difference across it. (Choking was experimentally 
observed with an initial pressure ratio as low as 


Looe) 
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A Review of Certain Analysis Methods 
for Swept-Wing Structures 


— 


M. L. WILLIAMS* 
California Institute of Technology 


SUMMARY 


Several proposed methods for swept-wing static structural 


lvsis are reviewed with the aim of consolidating the knowl 


inal 
edge currently available to the structural engineer faced with 


the problem of a swept-wing airct ift or missile design. In addi 
tion a fairly comprehensive bibliography on swept-wing analysi 
is included 


Three common types of construction——stiffened sheet, solid 


ind thick-walled shell structures—are discussed with respect to 
exact and approximate solutions for deflections and_ stresses 
The main analytical deficiencies lie in evaluating the effect of 


upport conditions, particularly for the thick-skin wing 


wing 


It is concluded that, if the structure is specified and is one 


of the above three types, it can in most cases be analyzed by 

plying or extending analysis methods now availabl Never 
theless, important work yet remains to be accomplished in the 
field of comparative design before the optimum swept-wing struc 
ture itself can be specified 


INTRODUCTION 


_— THE PAST SEVERAL YEARS, the practical 
importance of methods for the structural analysis 
of swept wings has been indicated by the growing 
While 
GALCIT has been primarily concerned with the funda 


number of papers relating to this problem 


mental problems in thick-skin and solid-wing analysis,t 
it has been necessary to examine and investigate to 
some extent the more prevalent stiffened sheet con- 
struction. On the basis of this combined experience 
it seemed desirable to review the published work in 
this field to date and attempt to point out the major 
analysis deficiencies that are believed still to exist. 
While the ensuing discussion is restricted to static 
malysis, it is impossible to avoid the dynamic effects 
that 


ture. In previous years, the structural engineer was 


are closely allied to the deformation of the struc 


primarily concerned with strength, whereas, at the 
present time, deflection criteria are becoming of such 
importance for the structures that operate in the re 
gion where swept wings might be used that, often, the 
importance of strength criteria becomes secondary. In 
some of the the most recent designs based on deforma 


tion design, the structure has been operating at as low 
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ing, I.A.S., New York, January 28 
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as 50 per cent of the permissible |.0g design load from a 
strength standpoint. Therefore, while the important 
field of dynamic analysis as such is omitted, this paper 
will consider static analysis methods for both deforma 
tion and strength with respect to three types of wing 
solid sections, and 


structures: (1) stiffened sheet, (2 


) 


3) thick-walled shell structures 


STIFFENED SHEET CONSTRUCTION 


One of the first authors to consider the problem of a 
skewed sheet element with stringers was S. Levy,! who 
in 1947 included such a case as one of his illustrative 
examples in a basic paper dealing with the application 
of energy methods to the finding of influence coefficients 
in stiffened sheet structures. In Levy's approximation 
of the energy in a skewed panel, it was assumed that the 
ixial component of the edge shear, when resolved to an 
equivalent rectangular panel, contributed to the sheet 
pointed out that this approximation 


energy Islinger 


violated the assumption that the thin sheet carrie« 


only shear stre sses: he proposed, for the « ise of sheet 
supported by stringers, that the axial force be divided 
between the stringers and, hence, contribute to the 
energy of the stringers. This modification was found to 
igree well with test data. 

Levy's method has been extended and applied by th 
McDonnell Aircraft Corporation in two reports by 
Islinger? and Goran and Goetsch*’ and by Bisplinghoff 
Tech 


Either method leads to 


ind Lang, ot the Massachusetts Institute of 
nology Aeroelastic Laboratory. 
sets of simultaneous equations which can be handled in 
a manner fairly satisfactory for practical applications 
to a specified structure. One of the more interesting 
results of the Lang-Bisplinghoff paper’ is shown in Figs 
| and 2, where the effect of sweep angle on tip defle« 
tion, rotation, and rear spar stress 1s indicated 
A second approach has been proposed by Bisch, 

who has set up the general equations for the consistent 
deformations of an arbitrarily shaped wing in terms of 
many tetrahedral panels. The solution hinges upon 


solving very large numbers of simultaneous equa 
tions by automatic calculating equipment and, in its 
present stage, is probably not suited for engineering 
application. 

Dixon”® has recently presented an analysis of swept 
and unswept wing structures which is based upon a 
the partial differential 


direct numerical solution of 


equations governing the elastic behavior of orthogonal 
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under a tip shear load and torque for the M.I.T. experimental 
wing. (See Fig. 18 of reference 5.) 
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shells. 
pends upon the number of elements in the idealized 
structure. In calculating stresses and deflections to 
compare with N.A.C.A. experimental data,* 


The degree of complexity of the solution de- 


Dixon 
applied a relaxation technique to the solution of the 
finite difference equations using desk calculators and 
considered a wing idealized to four, six, and ten spar 
caps or stringers. As would be expected, he found pro- 
gressively better agreement with the test data, some of 
While the re- 
sults were relatively easy to obtain using relaxation, 


which are reproduced in Figs. 3 and 4. 


especially for the four-element approximation, Dixon 


suggests the application of automatic calculating 
equipment. in actual engineering applications. 

Because the Dixon solution was based on the forinula- 
tion of the orthogonal shell equations given by Ben- 
scoter,? no cutouts in the structure were considered. 
The method itself—i.e., working directly from the dif- 
ferential equation—is not invalidated, however, if the 


equations contain the effect of cutouts. For the cases 


in which it can be used, the Dixon method is thought 
to offer a more efficient solution for a specified structure 
than either the energy or consistent deformatioy 
method. 

It will be noted that all of the three methods briefy 
discussed above have one common deficiency—they 
are analysis methods for a specified structure and thys 
have a relatively low degree of flexibility with respect 
to evaluating the effect of small changes in the strye. 
tural configuration upon the stresses and deflections 
For this reason it is desirable to have approximate 
methods available to evaluate the gross effects. One 
attempt in this direction has been proposed by D. B. 
Hall’ who has treated the wing outboard of the first 
rib as a cantilever, thus arbitrarily reducing the re. 
dundancy to the root triangle. D. B. Hall has empha- 
sized that his work was applied to a specific wing, and, 
hence, great care should be exercised in using this idea 
as a general approach for any arbitrary wing of any 
type of construction. Another simple analysis pro- 
cedure has been proposed by Bisch'! based upon a 
‘slice’ concept, in which the actual wing is repre- 
sented by a series of independent spanwise beams whose 
end conditions are determined by the rib or support 
restraint at these points. While the comparison with 
experimental data shows only fair agreement, the 
method may have applications in preliminary design 
before the wing structure is completely defined. 

There are some other experimental data available, 
besides those already mentioned, which may be used 
to check either approximate or exact theories. Zender 
and Libove*’ and Zender and Heldenfels'’ have re- 
ported N.A.C.A. stress and deflection data for a 45° 
swept box wing with a dog-leg; symmetrical and anti- 
symmetrical bending and torsion loading were used. 
Grossman and Forsberg'* obtained the stresses under 
symmetrical bending for a 35° half-scale swept wing. 

One of the conclusions of the latter investigators, 
when using various combinations of tip loads, was that 
they found no spanwise locus of points along which 
there was pure bending and no torsional deflections 
and vice versa. While this conclusion probably de- 
pends on the wing configuration, it is nevertheless 
interesting in view of a relatively common assumption 
that an effective length or apparent elastic axis exists. 

A somewhat diferent experimental approach has 
been taken by Redshaw and Palmer" and A. H. Hall,” 
who have tested plastic swept-wing models. The 
former paper is an enlightening exposition of the pit- 
The material actually 
used was Xylonite, a cellulose nitrate, while Hall's 


falls of plastic experimentation. 
work was based on a cellulose acetate. Two main dif- 
ficulties were in controlling the test conditions with 
respect to temperature and obtaining a satisfactory 
value of Young’s modulus and Poisson's ratio. Aside 
from the fact that plastic specimens may well prove 
to be a powerful experimental tool at comparatively 
low cost, certainly for qualitative investigations, the 








datz 
inte 
cha! 
vari 
func 
well 


wor 
buil 
add 
rib | 
tion 
veri 
exte 
pap 
stru 

I 
to | 
qua 
the 
(1) 
toa 
exa 
in ¢ 
cati 
at | 
she 
dev 


I 
and 
lish 
pla 
anc 
the 
apy 
As 
giv 
par 


Exa 


< 


con 
ten 
Th 
pat 
up 
fre 
det 
Fu 
cal 
ma 
Ne 
in 1 
Fig 
me 
in | 


Ought 
1cture 
lation 


riefly 
they 
| thus 
spect 
struc. 
tions, 
Mate 
One 
). B. 
first 
P te- 
pha- 
and, 
idea 
any 
pro- 
ma 
pre- 
hose 
Dort 
vith 
the 
sign 


ble, 
sed 
der 

re- 
45° 
iti- 
ed. 
der 











ANALYSIS 
data contained in these two reports are quite 
interesting. Redshaw and Palmer have found the 
changes in resonant mode position and frequency with 
variation in fuel load in the wings. Agreement of the 
fundamental mode and frequency for the model agreed 
well with the full-scale wing. A. H. Hall confined his 
work mainly to the static loading of the structure, a 45° 
built-up swept cantilever box section, and furnishes 
additional evidence (e.g., references 2 and 3) that wing 
rib orientation affects the stress and deflection distribu- 
tion very little. One of Hall’s main objectives was to 
verify the validity of plastic experimentation before 
extending his efforts to dynamic loadings. These two 
papers are strongly recommended to experimental 
structural engineers. 

In concluding this abbreviated discussion relating 
to stiffened sheet construction, it appears that ade- 
quate basic analysis methods exist for swept wings if 
the structure is specified. Two major deficiencies arise. 
(1) It is often difficult to apply a given analysis method 
toa practical structure, one with cutouts being a typical 
example. (2) A simple way of evaluating the changes 
in deflection and stress due to a small structural modifi- 
cation is needed. At the present time, it appears that, 
at least, methods of swept-wing analysis for stiffened 
sheet construction are now at about the same stage of 
development as those for unswept wing structures. 


SoLiIp SECTIONS 


In the field of solid wing sections of varying chord 
and thickness distributions, the majority of the pub- 
lished work has been concentrated upon the cantilever 
plate problem for a specimen having a constant chord 
and uniform thickness. As in most analysis problems, 
the solutions that have been tried vary from exact to 
approximate representations of the physical behavior. 
\ short discussion of the various methods will first be 
given and will then be followed by a quantitative com- 
parison of the various theories and experimental data. 


Exact Solutions 

Some of the first work was done by Fung’ " and 
concurrently by Young,®® with the former author ex- 
tending the formulation to the swept configurations. 
Their procedure, as reported earlier,‘! consisted of ex- 
panding the deflection function in a product series made 
up of a linear combination of the normal modes of free- 
free and clamped-free beams. The coefficients were 
determined following the Rayleigh-Ritz 
Fung concluded that the method was not practical be- 
cause of the slow rate of convergence of the approxi- 
mating series, especially at the larger sweep angles. 
Nevertheless, the calculated centerline deflections given 
in reference 17 for the case of an end load if plotted on 
Fig. 5 would agree almost exactly with the experi- 
mental data up to 40° of sweep. The chordwise slopes 
in the air-stream direction are, however, comparatively 
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stresses with the N. A. C. A. test data of reference 8 


low, and the stresses are probably in further disagree- 


ment 


A sector plan form has been investigated by the 
author?! following the suggestions of H. Wagner as 
reported by Zahorski.*® An infinite set of exact solu- 
tions, matched to a corresponding set of infinite char- 
acteristic loading conditions for a sector plate of uni- 
form thickness clamped along one radial edge, was ob- 
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Fic. 5. Tip air-stream rotation and tip deflection for a swept 


plate of constant thickness and chord ‘under a concentrated tip 
shear. 
tained. Unfortunately, from an engineering stand- 
point, the inverse problem is usually posed—.e., given 
some arbitrary loading, find the deflection and stresses. 
While it proved possible to obtain a representation of 
the deflection and stresses for this latter problem in 
terms of the infinite set of (nonorthogonal) character- 
istic deflection functions, the results, while of theoretical 
interest, are probably not easily enough applied to be 
of direct practical use. On the other hand, because 
each individual solution was exact (for the corre- 
sponding characteristic load), there was an interesting 
by-product believed to be of some practical importance 
in calculating stresses. There also seems to be some 
basis for extending the results to plates of variable 
thickness. 

The primary purpose in mentioning these re- 
sults at this point is to indicate the possibility of 
using this root-stress variation in conjunction with the 
Schuerch-Reissner-Stein results discussed below to 
save, in practical work, the necessity for a higher order 
approximation of the chordwise deflection. Some 
comparisons illustrating this point are given subse- 
quently. 

In another approach, which may possibly be termed 
an ‘“‘exact’’ solution (in the sense that an orderly method 
for obtaining successive approximations is developed), 
and a method that appears powerful for obtaining the 
deflection of cantilever plates, Schuerch*‘ derived from 
physical considerations a set of equations for plates of 
arbitrary thickness and plan form, subject to the 
conditions that the plate was thin and its middle plane 
deformed in such a manner that the chordwise de- 


| EXPTL(REF 18,20) 


flection varied linearly.* Reissner and Stein arriveq 
at the same equations in possibly a more elegant man. 
ner by first assuming the chordwise linearity and they 
determining the spanwise variation by using the varia. 
tional calculus upon the energy integral to obtain the 
governing total differential equations, which could, 
in many cases, be integrated directly. The Reissner 
and Stein approach can be extended with increasing 
ease to permit parabolic, cubic, or higher degree chord- 
wise representations, while the Schuerch method jg 
more amenable to considering a cubic variation directly, 
for example, without first carrying through the first. 
and second-order approximations. The latter physical 
approach is therefore of value if high-order approxima 
tion to the chordwise bending is required. 

For the cases in which the linear theory was com- 
pared with existing data for swept plates of constant 
chord and thickness, it appears that the linear approxi- 
mation gives an adequate representation of the deflec- 
with an order of magni- 
As will be illustrated 


tion up to approximately 40 
tude check on the root stresses. 
shortly, certain plan forms are more favorable to the 
method than others, and it seems likely that at least 
the parabolic term will have to be employed to obtain 
satisfactory results for deflections of variable thickness 
plates with sweptback trailing edges and even a cubic 
term if one wishes to obtain quantitative data on the 
stresses by this method.; This Schuerch-Reissner-Stein 
method (abbreviated to SRS in the ensuing compari 
son), with its flexibility in thickness, plan form, bound 
ary conditions, and application to dynamic loading, 
offers a most promising method for attacking swept 
plate structures. 

A numerical solution, and exact in the limiting sense, 
has been proposed by Fung, who has formulated the 
problem as one in relaxation techniques employing the 
ideas of Southwell.“© In the conventional biharmonic 
formulation of the cantilever plate problem, the de- 
flection function and its derivative are specified along 
the clamped boundary, and higher derivative require- 
ments (which cause the difficulty in the problem) are 
prescribed along the. free edges. The alternative 
formulation introduces two new “‘stress’’ functions, 
U and V, which assume specified boundary values on 
free edges and satisfy certain relations between their 
first derivatives on clamped edges. The relaxation 
technique is fairly rapid and can be applied to arbitrary 
plan forms and thicknesses. As in the case of Dixon's 
work for the stiffened sheet construction, it suffers the 
usual drawbacks of numerical methods but may be 
adapted to certain cases such as the dog-leg structure,” 
for example, which is not easily handled by the an- 
alytical method. 

* Schuerch’s formulation in reference 34 is not restricted to a 
homogeneous thickness distribution but is developed for hollow 
boxes and sandwich construction. The examples worked out in 
references 35 and 36 are, however, for solid hormogeneous 


plates. 
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ANALYSIS METHODS FOR 


Comparison of Theoretical Solutions 


As an indication of the degree of accuracy obtainable 
with present theories, three quantities have been 
selected for comparison with experimental data. These 
are the centerline deflection, the slope in the direction 
parallel to the support edge (which is designated the 
air-strearm rotation), and the stresses along the support 
edge (called the root stresses). It may be assumed 
that the degree of agreement in any of these three quan 
tities with the experimental data gives a measure of the 
overall acceptability of the theoretical prediction. The 
major comparisons have been assembled for two plan 
forms a swept rectangular plate of constant chord 
and thickness, and a delta plate. 

In Fig. 5, the tip deflection and rotation are shown 
as a function of sweep angle for a specimen of con- 
stant chord and thickness. As there is an aspect ratio 
effect, data are shown for two values of //c. Consider- 
ing Reissner-Stein’s statement that the assumption of 
linear chordwise deflection variation should be con 
sidered only a first approximation, the deflection is seen 
to be good up to approximately 40°; however, the ro- 
tations agree closely only up to 20 

As 60° of sweep is probably the highest value of inter- 
est, the deflections and rotations have been examined 
spanwise in Fig. 6. In this example, both the SRS 
linear) and Fung’s relaxation solution are low, with 
the latter, while yielding a considerably closer result 
still underestimating the deflection and slope by ap- 
proximately 20 per cent. The other curves on the 
graph are from approximate theories to be discussed 
presently. 

A further comparison of the SRS (linear) theory with 
experimental data from reference 23 is shown in Fig. 9 
for a delta plate having a varying trailing-edge angle. 
It is seen that, when the trailing edge is swept back, 
the theoretical and experimental data diverge. For this 
reason it is expected that at least a parabolic chordwise 
deflection will be required to represent the actual be- 
havior, 

One remark upon experimental technique may be 
mentioned at this point. It has been found that the 
deflection of ‘‘clamped ”’ specimens is strongly affected 
by the clamping pressure. As in the case of beams, it is 
difficult to match the theoretical boundary conditions 
that are imposed on the mid-plane of the plate with the 
experimental boundary conditions that are imposed 
on the surfaces of the plate. The two become more 
compatible as the thickness to chord ratio decreases. 
This conclusion is substantiated by Ashwell“ and Ger- 
ard* in their discussion of plate and beam rigidities, 
where it is found that plates develop their true rigidities 
only for large values of €(c/t)?. (€ is the outer fiber 
strains.) It will be noticed that the deflection data re- 
ported in references 18, 20, and 23 consistently lie 
above the theoretical predictions and, in the case of 
the l-in. specimen, give somewhat higher deflections 
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Fic. 6. Correlation of various centerline deflections and air 
stream rotations for a 60° swept plate under a tip shear load. 


than beam theory instead of (1 — yu’) less which would 
be expected from plate theory. It was furthermore 
found that, as the clamping pressure was increased or 
the plate thickness decreased, the deflections ap- 
proached the beam deflection values. The fact that, 
when the trailing edge is swept, the corner may yield 
also complicates the problem. As it may prove impos- 
sible to match the zero deflection and slope conditions 
exactly, theoretical data might be expected to have an 
inherent deficiency of up to 10 per cent of the experi- 
mental values, depending upon the thickness to chord 
ratio and the outer fiber strain. 


Approximate Solutions 


Having obtained some idea of the accuracy of the re- 
laxation solution and the range of validity of the linear 
chordwise deflection assumption by examining the 
extreme case of 60° swept rectangular plate, some ap- 
proximate theories will now be considered. The first 
stemmed from experimental evidence obtained by 
Martin and Gursahaney,** ** showing for various sweep 
angles the existence of an apparent elastic axis for uni- 
formly thick cantilever plates of delta and swept rec- 
tangular plan form. Recent unpublished data by 
Zender and Heldenfels, of the N.A.C.A., indicate that 
the position of the effective root varies not only with 
sweep angle but with aspect ratio. This latter conclu- 
sion is supported by the form of the theoretical solu- 
tion of Schuerch-Reissner-Stein. The main utility of 
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this concept lies in its application to high aspect ratio 
plates, inasmuch as the deflection in the vicinity of the 
root when using a straight elastic axis is undefined. 
Using the location of the apparent elastic axis deter- 
mined by a line parallel to the leading edge and inter- 
secting the root 18 per cent forward of the trailing 
edge—read from Martin and Gursahaney’s experi- 
mental curve (Fig. 1 of reference 43)—the deflection 
and air-stream rotations were calculated and plotted 
in Fig. 6. The results are seen to be in excellent agree- 
ment. As the method is so easy to employ, it will be 
particularly interesting to see whether or not an appar- 
ent elastic axis can also be found for variable thick- 
ness plates. The variation of its location with aspect 
ratio should also be further investigated. 


A second method of approximating the deflection of 
swept plates arose from an intuitive approach as- 
sociating the flow of stress with the flow of an incom- 
pressible fluid. That such an analogy may exist follows 
immediately from the two-dimensional similarity of the 
hydrodynamical continuity and irrotationality equa- 
tions to those governing the equilibrium of an elastic 
solid subjected to plane stress or strain when the latter 
equations are written in terms of the dilitation and ro- 
tation of anelement.* The connection between the ex- 

* See, for example, Hetényi, M., ‘‘On Similarities Between 
Stress and Flow Patterns,” Journal of Applied Physics, Vol. 12, 
p. 592, 1941, where existence theorems for the analogy for two- 
dimensional plane stress systems are discussed. 


tensional problems in plane elasticity to plate bending 
follows from the fact that the field equation is the same 

the Airy stress function and the deflection function 
are both biharmonic. This fact was first used in a prac. 
tical way by photoelasticians seeking experimental 
solutions to plate bending problems. Southwell has re- 
cently presented a compact association of the boundary 
conditions in the two cases. 

The basis for a consideration of the analogy stems 
from a somewhat different approach to the variables 
involved. Ordinarily, it is customary to consider two 
normal stresses and a shear stress on an element. In- 
deed, three quantities are required to determine the 
surface stress on a plate. Inasmuch as there are only 
two equilibrium equations, this means (in an approxi- 
mate sense neglecting compatibility) that some addi- 
tional condition must be used, such as minimizing the 
complementary or potential energy. Other common as- 
sumptions eliminate one of the variables by assump- 
tion, such as neglecting the cross stress (chordwise 
normal stress in the swept-wing case) for example. 
The alternative, which is suggested for certain cases, 
deals with three different quantities—the two prin- 
cipal stresses and their direction. An approximation is 
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then made to the principal stress direction. As a con- 
venient and intuitively correct way of defining the set 
of orthogonal lines of principal stress directions for cer- 
tain cases, the potential and stream lines of incom- 
pressible flow in a channel of the same shape as the 
plan form are chosen. The primary advantage gained 
as a direct consequence of dealing with principal stresses 
is that the shear stress vanishes on (curved) sections 
coinciding with the mesh, and, hence, the interaction 
of shear and normal stress is eliminated. The shear 
and direct bending stresses are thus separated in a 
convenient manner, and each individual stream line 
may be considered separately and treated as a curved 
beam. The assumption of this type of flow network is 
approximately equivalent to taking Poisson’s ratio 
equal to zero and, considering each stream line as a 
beam, neglects the normal stress between them. How- 
ever, the essential character of the curved stream lines 
which causes the elastic axis to move to the trailing edge 
has been preserved. Some additional comments and 
theoretical justification will be reserved until later when 
stresses are discussed. 

The center stream line, which, in the limit, checks the 
zero sweep case, is then chosen as an equivalent (curved) 
beam and the plate forces and moments established in 
relation to this curved beam, using the orthogonal 
network to give a spanwise moment distribution from 
which the deflection and rotation of the equivalent 
beam can be obtained. The results for the 60° swept 
plate of constant chord and thickness are shown in 
Fig. 6. These check the deflection but slightly over- 
estimate the rotation. A brief summary of the steps 
leading to the flow analogy solution is shown in outline 
form on Fig. 19. The coordinates of the center stream 
line for 45° and 60° swept rectangular plates and a 
30°-120° delta plate are given in Fig. 20. 

While not so clean a method as the SRS analytical 
presentation and somewhat more detailed than the ap- 
parent elastic axis of Martin and Gursahaney, it does 
give fairly accurate representations of stresses and de- 
flections, including the root region and at the higher 
sweep angles. The major drawback is the need for re- 
establishing a new orthogonal mesh for each loading 
condition, although for cantilever plates of constant 
thickness this is not believed too serious for reasonable 
loadings. Conceptually, the method is not restricted 
to rectangular plan forms* and may possibly be ex- 
tended to variable thickness and hollow plates and, 
as a final suggestion, may indicate the optimum ar- 
rangement of ribs and spars with respect to alignment 
in the stiffened sheet type of wing construction. It 
should be emphasized that, at present, it is suggested 
as a way of thinking and one possible means of at least 
semianalytical approach in the event difficulties occur 


7 . . . =* 
The symmetric, tip-loaded delta plate shown in Fig. 20, for 
example, may be treated as a diverging channel symmetrical 
about the root line, with a source and sink located at opposite 


vertices (points of load application ) 
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in other methods at the higher sweep angles. In con- 
clusion, it should be mentioned that it may be possible 
to modify the straight apparent elastic axis by curving 
it along one of the stream lines to give better deflection, 
and possibly even stress, results in the root region. 


Stress Formulas 


It is a well-known fact that, when the deflection is 
approximated, the stresses calculated by differentia- 
tion are in greater error than corresponding deflection. 
With one exception, the previous solutions have been 
obtained for deflections and slopes. As may be ex- 
pected, the stresses are somewhat in error. Inasmuch 
as deformations have recently been of considerable 
interest because of the aeroelastic implications, no 
particular emphasis has been placed on stress calcula- 
tions. Nevertheless, it seems desirable also to check 
their accuracy. 

Root stresses using the SRS (linear) theory have been 
calculated for a tip-loaded swept rectangular plate as a 
function of sweep angle and have been plotted in Fig. 8. 
In the notation of reference 37, the following formula 


applies where yu is Poisson's ratio, 
A = (l/c) V3(1 — pw)/2 
and the other quantities are obvious from Fig. 5. 


M,(o, vy) tan A 


P(l/c) l—y4p 


l y 
(1 — 4\ tanh 4A — )( — tan \) (1) 
cosh 4\/ \/ 
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Fig. 10. 
having a variable trailing-edge sweep under a tip load. 
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principal stress for a full-span, four-point-supported, variable 
thickness wing specimen under a uniform loading 


The case of a delta plate with swept trailing edge has 
also been considered, and the SRS (linear) results have 
been plotted and compared in Fig. 10 with experimental 
data.** For reference, the SRS formulas applicable 
to this case are shown on Fig. 9. 

The comparisons shown in Figs. 8 and 10 indicate a 
relatively poor correlation, and, aside from consider- 
ing a parabolic (or higher) variation in the chordwise 
deflection, a better method for approximating the 
stresses is indicated. 
mation to the stress field by writing stress equations, 


Fung” has given a direct approxi- 
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containing unknown coefficients, which satisfy the equa- 
tions of equilibrium of a swept rectangular plate 
and then determining the constants by applying the 
principle of minimum complementary energy.®  Thig 
procedure led to the stress rise curve given in reference 
17 or 41, which satisfactorily estimates the gross effects. 
On the other hand, it is known that, for plates in bend. 
ing, the stress tends to a mathematical infinity at the 
junction of a clamped and free edge if the enclosed 
angle exceeds approximately 90°.?! 
may be noted in Fig. 10. 


Such a tendency 
It seems reasonable, there- 
fore, that an estimate of the stress behavior should 
include this singular behavior. 

In the earlier discussion dealing with the deflection 
of sector plates, it was noted that an infinite set of exact 
solutions had been found. Any of these N solutions of 
the plate equation which satisfy the boundary condi- 
tion of zero deflection and slope along one radial edge 
6 = 0, and zero moment and shear along the other radial 
edge 0 = a, may be written in the form 


W,, if, 0) = a gue , 'F, (0; An(a) | (2? 


where is the distance from the vertex of the sector. 
Because the deflection and slope must be finite at the 
origin, A, > 0; the mome::t, which is directly propor- 
tional to the stress for constant thickness plates, is 
essentially the second derivative of the deflection, which 
means 


An 


a ~r"—* FO; Ana) 


It is then seen that, if any of the N values of \, are 
0 < A,(a) < 1, the stress can attain mathematically 


infinite values at the origin. The minimum values of 


Re d2(a) = ax(a) 


are plotted against the vertex angle, a, in Fig. 7, and 
. . . ~ . “- ~ 
it is seen that these infinite stresses can occur if a > 
90°. 


tion of the flow analogy arises. 


It is in this connection that further substantia- 
As can be seen in Fig. 
7, the curve of 7/(2a) lies close to the actual variation 
of d2 (a). The value 7/(2a), however, is the exponent 
of the conformal transformation of rectilinear flow to 
flow around a corner 2a. The velocity is thus propor- 
tional to 7 the 
r~~ Tt is then easily seen that the flow velocity 


A(n/2a)~1) Whereas stress varies as 
becomes infinite as a exceeds 90°, which is the same 
tendency present in the stress. 

The observation of this fact led to the postulation” 
of the following approximate expression for the stresses 
along the root line of a uniformly thick plate of fairly 
arbitrary plan form: 


6M T' [a2(ay) os 2( ae) | 


ol’, 0) = 
ci? Plaa(as)] Plas(as)] 


The notation is similar to Fig. 9, where 7 is measured 
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from the vertex of the angle, a; = A; + 90°, a2 = 90 
— Ae; Mo is the moment of the applied load on the 
plate about the root line; and a2(a;) is given in Fig. 7. 
The main restriction on the use of this expression seems 
to be that no large irregular loadings be imposed in the 
vicinity of the clamped edge. . The infinity at the origin 
may be dealt with by cutting off the curve at the yield 
stress. These latter two points are more fully discussed 
in reference 21. 

Eq. (3) 
swept rectangular plate and for a delta plate with 


was used to compute the root stress for a 


varying trailing-edge angle, and the results are shown 
in Figs. 8 and 10, respectively. In both cases it is 
seen that the SRS (linear) values are qualitatively cor- 
rect but do not give results so close to the correct dis- 
tribution as the product solution. With regard to the 
agreement with the experimental data of Hrebec*’ 
in Fig. 10, it should be noted that the product solution 
cannot change sign and, hence, cannot predict the 
negative stresses arising in the case of A; = 60°. In 
this sense, the formula gives a bound that is satisfactory 
only for single-signed stress variations. As the con- 
stant containing the Gamma functions is based upon 
equating the applied and reacting moments if the area 
under the moment curve is not made up entirely of 
positive contributions, an error proportional to the 
amount of reversed stress will result. For most load- 
ings, this error will be slight. 

In examining the experimental data, it is encourag- 
ing to see how the inclusion of the proper singular be- 
havior aids in representing the stress variation. In par- 
ticular, it should be noted that, for A; < 90°, the stress 


returns to zero at the origin. 


Variable Thickness Plates 


So far, the remarks upon solid plates have been con- 
fined to those of uniform thickness which, while of 
considerable theoretical interest, are not of particularly 
direct practical use. In this connection, the SRS solu 
tion should be most useful, because it is not restricted in 
thickness variation. On the other hand, it remains to 
be determined just what degree of approximation 
linear, parabolic, etc.—will be necessary to represent 
the plate behavior for various thickness distribu- 
tions. Some work along this line has already been car- 
ried out by Stein, Anderson, and Hedgepeth*® for delta 
plates of double-wedge cross section and unswept trail- 
ing edge. 

For comparison purposes, some of the N.A.C.A. 
data from reference 38 for the 45° delta has been re- 
produced in Figs. 12 and 13. In the first of these the 
plate is of uniform thickness, while in the second a 
double-wedge cross section of constant thickness ratio 
is considered. Inasmuch as the normalizing parameter 
in Fig. 13 is based on the average plate thickness, these 
two specimens are directly comparable on a weight 


basis. Regarding the deflections, it is seen, first of all, 
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plate of double-wedge cross section under a uniform load. (See 
reference 38.) 


that the effect of increasing the chordwise deflection 
approximation from linear to parabolic is small, al- 
though the effect on the rotations is somewhat more 
important. In the second place, the deflection of the 
double wedge is about two-thirds of the constant thick- 
ness section. Unfortunately, no test data are avail- 
able to check the validity of the theoretical results 
shown in Fig. 13, although a partial check may be 
made using some data recorded in reference 36. 

In the stresses, the change in approximation affects 
the distribution to a slightly greater degree, but in 
neither case does the approximation permit the stress to 
return to zero at the trailing edge. The pertinent fact 
to note in evaluating the effect of thickness variation 
is that the maximum stress for the double wedge is 
approximately one-third that of the constant thickness 
plate. 

In view of the possible change in character of the 
stress in variable thickness plates with trailing-edge 
sweep angle (which fortunately is not expected to be so 
serious as for plates of uniform thickness), it would be 
extremely valuable to make a check similar to that 
above for double wedges of the plan forms shown in 
Fig. 10. Considering the surprising success of the 
product solution for the uniform thickness case, a 
semiempirical extension of its applicability to variable 
thicknesses was suggested for qualitative work.” The 
bending moment about the root line, MM, (r, 0), as- 
suming the flexural rigidity, 
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(D = Et®/12(1 — p?) 


varied according to the cube of the thickness distribu- 
tion along the root line, was represented by 


MoV [m + ae(ar) + ne + de(ae) | 
Cc T [ny + 2(a;) | T'[ne of (ae) | 


r ni+a2(ai) —1 r n2+a2(a2) —1 
1 — (4) 
Cc Cc 


where, again, J, is the total moment of the applied 


M, (r, 0) = 


load about the root line. The values of 1, me are re- 
lated to the root thickness distribution in the vicinity 
of a, and as, respectively, such that 


th(r/c, 0) & 6 tar. (r/c)""/ 
te[(1 — r/c), O] & 6 tar. [1 — (r/c) ]" 
Having the relation that 
M,(r, 0) = (1/6)a¢ (r, O)é? (7, O) 
the stress at a point along the root is obtained by divid- 
ing Eq. (4) by one-sixth the square of the local thick- 
ness at the point. 

It was expected that Eq. (4) might not give so good 
results as Eq. (3), but, being a semiempirical relation, 
test data or other theoretical calculations were re- 
quired to test its validity. One set of experimental root 
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stress data* for a double-wedge section at zero Sweep 
angle*' indicated correlation with magnitude but not 
with position of maximum stress, while additional data® 
for the same specimen at 30° sweep indicated fairly 
satisfactory correlation in both magnitude and posi- 
tion. The distribution predicted by Eq. (4) was good 
only near the position of maximum stress, as shown in 
Fig. 14. 
tion has been superimposed on the N.A.C.A. data in Figs, 


For comparative purposes, the product solu- 


12 and 13 and is seen to give reasonably good agree- 
ment with the theoretical results. Experimental stress 
data for the double-wedge case would also be desirable. 

As a qualitative application of Eq. (4), the magnitude 
and location of maximum root stress for a cantilever 
double wedge of swept rectangular plan forms have been 
calculated, with the results shown in Fig. 14. It is 
important to note that double-wedge cross sections 
always have a finite maximum stress regardless of sweep 
angle. While the stresses for any thickness distribu- 
tion that can be represented as a power variation of r/ 
can be calculated, it is probable that Eq. (4) will yield 
closer results for the higher sweep angles. The main 
utility is expected to lie in its possibility of furnishing 
the approximate value of the maximum stresses in the 
plate or wing in the event proposed theories prove too 
unwieldy for rapid use in preliminary design. 

As a final example of dealing with swept wings ideal- 
ized as solid plates of variable thickness, Fung” has 
used relaxation techniques to obtain the stress distribu- 
tion in a 45° swept dog-leg full-span structure on a 
four-point support. The idealized wing has a constant 
thickness, /) from the leading edge back to the 50 per 
cent chord line and a linear taper to a thickness of fo/2 
at the trailing edge. The results calculated from the 
data in reference 24 are presented in Fig. 11 in terms of 
principal stress direction and magnitude. These 
numerical data indicate that the maximum stress occurs 
The stress at the trailing- 
edge corner is 40 per cent lower. 


at the rear reaction points. 
As this corner 15 
formed by the intersection of two unloaded radial 
edges, the stress is theoretically finite,‘ and, hence, 
the numerical value is representative of the actual state 
of stress. 

The relaxation technique employed by Fung is a 
numerical method that is not believed too cumberso:e 
for engineering investigation and is, indeed, believed to 
be the only method presently available, except for 
electrical analog methods described below, for handling 

* The data reported in references 31 and 32 were ol tained from 
a test specimen made from a gravity castable, thermosetting 
plastic, Paraplex P-43, manufactured by Rohm and Haas Com- 
pany, Philadelphia. 
over most thermosetting plastics because it can be bonded to 
itself without the use of a dissimilar bonding agent. The joint is 
as strong as the parent material. Rohm and Haas and GALCIT 
have been cooperating in an evaluation of the use of P-43 as a 


This material has a fabricating advantage 


structural test and/or photoelastic material. It shows promising 
characteristics for both purposes, and it is a pleasure to acknowl- 
edge the generous assistance and cooperation of the Rohm and 


Haas Company in this experimental investigation. 
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Fic. 15. Magnitude and location of the maximum stress 
along the root line of a swept rectangular plate of double-wedge 
cross section under a tip shear load. 


swept plates on isolated supports. Experimental data 
at GALCIT to check this example are only partially 
complete and are not ready for presentation at this time. 

One other method of solution should also be men- 
tioned, since it may prove the most versatile for rap- 
idly evaluating various loading conditions. It concerns 
the application of an electrical analog to plate prob- 
lems developed by MacNeal.'"® While the method is con- 
ceptually satisfactory for plates of fairly arbitrary 
thickness and plan form for both static and dynamic 
loadings, the investigation in reference 19 was re- 
stricted by the number of high-quality transformers 
For the plate specimens simulated, how- 
More recent and 


available. 
ever, the results were encouraging. 
extensive work has been carried out by the California 
Institute of Technology Analysis Laboratory on a delta 
wing. A comparison with experimental deflection 
tests made on a half-scale model showed favorable 
agreement for the static loadings investigated. Addi- 
tional work sponsored by the N.A.C.A. is under way 
for swept and delta wings of multicell construction. 


THICK-WALLED SHELL STRUCTURES 


In the arbitrary division of structural types, a thick- 
shelled structure was not defined. A precise definition 
would indeed practically solve the problem of which 
analysis method should be used for a given structure. 
At the present time, however, let it suffice to class a 
wing structure as thick-walled if the bending stresses 
through the skin thickness are believed sufficiently 
large to be included in the analysis. 

To the best of the author’s knowledge, three formu- 
lations have been published for thick-walled swept 
shell structures. The first, given by Wan,** is essen- 
tially an extension of Levy’s solution! for the stiffened 
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sheet construction. Wan includes the bending energy 
of the thick sheet and computes the influence coef- 
The results are 
As men- 


ficients using Castigliano’s theorem. 
theoretical, and no examples are calculated. 
tioned earlier, the formulation of Schuerch includes ap- 
plication to both thick-walled and sandwich structures, 
but no extended examples were calculated in refer- 
ence 34. The third solution is due to Fung” but has 
been applied only to a wing with cantilever boundary 
conditions. The hollow rectangular thick-walled sec- 
tion is analyzed as four plates rigidly joined together. 
The first step in the solution is to compute the stresses 
in the box by the conventional beam theory assuming 
the box to be clamped at a section normal to the spars. 
This fictitious root is taken norimal to and at the root 
of the leading edge for a sweptback wing. The bound- 
ary conditions that the slope and deflection at the root 
of the trailing edge be zero are thus violated. To satisfy 
all the boundary conditions, a system of self-equilibrat- 
ing forces and moments is applied at the root section 
(the beam is then free of constraint) and is chosen of 
such magnitude that it compensates the error introduced 
by the normal beam theory. The final stress state is 
obtained by a superposition. 

While it is believed the latter method admits of ex- 
tension, the most efficient way of handling different 
support conditions is not yet clear. The present formu- 
lation may, however, be adapted to multispar construc- 
tion containing isolated ribs. One of the most important 
stress effects that has appeared in the experimental 
work is the relatively high chordwise stress that occurs 
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for a tip shear load.” 
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Fic. 17. Comparison of various air-stream rotations and 
centerline deflections for a 40° thick-walled shell specimen show- 
ing the effect of internal structural changes. 4 
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(below), showing the effect of internal structural changes. 
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in the cover plate. In other words, theoretical soy. 
tions must consider bending rigidities in two orthogonal 
directions—not solely in the spanwise direction. As , 
typical illustration of the agreement with experimental 
data, Fig. 19 of reference 27 has been reproduced jp 
this paper as Fig. 16. It shows the maximum principal 
stress due to a tip shear load for a 40° and 60° swept 
wing along a section (designated the rib chord) per- 
pendicular to the leading edge and intersecting the 
trailing edge at the support. The basic specimen is q 
two-spar hollow rectangular shell with the cover plates 
‘40th of the rib chord. Two thickness ratios 


25 and 15 per cent overall 


being 
are shown again based on 
the rib chord. 

In concluding this section on heavy skin construc- 
tion, attention should be drawn to Figs. 17 and 18, 
which show some preliminary experimental data on the 
effect of internal construction changes that have been 
extracted from reference 32. As may be seen, neither 
the air-stream rotation nor the centerline deflection is 
particularly affected by the three changes considered. 
Furthermore, a comparison of these normalized rota- 
tion and deflection values at the tip with the similar 
values for the plate specimen of Fig. 5 (40-in. length 
and 40° sweep angle) shows good correlation. This 
similarity may be important in design applications. 
The principal stresses across the same section as in Fig. 
16 are also approximately the same in both distribution 
and magnitude, which is somewhat surprising, although 
the proper tendency is evident. One might conclude, 
however, that the stress rise is not reduced so much as 
might be expected by the addition of a spar or a rib 
perpendicular to the trailing edge at the root. 

It is interesting to note that this small effect of a rib 
was also noted by Grossman and Forsberg!* in the case 
of a swept wing of stiffened sheet construction. 

In evaluating the effect of internal structural changes 
or additions upon the basis of the data shown here 
for rectangular box sections, it is well to remember that, 
aside from adding material to take the load, there is a 
second alternative that in many applications to swept- 
wing design is more important—namely, the high stress 
can be avoided by putting less material at the point of 
concentration. Two of the possible solutions, for ex- 
ample, are to build essentially single-spar wings or re- 
duce the wing thickness proportional to the stress rise. 


CONCLUSION 


One of the obvious difficulties in writing a paper 
purporting to review a given subject is the sin of omis- 
sion, and the author regrets any methods that have not 
been discussed in their proper perspective. One reasom 
for the comparative brevity of the paper is the bib- 
liography and critical review of swept-wing literature 
prepared by Sechler and Dixon® in the summer of 1950, 
from which many of the unclassified references have 


been taken. The main intention of the author has been 
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to discuss various methods that, to him, seemed the 
more promising from the comparative basis, because, 
while it seems fairly safe to conclude that there are now 
enough basic analysis methods available to handle 
with engineering accuracy the stresses and deflections 
in stiffened sheet, solid, and thick-walled swept-wing 
structures, with the possible exception of thick-walled 
wings with arbitrary support conditions, much remains 
to be done in optimizing the structure for a given flight 
velocity. 

It is well known that sweep is a favorable aerody- 
namic device to employ at transonic speeds because of 
the drag reduction. On the other hand, after a certain 
amount of it on the wings of either aircraft or missiles, 
sweep becomes structurally inefficient from the weight 
and stiffriess standpoint, and it again becomes desir- 
At hypersonic speeds, 
Neverthe- 


able to build straight wings. 
wings as such may practically be omitted. 
less, in the transonic region it appears that an optimiz- 
ing problem should be set up, possibly as a variational 
problem with dual criteria—namely, that both the 
drag and weight be minimized simultaneously. One 
interim report has been written along these lines and 
an alternative formulation in terms of minimizing the 
overall work output has been prepared, but much re- 
mains to be done. The encouraging thing about the 
solid and thick-shelled structures is that they seem 
more amenable to analytical treatment, from the stand- 
point of investigating overall weight and stiffness varia- 
tions in relative skin and spar thickness for instance, 
than does the conventional structure. 

The foregoing review has been confined to three types 
of swept-wing structures, excluding sandwich con- 
struction, for example. To date, sandwich construc- 
tion has been used relatively little on aircraft, although 
the author believes it may prove extremely efficient, 
in the production and economic sense, for unrecover- 
able pilotless missiles. It is expected that, if a satis- 
factory method of comparing the structural efficiency of 
the three types discussed herein is obtained, the sand- 
wich-type construction will be bracketed and _ the 
proper region for its utility will be determined. 

In a sense then, a very real problem remains. 
if there are sufficient separate analysis methods avail- 


Even 


able to treat a specific type of construction, no analyti- 
cal method yet exists, except possibly one based on 
Schuerch’s work or the flow analogy, which permits a 
rapid evaluation of the relative efficiency of one type of 
construction against another for a given job. The 
proper compromise between the structural engineer's 
round tube and the aerodynamicist’s thin plate is still 
unresolved. 
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Fic. 19. A brief outline of the flow analogy method for cal 
culating the deflection of a uniformly thick swept plate under a 
tip shear load. > 

(1) After approximating the potential lines by straight lines, 
determine the deflection and twist of the center streamline by 
numerical integration of 

ElI(s) [d?W(s)/ds?] = M(s) = Pd,(s) (a) 
GJ(s) |d*y(s)/ds?| = T(s) = Pds) (b) 
using 
I(s) = (1/12)b(s)t8 
J(s) = (B/3)b(s)t8 
8 = B(b/t), torsion factor 


(2) Knowing W(s) and ys), the deflection elsewhere on the 
plate is calculated assuming linear variation of the deflection 
along the (straight ) potential lines. 

(3) Other loadings are handled by considering strips As wide 
and resolving forces and moments to the center stream line. In 
the case illustrated, the distributed tip shear is considered as a 
concentrated load at mid-chord 
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SUMMARY 





A method for developing cascade profiles for high subsonic 
potential flows to satisfy prescribed inlet and exit flow conditions 
and approximately prescribed cascade solidity is given. The 
method is based on the assumption that the compressibility of a 
gas can be represented by a linear pressure-volume relation (the 
Karman-Tsien approximation). It is shown that with the proper 
selection of this pressure-volume relation, the error in specific 
volume in the linearized gas as compared to an isentropic gas at 
the same velocity is less than 2 per cent of the isentropic stagna 
tion specific volume for Mach Numbers between approximately 
0.4 to 0.9. An example shows that this Mach Number range 
covers almost the entire length of the cascade profile for an inlet 
Mach Number of 0.7 and a deflection of 22°. The deviation of 
the fluid angle from the blade angle at exit and the influence of 
inlet Mach Numbers on deviation are found to be in substantial 
agreement with previous experimental results. 


NOMENCLATURE 


location of symmetrically situated complex sources 


a = 
a = [in Appendix (II) | local sound velocity 
a, @2 = location of complex sources 

ao = sound velocity at stagnation of linearized gas 
1,B = constants, strength of complex sources 
¢ = chord 

Cc = constant defined by Eq. (23) 

Ci, C2 = constants defined by Eq. (17) 

E = constant defined by Eq. (14) 

F = ¢ + ty, complex potential 

o(¢) = mapping function in Eq. (1) 

a 

J;,... = complex constants 

Ki, Ke, 

Ks, = complex constants 

k = ratio of specific heats 

is = leading edge 

M = Mach Number 

p = pressure, pitch 

P = singularity in mapping function 

q = velocity in the cascade plane 

s = singularity in mapping function 

ya = trailing edge 

v = specific volume 

w = transformed velocity, see Eq. (8) 

W = we *? 

z = cascade plane 
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The Development of Cascade Profiles for 


a = angle on the unit circle in ¢-plane 

B = see Fig. 3 

¥ = see Fig. | 

€ = profile angle at trailing edge 

¢ = unit circle plane 

7) = inclination of velocity in the cascade plaix 
p = density 

01, 0 = inlet and outlet flow angle, see Fig. | 

T = blade exit angle, see Fig. 6 


INTRODUCTION 


— THEORETICAL ANALYSIS of flow through cascades 
has been the topic of a considerable amount of jn 
vestigation. A recent survey” lists almost 50 articles 
on this subject. However, among the many papers i 
existence, few deal with compressible flows, so that 
although both the direct and the inverse problem of an 
incompressible cascade can be said as more or less 
solved, a theory of comparable scope for compressibk 
cascades is still lacking. 

That a theory of compressible flow in cascades is im 
portant is shown by the fact that many cascades of air 
foils in actual use in turbomachinery are approaching 
the sonic velocity, although with only a few exceptions 
they are still in the subsonic range. 

Lin! first suggested that, with the assumption of a 
linear pressure-volume relation for a gas, the problem 
of generating a cascade for compressible flows could be 
accomplished without too much difficulty. To carr 
out such a process, there can be two methods of aj 
proach. One approach is to try to solve the complet 
“inverse” problem—namely, the finding of a_ profil 
that will give a prescribed velocity or pressure distribu 
tion along the entire profile or, if this is impossible 
along a greater part of the profile. This approach 1s 
necessarily tedious, since it requires a great deal 
numerical or graphical integration. Furthermore, tt 
is by no means certain as to what is the most desirable 
velocity or pressure distribution. Carried to its logt- 
cal conclusion, this method eventually aims at de- 
veloping an ideal cascade profile.** 

A second approach, which is the one used in this 
paper, is to find the simplest possible mapping functiot 
that will generate a cascade profile with a reasonable 
velocity The labor 1s 
thus reduced to a minimum. 


distribution. computational 


tained by this method may not be the ideal one, still, by 


changing the location of singularities (see below), one 


** The complete ‘‘inverse’’ problem has been investigated by 


Costello, Sinnette, and Cummings.!'!~! 





Although the profile ob- 





x 


/ 
J 


can obt 
to appr 
The 
in com 
in dev 
in whic 
ence fo 
parison 
tion Cal 
the thre 


With 
tion (o! 
compre 
mation 
to a cor 
by the f 


t 
os 


In this | 
of the : 
(=a,4 
Fig. 1) 
a; and 

then th 





-ades 
f in 
ticles 
TS it 
that 
f an 


less | 


sible 


Fair 
hing 


tions 


this 
tio! 


able 


ob- 
, by 


one 











CASCADE PROFILES 





a 
= 
T s 
a r4 
\ £ 
a Ja \ ‘I x 
= - r 
\ a ‘ 6, ™y. 
5 an @ 
a - 
AS 
u\ ee es 
\ = ox T 
a ji— 7 6 
xX ? < ie: 
a a Y ~ ie 
¢€ a» 
C -— PLANE Z- PLANE 


NIT CIRCLE COMPRESSIBLE CASCADE 


Fic. 1. Mapping between cascade and unit circle 








¢-plane for one family of cascades 


can obtain a vast family of cascades that may be made 
toapproach very nearly the most desired one. 

The method given in this paper, because of its ease 
in computation and its flexibility, might prove helpful 
in developing profiles for high-speed cascade testing, 
in which a potential flow solution is desirable as a refe: 
ence for comparing actual test data, for only by com- 
parison with a known two-dimensional potential solu- 
tion can one evaluate the magnitude of the viscous and 
the three-dimensional effects of an actual cascade. 


ANALYSIS 


With the assumption of a linear pressure-volume rela 
tion (often called Karman-Tsien approximation) for a 
compressible fluid, Lin' has shown that the transfor- 
mation from an incompressible unit circle (¢-plane) 
toa compressible cascade (z-plane) can be accomplished 
by the following functional relation: 


1° == 


az = 2(C) (¢ — (1) 


(1/4a9?) F’? ¢) 2(¢) l (f—a,) ({ — Ge) dé (1) 
In this equation, F(¢) is the complex potential function 
of the incompressible flow of two complex sources at 
f= a, and ¢ = a, outside of the unit circle in the ¢-plane 


Fig. 1). If A is the strength of the complex source at 


a, and B that at a, (A and B are complex numbers), 


then the function F(¢) is represented by 


“OR SUBSONIC FLOWS 631 
F(¢) = A log (¢ — a) 4 A log [¢ — (1 ay) t 
B log ([ — do) + B log ¢ — (1/a) (2) 
Its derivative, or F’(¢) in Eq. (1), is 
‘ A 
F (Cc) = , T 
> — a c= (i/d; 
B B 
: t (3 
¢ de ¢ (1/ae) 


Since the strength of the source (or sink) at a; must be 
the same as the strength of the sink (or source) at do so 
that there is no net source or sink at ©, we must have 


Re A = —ReB (4) 


The function g(¢) in Eq. (1) must satisfy the follow 
ing four conditions: 


(a) Along the unit circle |¢) = 1, 


f dz = J (5) 
This is the condition for closure of the profile. 

(b) It is regular on and outside the unit circle, in 
cluding infinity. 

(c) It does not vanish on and outside the unit circle, 
including infinity, except at one point on the circle 
where, also, F’(¢) = 0; the order of the zero does not 
exceed 1. This point on the circle is mapped into a 
sharp trailing edge on the profile. 

(d) The following inequality holds for the region on, 
and outside of, the unit circle: 


(—1 


rmUvY/e ‘ee o. ra iui ? " 
F’(¢)/ao] \2(0) 457! (© — ai) (€ — ae)| <2 (6) 


The compressible velocity field in the cascade plane 
is completely given by 


i0 


W(t) = we = [F’(¢)/2(¢)] (€ — ay) (€ — ae) (7) 


where the ‘‘distorted”’ velocity w is related to the true 


velocity g by 


w = 2aog/(ado + V ao” + q?*) (5) 


For a complete background of this type of trans 
formation, the reader is referred to references | to 4. 

It is seen from the above discussion that, in order to 
complete the transformation represented by Eq. (1), 
our task is to find: (a) the strength of the two complex 
sources, or A, B; (b) the location of the two complex 
sources,.or @}, de; and (c) the function g(¢). 


Determination of A, B 


In many cascade problems, the fluid inlet and exit 
conditions (i.e., velocity or Mach Number, and direc- 
tion) are prescribed to satisfy other design requirements. 
Referring to Fig. 1, we are given the following quanti 
ties: oi, gi, o2, and g.. These quantities must, of 
course, satisfy among themselves the equation of con- 


tinuity, 


Gipi COS 0) = G2p2 COS a2 (9) 
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Fic. 3. Flat-plate cascade 
With gi, g2 given, we also have w, we, pi, and p» [see 
Appendix (II) |. 
Since the points ¢ = a, dz transform into s = — ©, 
+o respectively, the velocity at z = +o of the cas 
cade can be obtained by substituting ¢ = a, a. into Eq. 
(7), which becomes 
W_. = we “' = A [(a, — a2)/g(as)] (10) 
W4. = Wee *® = B [(az — a;)/g(ar) | (11) 
From these it can be shown | Appendix (1) | that 
Im A/Re A = —(po/p;) tan oy (12) 
Im B/Re B = —(po/p2) tan a» (13) 
Eqs. (12) and (13), together with Eq. (4), determine 
the complex source strengths A and B except to a com 
mon factor. This common factor can be obtained Fic. 4. Solidity for flat-plate cascade Numbers on cur 
finally by using either Eq. (10) or (11), but only after represent ar in degrees. ) 
we decide on the values for a), a2 and the function ¢(¢). ae ; 
It can be seen that, if we divide Eq. (10) by Eq. rhe right-hand side of this equation consists of know! 
quantities; hence, let it be called H. This equatior 


(11), we obtain i ‘ Rig : . 
therefore imposes a condition on the function g(f 
8\d2) wi eile F —_ #ho/p2) tan “| =E (14 It can be shown that Eg. (14) is identical with th 


g(a) Ws 1 — 2(po/p;) tan oy condition for closure or Eq. (5). 


Stagnation Points 


Let the leading and trailing stagnation points be located at ¢ = e'*” and ¢ = e'*’, respectively, on the unit cirel 
(Fig. 1). From F’(¢) = Oat these two points, we have, by Eq. (3), 
A A B : B 
ta 1] ta = + 1a 7 ta 0 7 
e —-aQ e — (1/d)) e — as e* — (1/az) 
a= a,,a,. Through some manipulation, Eq. (15) can be written in the following form, which has the advantage 
of involving only real quantities: 
sis (=<) : ‘ ) §1 + Ja)? , 
- SM AQ -, AL, 21) ¢ * COS (@ — AZ a2) = 
( 2a, Re A $2 ay f 
jl aol? (Im =) i \ pag d se ey . \ ' , 
. - sin (a Arg a2) 7 COs (a Arg a) 0 10 
( ia \ Re B ft la, ? f 


This equation can easily be solved for a with a few trials. One could also obtain an explicit solution for a from 
Eq. (15) by simply adding up the four fractions on the left side; the resulting numerator is a quadratic function 0 
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/ which can be explicitly solved. Such a procedure, however, requires much manipulation with complex quan 


tities and does not appear to be so simple as Eq. (16). 


symmetrically Situated Complex Sources 

From now on, to simplify calculations, we shall restrict ourselves to a case in which the two complex sources 
are symmetrically located with respect to the unit circle in the ¢-plane. For such a case there is no loss of gen 
erality if, by a rotation of the axis, we further make the two complex sources fall on the real axis—that is, a; - 
ete = a, witha being a real, positive quantity (Fig. 2). 

Bg. (16) can now be reduced to a simple quadratic equation for sin a, 


(1 + C,?) sin? a + 2C, sin a + (Ci? — C2?) = O (17 


» are known constants, 


| — “(= A Im = ) — 1 — a? (* A 4 Im = 
ReA ReB/’ ~ 21+ a?)\ReA° ReB 


where Ci, ¢ 


C 
la 


The two roots of Eq. (17) are a, and ap. 


Determination of a 

The value of a determines the solidity (chord /pitch ratio) of the cascade. However, the solidity is also influenced 
toasmaller degree by the mapping function g(¢). It is therefore impossible to predict exactly the solidity before 
we have the entire mapping function. Fortunately, in most cases in practice, the solidity, unlike the inlet and 
exit conditions, need only have some approximate values based on previous design experience. In the present 
theory, an approximate solidity can be predicted by comparing the compressible cascade with an ‘equivalent’ 
incompressible flat-plate cascade whose solidity can be exactly obtained. 

The transformation from a unit circle (¢-plane) to an incompressible flat-plate cascade (z-plane) is® (Fig. 


> 
> 


Db ad T c 


. ( 
g= - log + e” log IS 
2r a ¢ ¢ (l/a 
with the solidity (chord /pitch ratio) 
( | 1 + a’? + 2a cos az : 2a sin ar| 
cos 8 log + 2 sin 8 tan 19 
§ : 9 ; > ( 
p T 1 + a’ — 2a COS az a- — | 
and the location of a, 
tan a7 = tan B |(a° 1)/(a> + 1 20 


Eliminating 6 in Eq. (19) through Eq. (20), we obtain the solidity as a function of a and a@;. This is plotted in 
Fig. 4. 

The choice of a is therefore as follows: First assume a value of a. Find a7 from Eq. (17 
10°. Find solidity from Fig. 4. If the solidity is not what we desire, choose 


Experience shows 


that ay is usually between 4° to 
another a (smaller a corresponds to larger solidity and vice versa) and repeat the process until the desired a ts 
found 

Actual calculations show that, for ordinary compressor cascades, the solidity is only 70 to 90 per cent of the 
solidity of the incompressible flat-plate cascade with the same a, and that the higher the Mach Number, the lower 
the solidity of the actual cascade, again with the same a. With this knowledge we can come closer to the desir 
able solidity of the actual cascade by picking an a corresponding to a solidity of a flat plate cascade that is 15 to 


1) per cent greater than the desired solidity of the actual cascade. 


Function g(¢): One Family of Cascades 


A simple g(¢) function first suggested by Weinig* in connection with incompressible cascades is as follows: 


| 1 | ae Ss e'@ ] én P 2+(e/r 
g(¢) -(1 + — (1 : (1 ‘)( -) (1 2 
\ ag ag ( 4 4 


where —a and +a are the locations of the complex sources outside of the unit circle. 


5 is located on the radius to L (Fig. 2) at a distance of from 0.0125 to 0.15 from L. 
P is to be calculated in each case by Eq. (14), and ¢ is the profile 


Larger distance corre 


sponds to a more blunt nose, and vice versa.) 
angle enclosed by the trailing edge. 
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For this class of cascades, the mapping function is, after substituting Eqs. (21) and (3) into Eq. (1 


l €e/1) 9 





(¢ — S) (¢ — e’*7) Whe 
dz = : ; : = as “ee =e 4 — 
(¢ + a) (¢ — a) [§ + (1/a)] [F& — (1/a)] (¢ — P)* ; 
( ( :) (¢ et! \2 (¢ -_ e'%! 1 t+ (e/a (¢ vee P)* € 
fay” iC =F a) (E = 8) [¢ + (l/a)| |¢ — (1/a (¢ As OP a 7 ik: 
where 
C = (B — Ala - (B-=.A) 0 a) 93 
The calculation of the cascade profile from Eq. (22) K3J3; = B/C? 
requires, in general, a stepwise integration. For the 72 : 
ee we a Ji + Js + (c?/4ao?) (Ki + Kz) = 0 
special case of « = 0, corresponding to a cusp at the ; 
aa ae ale 7 Be: : ar + a, = 180° — 2 Arg (C/a 
trailing edge, one can easily integrate Eq. (22) as fol 
lows: Let Eq. (22) be rewritten as The mapping function can, for « 0, be integrated 4s 
C? oma (C2/4ao2)z 97 
dz = H,(¢) dt — ——. Hx(¢) dt (24) _—* mn “1 
rly 
where 
where ae aan 
% = J; log “a + @) + J» log ¢ + (1/a | 
IT, (¢) = Js log (¢ — a) 4 Js log |¢ 1/a)| 4 
(¢ — S) (€ — e@T)e? 5 Js log (¢ — P) — Jel€ P 
(25) 
l l : : 
(¢ + a) (¢ — a) (« ++ ¢(- (¢ — P)? Zo = K, log (¢ + a) + Ke log [¢ + (1/a 
a a Ks log (C - a + K,4 log ¢ : | a + 
H(t) = Ks log ¢ — Acf—' + K; log (¢ -— S$ 
(eg — @*L)? (¢ — @*T)\(t — P)? For the family of cascades represented by Eq. (21 
: : ; G5 ; : 7 : ag : 
| | ‘-O) the velocity at every point is given by Eq. (7), which 
Ga a —@) (+ 3 aG - ) (¢ » i now becomes, for any e, 
( a 7 
? a ; ' ou " (fp a e@Ly(¢e — 97 €/T (% P)? € 
are known functions of ¢. We now write en ee an - : 28 
a ay (Cc S 
J; Ie A 
IT(¢) = : + : ; + F T 
(¢ + a) i¢ + (1/a) (¢ — a) Procedure 
Js 4 Js 4 Jt A step-by-step procedure to generate compressibl 
[¢ — (l/a) (¢ — P) (¢ — P)* cascades with prescribed trailing-edge angle € and 
K K I prescribed flows at + ~ can now be outlined as follows 
AN \2 \3 “* ° ° , ° 
H(t) = af 4 Given inlet velocity g:/do, inlet angle o;, outlet angk 
b ‘iss l [ % ( )] — ege 4 ~ _— . 
SF @) f+ (1/a)] = a, and trailing-edge profile angle «. [The relation be 
Ky fl Ks 4 Ke 1 K; tween the Mach Number of an isentropic gas and 4 
[¢ — (1/a)] ¢ ¢? (¢ — §)  ofalinearized gas is given in Appendix (II). | 


where the constants J;, Js, J3, /; can be evaluated by’ ; (1) one PAE Hoe Pal Phe Sei ee, Wa Se ware mie 
; tinuity Eq. (9) and linearized gas relations of Appen 
J, = (+a) A), (¢ > —a) dix (II). 
(2) Find ImA/ReA, Im B/ReB from Eqs. (12 
and (15). 
and so on, and the constants J;, /s are determined by (3) Determine a from the approximate solidity de 
sired (discussed under the heading ‘‘Determination 0! 


Jo = [¢ + (1/a) |JX(8), (¢ > —l1/a) 


J; = (d/dt) [((¢ — P)*(¢)], (¢ > P) 
a’’). 


Je = (§ — P)? MS), (¢ — P) (4) Find ar, a, by Eq. (17). 
(5) Locate P in the mapping function g(¢) of Eq 


similarly for the constants A,, Ao,..., Az. The follow- ; ao 7 : 3 
sean sbar scema gts . (21). Substituting Eq. (21) into Eq. (14), we have 


ing identities can serve as checks to one’s calculations: 


-_ _— pP\2-(e/z) a oe _ ~sag l e/n 
J, = g(—a)/—2a ( a -) & . )( a — € ) ai 
Iz = g(+a)/2a =F —a-— 35/\-e -—¢" 
Jy = | 2 
KiJ, = A?/C? This can be explicitly solved for the only unknown P 











This 
limit 
8 


a] 


direc 
Q 
It 

appli 


anes 


Exan 

Tv 
meth 
com] 
let a 
The | 
€ 1s | 
for t] 
an in 


mari: 


ated " 


P 


which 


ssible 

and 
llows 
angle 
n be 


| g/a 


con- 


pen 


r de- 


yn of 


" 
a 


if 

















CASCADE PROFILES 





4 4 PROFILE NO. | 


Vi 
/ Lo FLUID INLET c,; 45° 
FLUID EXIT de 23° 
BLADE EXIT Tv 20.365° 
ot SOLIDITY c/P 1.305 
STAGGER 32.7° 
SEE FIG 6 FOR EXPLANATION 
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Fic. 5 Incompressible cascade 

6) Find |A, B by Eqs. (10) and (11). A and B 
are now completely determined, since we already have 
Im A/Re A and Im B/Re B. 

7) Check for the fulfillment of the inequality Eq. 
6). This inequality becomes, for the given family of 
2(¢ 
(C do) (¢ = e’@L) (¢ a Pim: \é “( 2 P)* é/n x 

(¢ —S)¢-| < 2 


This condition is, in general, satisfied when the flow is 
limited to subsonic. 

(8) The profile is now obtained by integrating Eq. 
22). For the special case of « = 0, the profile is given 
directly by Eq. (27). 

(9) The velocity distribution is given by Eq. (28). 

It might be mentioned that the above procedure also 
applies to the incompressible cascade, which becomes 


an especially simple case. 


Example 


Two cascade profiles have been calculated using the 
methods developed above. In order to form a basis for 
comparison, both profiles are designed on the same in- 
let angle ¢; = 45° and the same exit angle 0. = 23 
The deflection for both profiles is therefore 22°. Also, 
€is taken as zero for both profiles. Profile No. 1 is 
for the incompressible cascade and Profile No. 2 is for 
an inlet Mach Number of 0.7. The results can be sum- 


marized as follows: 


FOR 


SUBSONIC FLOWS 635 


Profile No. 1 (Fig. 5) 


Solidity (chord/pitch ratio) = 1.305 
Blade angle at exit = 20.365° 
Deviation (fluid exit angle minus blade angle at exit) = 
2.635 

Stagger = 32.7 

Im A Im B 

= 1.0000, = ().49447 
Re A Re B 
ar = 4.00°, az = 189.53 


P = 0.113764 0.080675: 


Profile No. 2 ( Fig. 6) 


Solidity = 1.086 
Blade angle at exit = 20.056 
Deviation = 2.944 
Stagger = 32.1 

Im A Im B 

= 1.37144, = ().49899 
Re A Re 
aT = 4.677°, az = 193.024 


P = 0.123877 0.080825: 


For both profiles, a = 1.183, and S = 0.96 e'*”. It 
can be noted that, with the same a, the solidity of the 
compressible cascade is considerably smaller than that 
of the incompressible cascade. One could, if one de 
sires, bring the compressible profile nearer in solidity 
with the incompressible one by using a smaller value of 
a for the former. 

The velocity distribution around the profile is shown 


in Figs. 7 and 8. 
CONCLUSIONS AND SUGGESTIONS FOR FUTURE WORK 


Although no definite conclusion can be drawn when 
only two cascades have been calculated so far, there 
are several interesting aspects worthy of mentioning. 
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Fic. 7. Velocity distribution Profile No. | 


The deviation, or the difference of the fluid angle 
from the blade angle at exit, is around 3° in both cases. 
Howell’s* experimental data indicate, for this particular 
The 
slight discrepancy between theory and experiments on 


cascade arrangement, a deviation of about 5 


the value of deviation may be partly due to the vis 
cous and the three-dimensional effects in experiments 
and partly due to the somewhat unconventional upturn 
of the trailing edge of the developed profile. 


For the same fluid exit angle, the blade exit angle of 
the compressible cascade is only 0.309° larger. It 
appears therefore that, for a given cascade and inlet 
fluid angle, the influence of Mach Number on the fluid 
exit direction is negligible. This is also the result of ex- 
periments by Bogdonoff.® 

The velocity distribution does not seem out of line 
with those of actual blades used in practice. 


Future work could extend the use of this method to 
Mach 
There is also great flexibility in modifying the blade 


different .flow angles, solidity, and Numbers. 
shape and therefore in improving the velocity distribu- 
tion by varying the location of singularities in the func- 
Specifically, one can consider the following 
(1) Shift the location of Sin Eq. (21) 
(Q/¢) |" 


in Eq. (21), where Q is slightly away from the center 


tion g(¢). 
three variations: 
by aslight amount; (2) add another factor |1 - 


of the unit circle and the exponent n is preferably + | 
or +2 to facilitate computation ; (3) for greater gener 
ality it might be desirable to have the two complex 
sources located nonsyminetrically. Other variations of 
g(¢) function are, of course, also possible. 


APPENDIX (1) 


In the cascade plane in Fig. 1, the fluid angles at 


to, or 6, O, are given by Eqs. (10) and (11) or by 
4, Arg A(a, do) £(a, 
A, -Arg |Bl(ao — a))/g(ae 


SCIENCESSEPTEMBER, 


1952 


The angle of the cascade axis, or y, is, by the Caye) 
integral formula and by virtue of Eqs. (1) and (3 


cs 
7 Arg [¥ (a) dz] = 5 + Arg X 
g(a) 1 A*(a, — ay») : 
+ : \rg F ia. dz 
ay ad» lay g(a) 

7 g(a») | Ba 

- + Ag t 

Z ay = Ar hay ela 
where $,,, and §$,,) represent a small closed py 


around the singularities a, and de, respectively. 
If the fluid angles are measured from the casea 
axis, then (Fig. 1) 


Oj 6; + (7 2) Y 


a2 0. + (m/2) a 


Substituting the expressions for 6,, #2, and ¥ 


| 
0; = Arg 4 T , 
ay 


| Bla as)|* 
0» Arg B B 
by" Y\ dy» 








However, from Eqs. (10) and (11 
A (a, (le Bla» (1) 
Wy > (a 
2( a2) Y((» 
Hence, 
l (27/4a 9?) (Im A 
—tan o; ‘ ; ( 
1 + (w,"/4ao° Re A 
l (W»"/4ay?) (Im B 
tan oa» F , 
1 + (we? 4a Re B 
> q=VELOCITy 
° Q,=STAGNATION SOUND VELOCITY 
L247 
UPPER 
SURF ACE 
[0 4 
\ 
8 4 \ 


LOWER 


SURFACE 
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Fic. 8 Velocity distribution Profile No. 2 
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TABLE | 





Linearized Gas 


Error in 














MM. q/a q/a uv /r% w/a v/Vo Us /U 
() 0 0) | 0 0.93178 1. 00000 0. 06822 
(1 0.14085 0. 13899 1 OO980 0.13967 0.94091 1. 00500 0. 06409 
ds 0.2 0). 27987 0.26951 1.088438 6.27460 0.96759 1.02013 0.05254 
sa 0.3 0.41774 0. 38546 1. O8375 0.40095 1.00982 1.04561 0.03579 
(4 0.55319 0.48405 1.1428 0.51632 1.06484 1. O8192 0. O1711 
‘a, = all 05 0.68553 0. 56542 1.2124 0.61972 1.12973 1.12973 0 
. i, 0.6 0. 81414 0.63135 1.2895 0.71119 1.201538 1. 18984 +0.01169 
(a 0.7 0). 93854 0 68437 1.3715 0.79152 1.27794 1. 26330 +0.01464 
Rds 0.8 1. 05826 0.72685 1.4457 0. 86540 1.34710 1.351385 0.00425 
| 0.9 1.17299 0.76098 1.5414 0.92311 1.43625 1.45552 0.01927 
osed pat 1.0 1.28252 0.78861 1.6263 0.97667 1.515385 1.57744 0.06209 
C Casca ¥ " . F f , 
From the relations of a linearized gas, q/ay $(w/ado)/|4 (w/do) (II-4) 
v p 1 — (w*/4ao") In the above equations, g is the velocity, a is the ve 
? po 1 + (w?/4ay?) locity of sound, and the subscript 0 refers to the stag 
nation state. 

Hence, It has been customary in dealing with single profiles 
tan oy (p1/po) (Im A/Re A) to use the state of the free stream—i.e., the undis 
tan o» (p2po) (Im B/Re B) turbed flow at + o~—-as the reference state p), 7 in 

Eq. (II-1). In an infinite cascade, however, there is 
‘hare Eqs. (12) and (13). é : ; k 

i which ‘ no such thing as a state of undisturbed flow. For best 
approximation for high subsonic cascades, it is suggested 
APPENDIX (IT) here to pick the point corresponding to isentropic Mach 
Number (/,) of 0.5 as the reference point. Let the sub 
Relations of a Linearized Gas for Cascades script s denote the properties of an isentropic gas, as 
For a linearized gas satisfying the following pressure distinguished wom those bas ” linearized gas which use 

volun relation no subscript. Then the isentropic relation, 

Pi p Div — Vv) (II-1) (p Pos) (v./¥o.)* = 1 

where D is a constant and the subscript | represents any js being approximated by a chord (Fig. 9) passing 

chosen reference point. It can then be shown that the through M, = 0.5, with the slope equal to the slope of 

lollowing relations hold : another line joining the two points J, = O and /, = | 
> , 9 o > 1c °° RP _ 7s ( . P x. ; 
g?/a 1 — (o92/p*) (11-2) on the isentropic curve. Fig. 9 shows that, with this 
straight-line approximation, the error introduced is 

4 2 J 2 ) . . - e ° 
av/g? = (ao°/q°?) + I (II-3) kept to a minimum for the range of /, covering almost 
a, - : the entire length of a calculated cascade profile for high 
w/a 2q/ay)/{1 + V1 + (g/d) : : “ 

subsonic flows. 
or At the reference state (1), WW. = q:/ay, = 0.5 and 
Pi/ pos = 0.84302, 01/09, = 1.12973 when the ratio of spe 

——— — — —___—-— — cific heats is taken as 1.4. The slope of the line join 

> ; 

“Pre MACH NO ing the two points /, = O (where p,/ po, = v5/to, = 1 

Ms > os : ——— 

105 <0 and M, = 1 (where p,/ pos = 0.52828, v,/v0, = 1.57744 

\ is —0.81692. The linearized gas is therefore repre 
oh 3 
% sented by 
94 A 
\ - =A =O) __ 5) /e) /e 
\ nerenence nr p/ pos = 1.76592 O.S16092(v/ Vo, 
\ 

8 ‘NS The velocity of sound is defined as 

> 6’y ; 
NO (suoee = p1ee2) . ap ( v ) : a(P/ Py 
hes a~ v ry ) . 
74 ” dv Vos d(v/v 
ISENTROPIC . , , 
(K=1.4) rN For the isentropic gas, 
6 ee 
e t. : kp. bp, 
a ( VosP. ( ; 
1.0 9 ' 
= . Vos v,/t 
100 ai y T no T T — _ . : , 
10 1 1.2 1.3 1.4 5 1.6 For the linearized gas, 
SPECIFIC. VOLUME VA Vos 
a” (v/tos)" VosPos (O.S1692 


Fic. 9. Approximation by a linearized gas 
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At the reference point (1), v1 = v,; hence, we have, 
by taking the ratio of the two equations above and 
substituting the known value of p,/po, and v,/v, at 1, 


Q1/Q,; = 0.88429 


Since M,, = g:/a:, = 0.5, therefore, 


n/a = (a4, a) (M1 Qs) = 0.56543 


Eq. (II-3) gives g:/do = 0.68553. 

A table for this linearized gas can be compiled as 
follows: For a given isentropic Mach Number J/, be- 
tween 0 and 1, we can find the ratio of the velocity at 
that VW, to the velocity at W/, = 0.5. 
Then, 


This is g/q. 


q/do = (4/1) (Gi/do) = 0.68553 (g/qi) 
Having obtained g/d, we can find g/a by Eq. (II-5), 
w/ay by Eq. (II-4), and v/v by Eq. (II-2). The re 
sulting table is given herein as Table 1. 

In Table 1, the last three columns are for estimating 
the error involved in the linear approximation. At our 
reference point (1) where M, = 0.5, we have v/v = 
1.2124. 
ard gas tables. 


At the same point, v/v, = 1.1297 from stand- 


Hence, for this linearized gas, 


Uo/Vos = 1.1297/1.2124 = 0.93178 


Using uv, as the basis of comparison, we can express 
v of the linearized gas by the ratio 


o)/9 
U/U0s — 


0.93 l 7Sv ‘Uo 


Comparing this ratio with the isentropic v,/v, from 


standard gas tables, we find the error in v/v; shown in 
the last column of Table 1.* We see that the error is 
small for 17, between 0.4 and 0.9. 

* This error is an indication of the slight difference in com- 
pressibility of the linearized gas as compared with that of the 
isentropic gas for the same velocity. 
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RIEF REPORTS of investigations in the aeronautical sciences and discussions of papers pub 


lished in the JOURNAL are presented in this special department. 
completed as soon as possible after receipt of the material. 
itself responsible for the opinions expressed by the correspondents 


The publication will be 
The Editorial Committee does not hold 
Contributions should not exceed 


800 words in length 


Comments on ‘‘Remarks on Dynamic Loads in 
Landing’’' 


H. Norman Abramson 
Analytical Engineer, Dynamic Analysis Group, Chance Vought 
Aircraft Division, United Aircraft Corporation, Dallas, Tex 


May 9, 1952 


PROBLEM DISCUSSED by Mr. Zahorski! has received a 


_ 
great deal of attention in recent years; 


ysis of the problem is not yet in a completely satisfactory state 


however, the anal 


The writer would like to point out that the method of ‘‘normal 
modes” used in the subject paper has several limitations imposed 
on it which tend to hamper its accuracy. Among these the fol 
lowing may be specifically mentioned : 

(1) The duration of impact must be large compared with the 


period of the lower modes of vibration so that only a reasonable 


number of modes need be considered. This was indirectly 
brought out in the classic paper dealing with the landing impact 
problem by Biot and Bisplinghoff.2 It is shown that the mo 
ment and shear of a simply supported rectangular beam sub 


jected to a unit impulse is given very poorly by the normal mode 
procedure, while the same quantities resulting from a unit step 
function are given relatively more accurately. 

2) Obviously, the distribution of the loading on the struc 
ture will be of great importance, for, if the distribution were of the 
shape of one of the normal modes, only that single mode would 


be ex 





ited. In the problem under discussion, we normally have 
the load applied at two discrete points along the wing span, and 


If the 


landing is an unsymmetrical one, then we have even more modes 


therefore we will expect that many modes will be excited. 


to consider, since the unsyminetrical modes will also be excited. 

3) Finally, one must be aware of the fact that resonance be 
tween slightly coupled components of the structure may give rise 
to a transfer of energy from one component to another. 

+) As the author has pointed out, the initial shape of the 
forcing function is also of considerable importance. This is, of 
intimately related to item (1). 


course, 


In an experimental investigation, Ramberg* has compared the 


mode method (using the analysis of Biot and Bispling- 





hoff but 


modified to account for phase relationships between 
des) with two other methods It was found that both of 
the other methods gave more accurate results than the normal 


mode method for the case of unsymmetrical impact of large dura 
tion compared with the fundamental flexural period of the model 


wing under test. Foran unsymmetrical impact of short duration, 


one of the methods gave particularly good results; however, 
gain, the methods of references 4 and 5 gave results that wert 
etter than the normal mode method. In all fairness it should b« 
mentioned that the analysis by normal modes employed only 
three unsymmetrical mode 


O50 


Robinson® has also felt that . the method of normal modes 
may be unable to describe contemporary impact load problems 
adequately”’ (with regard to landing impact). Perhaps his anal- 
ysis of beam problems will lead the way to a traveling-wave-type 
analysis of complex structures under impulsive loading 

In connection with this discussion, one more paper should be 
mentioned. Kroll and Levy’ have taken the paper dealing with 
gust response by Houbolt* and adopted his suggestion that a 
similar analysis might prove of value in the landing problem but 
with some modification. Kroll and Levy studied the experi- 
ments of Ramberg’ and found, for the long-duration impact, a 
good correlation with the experimental results. The time re- 
quired for analysis by this method was slightly less than for the 
other methods 

This writer feels that the method of normal modes may be used 
with good results if the modes used (both number and shape) are 
carefully selected and if the loading situation is not severe |se¢ 
items (1) and (2) above 

Finally, the writer would like to refer to the question of treating 
the landing impact force as independent of the response of the 


Mr. Zahorski feels that information obtained 


structure 


from an arbitrary free drop of the landing gear is practically 
useless so far as the dynamic loads are concerned.”’ The writer 
concurs with this statement as it stands. However, he would 
like to point out, in view of several references,® that a good 


analysis is still possible when assuming the landing impact force 
to be independent of the airplane structural flexibility, provided 
that the force-time relations for the shock strut are obtained by 
replacing the flexible structure by a rigid one of the same mass 
Then, perhaps, we can use drop tests, but for the complete as- 


sembly and not just the landing gear itself 
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A Third Kind of ‘‘Incompressible’”’ Fluid* 


A. Martinot-Lagarde 
Director, Institute of Fluid Mechanics, Lille, France 
May 16, 1952 


gestae KOVASZNAY AND CoORRSIN! have considered the 


appropriate approximation to the heat term of the energy 
equation for two kinds of ‘‘incompressible’’ flow: (a) flow of a 
“perfect liquid”’ (density, » = constant, exactly) and (b) flow of 
a gas with Mach Number WV <1. I should like to examine 
the question in slightly different terms and to point out a third 
case: (c) flow of most real liquids. 

First, I should like to propose the term ‘‘isovolume fluid’’ (mete 
orologists often use ‘‘isopyenic’’) instead of “incompressible 
fluid’; the latter expression overemphasizes the pressure at the 
expense of other thermodynamic variables. 

In purely mechanical problems, we need only write that, for 
nearly isovolume fluids, dV == d(1/p) is zero or negligible. In 
thermodynamics, things are more complicated. For studying the 
approximation in question it is convenient to choose alternative 
pairs of independent variables. For example, with T and V, 
the heat received is 

dQ =c¢,dT +/1dV Sc, dT (1) 
while with 7’ and p, 


dQ =c,d7T +h dp =~ Cp dT (2) 


These are the two “‘isovolume’’ approximations in question, 
and we shall see that in most cases Eq. (2) is closer than Eq. (1) 
cp and ¢, are the customary specific heats, / (OU/OV)r + Pp, 
h==1(OV/dp)r, and U is internal energy. 

Let us denote the errors in Eqs. (1) and (2) by 


e’ = 1dV/c, dT; e = hdp/c, dT (3) 

and represent the equation of state over a small range by 
dV/V = adT —xdp (4) 
where a = (1/V) (OV/O7), and « = —(1/V) (OV/dp)r, the 


isothermal compressibility. 
From the two laws of thermodynamics, we can deduce some re 
sults of Clapeyron. 


h = —TVa; l= ; Coma = (5) 


These will permit numerical estimates of «and e’. We have first 
e = -—(TVa Cp) (dp dT) (6) 


As a typical real liquid, consider water near the boiling point 
flowing through a radiator: a@ = 0.00073, V = 0.00105, « = 
0.5 KX 107%, and c, = 4,210, all in m.k.s.a. units. For conserva 
tism in estimating the negligibility of «, take Ap relatively large 
and AJ small: Ap = 10° (about 1 atmosphere) and AJ = 10 
C. These give —e = +0.0007 < 1.0. 

We pursue ¢ a bit farther. If Ap corresponds to a velocity 
change from 0 to u, 

VAp = (1/2)u? 
whence, 
e = Tau?/c,AT ( 
For a gas €« = u?/2c, AT; for a liquid it is much smaller. 

In general, € can be interpreted in terms of temperature changes 
alone. If (d7), denotes the temperature change that would 
occur if the pressure variation dp were produced isentropically, 
(dT); = —hdp/cy, so that 

«= (dT ),/dT (8) 

* Adapted from a note in the Revue Générale de Mécanique, Vol. 40, p. 110 


1952 


Now we turn to e’. With the equation of state [Eq. (4 
lin Eq. (5), 
e’ = (TVa/kxc, dT) (a dT x dp) 


Although dV/V = (adT — «x dp) is small in isovolume floy 
multiplying factor may be large. Introducing a and x 4 
Eq. (5), Eq. (9) becomes 
o&—ec h dp 
<« = + a ety et) + ye 

é ae | 
where y = ¢p/¢y. The second term becomes negligible relatiy 
io unity as ¢ does, but the first term is negligible only in cer 
cases. For example, for water at 4°C., a = 0. However ¢. -| 


c, and approximations (1) and (2) are identical in this case any 


way. 

For gases, (y — 1) is not < 1.0. Consider a typical liqui 
For water at 100°C., y = 1.11; for mercury at 0°C., y = 113 
for ethyl ether, y = 1.40. It is clear that for gases and mo 


liquids, the c, dT approximation is permissible, while the , 
is not. 
Professor Corrsin has suggested to me that the e'/e ratio 


nits simple presentation. For example, 


é ldV ¢, y dV Kact 
€ hdp « «kV dp k 
where kact, is the compressibility in the actual process 1/I 


(dV/dp)| and x, is the ‘isentropic compressibility” («/; 





Eq. (11) shows that approximation (2) is better than (1) for «| 
process far from isentropic but that both are bad-—and roughly 
the same—for nearly isentropic processes 
The principal conclusion of this note is that dQ = c, dT whe 
ever |(d7),| < |d7|, even though it may seem paradoxical t 
equate actual heat received to constant pressure heat received j 


a process with varying pressure. 


REFERENCE 


' Corrsin, S., and Kovasznay, L.S. G., The Energy Equation for Two Ki 
f “‘Incompressible’’ Flow, Readers’ Forum, Journal of the Aeronautica 
Sciences, Vol. 18, No. 12, p. 843, December, 1951 


Computation of Possio’s Integral for Linearized 
Supersonic Flow 
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| {HE POTENTIAL,! ¢(x, y, + 0, t) for harmonic motions, 


x be 


de! = e — (vM/B?a)(x ©) dé wit, n) X 
Br fi J0 


4 
cos (x £) sin 6 | dé 
Ba 


is well known.?~‘ In Eq. (1), @ is defined? by 
n = [(x — &)/B] cos6+ y - 


The pressure on the +0 side of the wing is 


OP ly ; 
p(x, y, +0, the’! = ov ( -— ) e” 
Ox J 


The regions of wing plan forms for which Eq. (1) is applicat 
is, in part at least, well known. Evvard’s theorem* ! and 
theorem due to Miles and Froelich®:* are useful in this fiel 
The purpose of this note, however, is to outline a method wit! 
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which, by means of high-speed digital computing machines, the 
essure distribution and the flutter aerodynamic coefficients due 


pr a ° ° 
to Eq. (1), where applicable, may be computed for suitably arbi 
trary downwash functions, w(£, 7). 

By definition, the flutter aerodynamic coefficient is 


Qasp = { { Za (x, ¥) pa(x, y) dx dy (4) 
Jos 


where Za(x, y) is the deflection form a and pg(x, y) is the pres 
sure due to the deflection form Z(x, y) and the integration is over 
the top and bottom surfaces of the wing. On substituting Eq 
3) in Eq. (4) and integrating by parts, 


‘ - ole ly 
Zap t = 


Qas = eV Za) 3 dx | dy 


} Ox V 


rl / 


(5) 


where x, and x, are the values of x at the leading and trailing 


edges, respectively. vanishes at the leading edge, the limit 


x, is shown for completeness. 


Also, 
OZ iv, 
w(ix,y) = | te — 2 (6) 
Ox J 
On substituting Eq. (6) in Eq. (1), 
I us ’ 
6 = e— Wy M/Ba(z—§) de (~ + z) x 
Br St J0 Ox | 
v . 4 
cos (x £) sin 07 dé (7) 
\ 62a { 


The plan is to compute the potential © according to Eq. (7) at 
a set of points on the wing; then by means of these values the 
pressure distribution according to Eq. (3) and the flutter aero- 
dynamic coefficients according to Eq. (5) will be computed. In 
these computations and results, the deflections zg and 2g will re 
main as explicit symbols, which may be assigned values as re 
quired in the use of the results 

For the computations, let the points rs be on a rectangular 
lattice, and at the points of this lattice on the wing let the de 
Define the wing plan form by its 


flections sa, zg be known 


intersections with the lines of this lattice. A rectangular lattice 
is chosen in order that the digital computing machine programs 
may be arranged in part, once and for all, for all plan forms to 
which this method applies. Devise interpolation formulas so 
that, with functions known at the points rs, values of the func 
tions may be computed at intermediate points. Intermediate 
points will be introduced as required for the numerical integra 
tions 


In the following, let | |, [°|, | {, | | denote rectangular, diagonal, 


column, and row matrices, respectively. Construct a matrix 
w| defined by 
joZ lv ' ee 
T = = [w]iZrs} (8) 
lox ye, 


Consider now the computation of the potential at the point 


rson the wing. First, cousider computation of the integral 
“(oZ iv ,) fv ig 
Z > cos (x £) sin a 


+ 10 
4 {dx vo§ | 82a , 


(9) 


for a given point rs at vy and for a particular within the forward 
Mach lines through rs. On the line £, choose points at equal 
intervals of @ and interpolate |(0Z/d0x) + (iv/V)Z} to these 
points. Then compute the integrand of Eq. (9) at these points 


and finally the integral (9). The computations may be arranged 


as follows: 


Let 


an ae 


u —— 


the part of [w| which gives values of |(0Z/dvr) 4+ 


(iv/V)Z} on the span-wise line & 


FORUM 641 

[7e"* | = a matrix for interpolating from the point on the 
line — to the points equally spaced in 8 relative 
to the point rs 

(Tks j = a matrix for integrating over @ on the points 
equally spaced in @ relative to the point rs 

[°cos **| = matrix of the quantities cos |(v/B%a) (x f) 
sin @{| at the points equally spaced in @ 

Then, 


loa ] { 


irs irs : . 
[ Je®”* | [cos ** | [we] {Zg} 


J [3% iv | fv ) ye 


+ -~Z> cos < (x £) sin 06> d@> = 
| 1  B2a f 


(10) 


When }(0Z/O0x) + (1v/V)Z} is constant over the interval of inte 
gration from 0 to z, the integral (9) is 


wJ ot (v/Bta) (x — E)} [w'] 1 Zp} (11) 
At the point — = x, the integral (9) is 
r[w*] 1Z3} (12) 


Having computed the integrals (9) relative to rs on segments of 
spanwise lines at £ intercepted by the forward Mach lines through 
rs or wing plan form, arrange them in a matrix 


la") {Zp} 





»M 2a > ‘ 
(vM/B°) (~ — £) for the corresponding span 


ETS). 


Compute the e 


wise lines and arrange them in a matrix [" 


{0 


construct a 
matrix [J:**] for integrating over all the — on the wing in the 
Mach cone forward of the point rs. Then, finally, for the point 
rs 


o” = [Jee] [| [e*] |Z} (13) 
Br 
Arrange all the ®” in a matrix form 
[@*] = —(V/Br) [P|] 'Z~,} (14) 
If the pressures at the points rs are required, then 
= [pV?/Br) [w] [P|] (Ze! (15) 


f 
1 Prs} 


For computing the Qag, construct the matrix 


LS] = [fw] (Fe"] — fared} [Jy (16) 
Then, 
Vas nities (pV? Br) |Za| |S] [P| \Z~} (17) 
where 
|'w | = transpose of the complex conjugate of [w| 
drt | = matrix corresponding to |'w] but with zero 
elements, except for the trailing edge points 
where the element is unity 
°J,|, |°\J,| = diagonal matrices of terms which, if multiplied 


as a prefactor with a column and summed, 
y of the function for 
the the 


give the integral over x, 


which the column gives values at 


points of x, y 


The matrices [S| and |w] are, except for iv/V in the diagonal 
elements, constants independent of Mach Number and frequency 
but depending on plan form. The matrices [J] are indepen 
dent of Mach Number and frequency and depend on plan form 


are constants but are dependent on the range 


The matrices | J” 
of @ and are proportional to the range of n corresponding to range 
of 6 Thus, al 
though the task of preparing the problem for the computing ma 
chines is formidable, much of it, when once done, applies for all 


The sets of numbers x — are also constants 


Mach Numbers and frequencies for which Eq. (1) is applicable 


for a particular plan form 


A major task in proceeding from the statement of the problem 
given here to programs for digital computing machines lies in de 
vising integrating means for equally spaced data but oddly placed 


limits. The number of cases for a particular Mach Number and 
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plan form is large; all of the cases change for different Mach 
Number or plan form. 
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Remarks on ‘‘Correction for Aspect Ratio of 
the Indicial Lift Functions of Wagner and 
Kiissner”’ 


S. F. Shen 


Department of Aeronautical Engineering, University of Maryland, 


College Park, Md. 
May 29, 1952 


IX A RECENT NOTE, Scanlan! pointed out that, knowing the 
aspect ratio correction on the Theodorsen function C(k), such 
as Reissner’s o(k), one can deduce, by the well-known transform 
relationships first proposed by Garrick, the appropriate indicial 
lift functions of Wagner and Kiissner including the aspect ratio 
effects. 

Attention might be called to the fact that Jones? had previously 
employed exactly the same resultant expressions in his theory of 
lifting surfaces in nonuniform motion, although the aspect ratio 
correction o(k) was derived in a different manner. In fact, 
Jones went slightly further to include a correction on the appar- 
ent mass term, which was omitted in reference 1. 
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Hooke’s Gas in Aerothermodynamics 


ing. Ennio Mattioli 
Scuola di Ingegneria Aeronautica, Rome, Italy 
May 19, 1952 


SUMMARY 


It is shown that the use of Hooke’s gas (4 1) is inopportune and that 
it is preferable to limit oneself to Karman-Tsien’s ‘‘tangent gas One 
observes that the ‘‘correction coefficient”’ of the latter is nearly exact in the 


entire flowif Mmaz ~0. 5. 


DISCUSSION 


M*** TRANSFORMATIONS OF AEROTHERMODYNAMICS'’ are as- 
sumed to be adiabatic. This gives, for perfect gases, the 


relation 





pe *% = const. =k, k>O 
One can construct an aerothermodynamics assuming as 
equation of the adiabatic transformation any relation 
f(b, p) = 0 
of a type different from Eq. (1). It will be sufficient to set a¢ 
some of the hypotheses of the perfect gas. In particular, ox, 
consider the gases with straight adiabatic in the plane (1 /p, ; 
pP-fpr = m((1 p) — (1 pi), m = const 
1/pi, p: being value at a state of reference of the gas 
Preserving, also, the definition of the velocity of sound, 
a? = dp/dp 
(differentiation along an adiabatic) one will have for such gases 
a*p? = —m 
and, therefore, we must have m < 0 
Bernoulli’s equation, 
(dp/p) + VdV=0 
integrated along a stream line gives 
a? — V2 = const. 
Dividing Eq. (6) by Eq. (5), one obtains 
(1 — M?)/p? = const 7 
with 14 = Mach Number. This result was found for a “tanger 
gas”’ by Tsien? and von Karman,’ who substituted for the adj 
abatic line of the perfect gas its tangent at the point 1/p,, ; 
corresponding to conditions at infinity. They utilized the co 
stancy of the first member of Eq. (7) in order to obtain thei 
noted “‘correction coefficient.”’ 
For this ‘“‘tangent gas,”’ therefore, Eqs. (5) and (6) are satisfie 
with 
m = —p;a;?> = —pi*(dp/dp); = yPipi < 0 8 
For the perfect gas, one has, in place of Eqs. (5) and (6 
noted# 
2,1 on 
a“p = const 


a? + [(4 1)/2]V?2 = const f 


7 1 2/(1 Y 
(: t ur) p” = const 7 
2 


Hooke’s gas is defined as a hypothetical perfect gas having 


1 5, 6 


i 


With y = 1, Eq. (1) becomes a straight line, 
p = R(1/p) 


Therefore, for this gas, Eqs. (5) and (6) hold true with m = 
However, since Eqs. (5’) and (6’) must hold true also, wit 
y = —1, it follows that the relationship of Eqs. (5) and (6) vali 
for the gases with straight adiabatics, and in particular for Kar 
man-Tsien’s tangent gas, can be obtained formally from thos¢ 

the perfect gas by putting in the latter y = 1. However, tl 


wit 


adiabatic equations of these gases are of the type of Eq. (3 


m < 0 and cannot result in any way from the corresponding 


formula for the perfect gas 
For Hooke’s gas, one has 
a= (p/p) 0 
therefore, one cannot speak in this case of the velocity of sou! 
In the numerous problems in which the speed of sound figures 
the square, one can use formally Eq. (10), but it is rather io} 
portune to introduce in the calculations an expression a? to 1! 


cate a negative number 


The adiabatics of Hooke’s gas pass through the point of inter 


section of the axis with the slope k > 0; therefore, such 


can never be used as a “‘tangent gas’’ of a perfect gas 
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In addition, Eq. (9), since k O, shows that in an adiabatic 


expansion of Hooke’s gas the pressure and volume increase simul- 
sneously ; this is in striking contrast to our usual concept of 
hysics in this connection. 

For these reasons it seems inadvisable to use Hooke’s gas even 
asasimple mathematical model. 

If one develops in series the numerator of Eq. (7’) and sup- 
noses the powers of ./ with exponent > 4 to be insignificant coin- 
pared with unity, one obtains again Eq. (7) as an approximation 
to Eq. (7’). For M = 0.5 and ¥ 1.4, the error committed in 
stopping the development in series at the second term is about 
} per cent and decreases to about 1 per cent for (0.4. 

Therefore, if the Mach Numbers of the entire field of speed do 
0.5, the Karman-Tsien coeffi 


not exceed much more than 
ient deduced from Eq. (7) can be considered sufficiently exact 
for all the points of the flow, from those at stagnation to those 
of maximum velocity. This represents a sure advantage over 
the linearized theory which is in error near V = 0 
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Reply to Comments on ‘‘Horizontal Tail Loads 





laving 





in Abrupt Pull-Ups from Level Flight’’ 


L. Decker 
The Glenn L. Martin Company, Baltimore, Md 
May 26, 19592 
in 


Loads 


Tail 
indicate that 


, COMMENTS ETKIN! on ‘Horizontal 
Abrupt Pull-Ups from Level Flight’’? 


points relative to the subject require clarification 


OF 


several 


The general tail-load problem involves the determination of 
the air loads that are critical in the design of the structure. The 
critical design loads fall into four broad categories: (1) maxi- 
mum overall tail load, (2) maximum torque load, (3) maximum 
load on elevator, and (4) maximum load on stabilizer. Of these 
four critical load types, only the first is critical in pull-up or 
pitching The 
maximum horizontal tail torque, for instance, is usually critical 


pushover maneuvers involving acceleration. 
ina steady turning maneuver at limit load factor, with a forward 
enter of gravity location 

Etkin raises three objections to the fourth-degree polynomial 
n vs. t relation used in reference 2 to represent the pull-up that 
gives the maximum tail load. These objections are: 


] di,/d, + 


2) dé,/d; during the entry to the pull-up is too small in ab 


0 at ¢ = O, the start of the maneuver. 


solute magnitude relative to dé,./d; at Mmaz 


6) 


hick 
ign 


If the solution is extended beyond the time for nmar., 








er tail loads are computed. 





With respect to the first objection, it is obviously theoretically 
true that, if a pull-up maneuver is initiated from stabilized level 


flight, d5./d, = 0 at the start of the maneuver. However, as may 


Se seen from flight records, of which those in reference 3 are typi 


' ; : : ; 
cal, the time interval required to develop an elevator deflection 
rate of the magnitude required by Eq. (13) of reference 2 


is ex 


ly small and, for engineering purposes, is negligible. Con 
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siderable success has been achieved in the calculation of maneuv- 
ering tail loads in the past‘! by assuming the elevator motion 
producing the maneuver to be made up of several constant dé,/d; 
segments. This implies infinite acceleration of the elevator at 
each end of a segment. 

Concerning the relatively low elevator rates required by Eq. 
(13) of reference 2 during the early phases of the maneuver, it 
was realized by the writer that larger values of positive pitching 
acceleration could, and sometimes would, exist in the early part 
of the pull up This is not considered objectionable, since this 
phase of the maneuver does not result in a critical loading condi 
tion of the horizontal tail. Flight records, such as are shown in 
reference 3, indicate that the tail load in the pull-up entry is but 
This 


is particularly true at aft c.g. locations where the maximum tail 


a fraction of the load existing at the maximum load factor 


load occurs. 

It was pointed out by the author in reference 2 that the quartic 
polynomial vs. ¢ relation resulted in a tail-load variation that 
was increasing at the time for ”,»,¢-. and that this relation was not 


to be used after mar, Was attained. It was also pointed out that, 
(which is the largest value of 


the 


if the pitching acceleration at My« 
@ occurring during the maneuver) is held throughout re- 
the tail load diminishes during the recovery. It 
is Maintained throughout the 
at If 


is to be devel 


covery, is also 
true that, if d6,/d; existing at 7, 
recovery, the maximum tail load will stil! occur a a 
tail load that exceeds the value existing at n, 
oped, it is necessary to have available a value of d6,/d; which ex 
at variation 


of the maximum value of dé, 


ceeds that existing : Since, in reference 2, a 
d, for use in the tail-loads calcula 
tions has been proposed, it is inconsistent with the general method 
to use a higher value of d5,/d; in the recovery. From a mathe- 
matical standpoint, the type of pull-up recovery may be defined 


by the condition: fort 2 tmaz., 


(d8,/d¢)auring recover dé 


rhe type of recovery so defined will be consistent with the pull-up 
of reference 2. 

Inspection of considerable flight-test data, for which reference 
3 shows representative results, indicates that the maximum tail 
load occurs simultaneously with the maximum load factor in an 
abrupt pull-up maneuver Normal piloting technique results 
in a type of maneuver for which the tail load decreases rapidly 
after Mma It should be noted that there is 


an infinite number of maneuvers in which an airplane can reach 


has been attained 
a given load factor in a given time as the piloting technique 1s 
altered. From a practical standpoint, the prime requirement of 
the method to compute the maneuvering tail load is that it pre 
dicts the maximum positive pitching acceleration that is likely to 
occur at the limit factor 
checked by the writer, the method of reference 2 has yielded ex 


load For airplanes that have been 


cellent results 
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Fig. 1 shows the variation of load factor, change in elevator 
deflection from trim, and change in tail load from the level flight 
trim condition versus time for an abrupt pull-up performed with 
an SB2C-1 airplane. Flight data from reference 3 and calcula- 
tions by reference 2 are compared. The maximum load factor 
of 5.8 and time to attain it of 1.75 sec. were selected for the calcu 
lations to duplicate the flight example. A maximum increment 
in tail load during the pull-up of 3,330 Ibs. was computed from 
reference 2, and this shows nearly perfect agreement with the 
flight results, which indicate approximately, 3,300 Ibs. delta tail 
load. The constant value of dé,/d; in the recovery produces a 
load variation that closely checks that in the flight results 
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Note on the Calculation of Transport Properties 
of Gas Mixtures* 


Eldon L. Knuth 

Guggenheim Fellow in Jet Propulsion, Daniel and Florence 
Guggenheim Jet Propulsion Center, California Institute of 
Technology 


May 23, 1952 


r. w A RESULT OF RECENT THEORETICAL STUDIES on transport 
properties of gases by Hirschfelder, et al.,! and others,” the 
dynamic viscosity (u), the thermal conductivity (), and the me- 
chanical diffusivity (D),as well as the dimensionless ratios Prandtl 
Number (Pr = C, u/d, where C, = specific heat at constant pres- 
sure), Schmidt Number (Sc = u/pD, where p = density), and \/- 
pC,Dt may now be calculated rapidly and with a reasonable degree 
of accuracy both for pure gases and for mixtures of gases, provided 
that the temperature of the gases is not too low. This happy 
state of affairs does not seem to be realized by aeronautical engi- 
neers, and it is the purpose of the present note to call attention to 
it. 

The greatest error in these calculations arises in the determi- 
nation of A, where the Eucken approximation is used. This ap- 
proximation involves the assumption that thermal equilibrium 
is established between the translational and internal degrees of 
freedom of colliding molecules, an assumption that describes 
physical facts most closely at elevated temperatures. How 
ever, even at moderate temperatures, the calculated values of 
\ are not in error, in most cases, by more than 2 per cent accord- 
ing to reference 1. 

As an illustration, calculations have been made for Pr of dry 
air (SO per cent nitrogen and 20 per cent oxygen) and for Sc and 
\/pC,D for water diffusing into air, assuming that C, (see refer- 
ence 4), uw, p, and Xd are those of dry air and that no diffusion by 
temperature gradient occurs. The results of these calculations 
are presented in Fig. 1. It is expected that the methods of cal- 
culation to which attention has been called in this communica- 
tion will be of value, for example, to engineers interested in 

* The calculations required for Fig. 1 were made with the assistance of 
E. Zukoski, G. McLeod, F. Hartwig, and J. Schroeder 

¢ Although a Klinkenberg and H. H. Mooy have given the name ‘‘Lewis 
Number”’ to \/pCpD in reference 3, this nomenclature is not generally used 
in the present literature 
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boundary-layer calculations or the determination of transport 
coefficients in combustion chambers 
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The Prandtl Membrane Analogy for 
Temperature Fields with Permanent Heat 
Sources or Sinks 


Paul J. Schneider 
Department of Mechanical Engineering, The State University of lowa 
June 4, 1952 


i pur PRANDTL MEMBRANE ANALOGY is frequently employed 
for experimental evaluation of mechanical stress fields as 
exemplified by the well-known Saint-Venant torsion problem for 
prismatic bars.'~’ The analogy has a number of logical exten 
sions, one of which includes experimental solutions for steady 
state two-dimensional temperature distributions satisfying the 
Poisson partial-differential equation—i.e., temperature fields m 
regions with permanent internal heat sources or sinks. Approx! 
mate numerical (iterative) solutions and exact eigenfunction 
solutions for these temperature fields are possible to obtain only 
when the boundary geometry of such regions is of elementafy 
form—circular, elliptical, and rectangular. On the other hand, 
membrane solutions encompass all cases of nonregular boundary 
geometry which are otherwise numerically awkward and mathe 
matically intractable. 

The ordinates 2(x,y) of the surface of a homogeneous mem 
brane (such as the soap film in Fig. 1) spread between appropti- 
ately designed boundaries and dilated by a uniform static pres 
sure are to be accurately measured, and these measurements 
are then to be reinterpreted in terms of the temperature field 
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Side view of a dilated soap film on a circular boundary 


Fic. 1 


(x,y) in a region having the same two-dimensional boundaries. 
The partial-differential equation of this membrane surface is 
given in reasonable approximation by 
oO? °? p 
=~ A229) + -~ Ax) + == () (1) 
Ox? oy? T 
providing the maximum slope of the film does not exceed say 
30°. The constants p and 7 are identified, respectively, as the di- 
lation pressure and film tension. Now, the characteristic par- 
tial-differential equation satisfied by the steady-state temper- 
ature ¢(x,y) in a region of constant thermal conductivity kp gen- 
erating or absorbing uniform heat q is given by 
oxy) +—Hxy) +2 = 0 (2) 
x,y x,y = 2 
Ox? oy? ko 
The apparent mathematical identity exhibited between Eqs. 
(1) and (2) leads directly to an analogic relation between the 
analogs 2(x,y) and ¢(x,y) of the form 
[t(x,y)] = (Ro/qo)t(x,y) = El2(x,y)] = (7/p)2(x,y) (3) 
or specifically 
U(x,y) = (qo/ko) (1r/p)2(x,y) 
Identifying go and ko for the region gives the algebraic relation in 
terms of the ratio +/p for converting experimental membrane 
data 2(x,y) to analogous temperatures (x,y). 
In all cases the experimental ratio p/r can be evaluated with 
sufficient accuracy from an auxiliary dilated soap film on a circu- 


lar boundary. The ordinary differential equation in cylindrical 
coordinates for the surface of this auxiliary symmetrical film of 


radius 7; is 
d? ld p 
- Za + — Za + = 0) 
dr? r dr 
which has the particular solution 
Sa = (1/4) (p/r)a (%a,? — 1a?) 
satisfying the required boundary conditions 2,(r1) = 0 = (dz + 
dr);.o. Herewith in terms of the easily measured maximum 
ordinate z,,at r = 0 
(b/r)a = (4/ra,?)Zap (4) 
Then, simply, (p/r) = (p/r)aif p = pa. In this way any bound 
ary temperature may be duplicated by building up z on the main 
boundary without introducing additional difficulties in evalu- 


ating p/r. For the most elementary case with ¢(n,4) = 0 on 
the boundary (7,4), then also 2(7,4) = 0 for the membrane 
boundary. 


These techniques encompass all cases with known distributions 
of boundary temperature regardless of the geometric complexity 
of the boundary. In support of this, experimental evidence has 
been obtained for a number of elementary examples to verify 
that such analogic solutions are in substantial agreement with 
corresponding exact analytical solutions. 
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Note on the Theory of Hurricanes 


Theodore Theodorsen 
Hq., Air Research & Development Command, Baltimore, Md 
April 30, 1952 


SE pnees ARE CERTAIN INTERESTING CONSEQUENCES of the new 
theory of turbulence* applying directly to the theory of hurri- 
canes. These are: Hurricanes are created by the ‘‘collapse,”’ 
in a universal pattern, of a large-scale boundary layer of consider- 
able horizontal extension. Contrary to the situation in a small 
scale boundary layer, where the more or less symmetrical ‘‘horse- 
shoe”’ develops, only the counterclockwise branch of the vortical 
motion will generally appear in the northern hemispheres, since 
there is no longer complete symmetry because of the action of the 
Coriolis forces. Evidently, the required boundary layer develops 
by the passing of an air mass over extensive smooth surfaces such 
as the oceans or plains. The boundary layer becomes superla- 
tively unstable, or, one may say, there obtains under such cir- 
cumstances a greatly delayed transition. Then, let us say that 
somewhere on the surface below the air mass an obstruction 
occurs; the transition takes place violently, since the lower criti- 
cal may have been passed by a thousand-fold. Turbulence de 
velops behind the obstruction and produces an extensive retarded 
air mass. This situation is then mathematically equivalent to 
the existence of a portion of a large vortex extending from the 
ground over the upper boundary of the disturbed air mass and 
down to the ground on the other side of the retarded air mass 
As the Reynolds Number is far above the critical, the vortex de- 
velops with great intensity by the universal turbulent mecha- 
nism; the horizontal portion of the vortex is moved downstream 
and upwards as shown in the new theory of turbulence. Clearly, 
the counterclockwise arm will develop most easily and, at least 
for large dimensions, will appear alone. Thus, the creative 
mechanism of hurricanes is entirely similar to the mechanism 
of turbulence. 

* The Second Midwestern Conference on Fluid Mechanics, March 17, 
1952, Theodore Theodorsen, ‘‘Mechanism of Turbulence.”’ 


On Boundary-Layer Temperature Recovery 
Factors 


B. des Clers and J. Sternberg 
Ballistic Research Laboratories, Aberdeen Proving Ground, Md 


May 15, 1952 


| A PAPER IN THE January, 1952, JOURNAL OF THE AERONAUTI- 
CAL ScIENCES, “Recent Investigation of Temperature Re 
covery and Heat Transmission on Cones and Cylinders in Axial 
Flow in the N.O.L. Aeroballistics Wind Tunnel,’ Eber presents 
some experimental data for boundary-layer temperature re- 
covery factors. The results of a series of recovery factor meas- 
urements in the B.R.L. flexible nozzle tunnei using cones and 
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Fic. 1. Spark shadowgraph showing transition at Re = 


3.2 & 108 Surface Mach Number = 3.40. 
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Fic.2. Cone recovery factors — 10° included angle cone surface 

Mach Number = 2.18. R, conditions at an insulated wall; 

©, conditions at outer edge of boundary layer; 0, conditions at 
stagnation 


cone cylindets do not agree with Eber’s results for the transitional 
and turbulent recovery factors 

Laminar, transitional, and turbulent boundary-layer flow was 
obtained on a 14.5-in. 10° included angle cone at surface Mach 
Numbers of 1.40 to 3.40 
and thermocouple wires attached to small plugs spaced along an 


The model was constructed of wood, 


element of the cone were used to measure the recovery temper- 
ature. The tunnel was operated at the test condition for 20-30 
min. to allow the model to reach its equilibrium temperature 
Shadowgraph pictures, such as Fig. 1, corroborate the identifica- 
tion of the transition region made from the temperature meas 


A typical set of measurements is shown in Fig. 2 


urements. 
The absolute accuracy in r is estimated at +0.007; the relative 
accuracy can be seen from the graph. Our conclusions are as 
follows: 

(1) The turbulent recovery factor has a value of about 0.88 


and does not significantly change with Mach Number, at least 


1952 
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up to .M = 340. Eber obtained a value of r = 0.92 at y ~ 
2.86 andr = 0.97 at M = 4.25. 

(2) The recovery temperature rises smoothly from the am; 
nar value to a maximum not more than 1 per cent higher thay 
the turbulent value. At M = 2.87, Eber obtained peak valye 
about 5 per cent higher. We believe that the peak that we ob 
tained (which is not present in all cases) may be due to non 
equilibrium turbulent profiles at the end of the transition region, 
since the temperature and velocity profiles at the end of the 
transition region should depend on the mechanism of transition 
in any given case. 

The measurements on a cone cylinder where transition 
curred on the cylinder agree with the above results. * 

It is possible that Eber’s measurements do not represent 
equilibrium conditions. The recovery temperature at the meas 
uring station may be affected by heat transfer from the model t 
the boundary layer ahead of the measuring station, since with the 
short running time available it is probable that the large thermal 
inertia of ‘his model prevents it from coming to temperatur 
equilibrium. As noted by Eber, his schlieren photographs shoy 
that the boundary layer is fully turbulent where the recovery 
factors are constant. Since the boundary layer is laminar op 
the cone, it is possible that, over most of the region downstream 
of the shoulder which he labels transitional, the boundary layer 
may still be laminar. 


* Measurements just published at the N.A.C.A. Ames Aeronautica 


Laboratory by H. A. Stine and R. Scherrer (N.A.C.A. T.N. No. 2664) are ir 


complete agreement with our concllisions (1) and (2) 


Secondary Flow in Cascades 


H. G. Loos and J. Zwaaneveld 
National Aeronautical Research Institute, Amsterdam Holland 
May 26, 1952 


IT A NOTE BY STEPHENSON,! the vortex field behind a cascade 
with an infinite number of blades in a nonuniform flow is in 
vestigated. Stephenson calculates the secondary rotation be 
hind the cascade and adds the trailing vortices. 

In a later note, by Eichenberger,? it is suggested that the add 
ing of the trailing vortices is incorrect. 

Eichenberger’s arguments are the following: The flow through 
the cascade can be idealized as a flow through a curved channel 
The secondary vortices then give a flow field as shown in Fig. | 
and on the two sides of the wall A-A there is a different velocity 

So, the wall A-A can be replaced by a vortex sheet. On the 
other hand, a vortex sheet with vortices along stream lines can 
be replaced by a thin wall. So, the resulting velocity field due to 
the nonuniform inlet flow can be said to be induced by the second 
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ary vortices alone if the trailing vortex sheets are replaced by 
walls 

Eichenberger concludes that the velocity field of the trailing 
vortices must not be taken into account. 

Making Eichenberger’s conclusion, one must say explicitly 
that the velocity field of the secondary vortices must be deter- 
mined with the boundary conditions on the walls, which replace 
the trailing vortex sheets. However, the shape of these ‘‘walls”’ 
is not yet known, since, therefore, the complete flow field must 
be known. Therefore, this point of view seems to be fruit 
less. 

To obtain an approximate solution, the secondary vortices 
and the trailing vortices are considered to be transported by the 
mean flow, so that the induced velocities can be calculated. 

Instead of the direct connection (suggested by Eichenberger) 
between the secondary vortices and trailing vortices, there exists 
another, more complicated connection via the downwash field of 
the secondary vortices and the Kutta condition at the trailing 
edges of the blades. 

As Eichenberger remarks, the vorticity vector does not re 
main constant through the cascade, in general. However, for a 
cascade with an infinite number of blades, it is impossible to de- 
termine the secondary vortices with Squire and Winter’s‘ as- 
sumption for the basis flow (¢/r = constant). 

In problems like this, where the secondary rotation in lightly 
nonuniform flows is investigated, the relation of Hawthorne* 


» 
—2\m af, (d6, gq?) 
1 
appears to be useful. 


In this relation, g and 6 are the modulus and the angle of the 
velocity vector in the mean flow, while — and 7 are the compo 
nents of the rotation in the mean flow direction and perpendicular 


(é q)e (é qu = 


toit (Fig. 2). 

The integral must be evaluated along the stream lines of the 
mean flow 

Using the fact that the velocity component u (Fig. 2) is con- 
stant in the whole field, the integral can be calculated directly 


“ d0 l , 
_= cos? 6 dé = 
vi ¢ u> J; 
l : 
sin 20, ) 


1 
|. A) + — (sin 262 
2u? 2 
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and is independent of the distribution of the circulation along the 
blades. Hence, 


Bc. ; 
0) + > (sin 26. sin 26;) 


2 m 
oe = (6 
cos 6; COS 62 
There is an interesting interpretation for the integral 
d0/q* 
1 
in Hawthorne’s relation. Mapping the (x, y) plane (Fig. 2) on 
the plane with polar coordinates (1/g, @), a figure for the stream 
lines is obtained, where areas with the origin as pole are corre 
sponding with the value of the integral 
d0/q* 
1 ; 

Therefore, the secondary vorticity £ is, apart for a factor, equal 
to areas in the plane of the reciprocal hodograph. This plane is 
represented in Fig. 3. 

The curve (a) corresponds with a typical stream line in the 
configuration of Squire and Winter‘ (¢ = constant), while curve 
(6) corresponds with a stream line in the Stephenson’s cascade 
with an infinite number of blades (u = g cos @ = constant). 
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The von Karman Integral Method as Applied 
to a Turbulent Boundary Layer 


Nguyen Van Le 
Gas Turbine Laboratory, Massachusetts Institute of Technology 


May 2, 1952 


2 pm PRESENT two-dimensional 
layer equation 


T) dé ( é* ) 6 dl 
is 1. 2 : = 0 (1) 
pU? dx 4 U dx : 


was obtained by Gruschwitz (reference 1, p. 118) by integrating, 
across the boundary-layer thickness, the Navier-Stokes equations 
The same result was later obtained by 


incompressible boundary 


simplified by Prandtl. 
von Karman by considering a control volume 

As derived, this equation is applied to laminar, as well as turbu- 
However, recent investigations® * show 


lent, boundary layer. 
that the equation fails to give suitable results in cases of positive 
pressure gradients. It considers only the mean flow. The ob- 
ject of this note is to present a simple derivation of the boundary- 
layer equation in a turbulent incompressible two-dimensional 
flow, taking into consideration the time components of the velo 
It will make use of the von Karman integral method and 


city. 
(See 


of the results available from the derivation of Eq. (1) 
Fig. 1.) 
Application of the momentum theorem to this control volume 


gives 


° 
5 Op ) o o 
~t > u*® dy l / u dy 
p Ox p Ox Jo ox Jo 


Upon taking the time average, one obtains 
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a) 


6 Op T) 0 7 ma?) 
“a is = u2dy — l 
p Ox p Ox 0 Ox 0 


ua? + u’u’, Eq. (2) becomes 


u dy 


ui +u’ and x42 = 
6 Op T) ce) . " z . 
-— — = u?dy — l 
p Ox | Ox 0 
re) 


ox Jo 


With u = 


u dy + 


p Ox 


Because of the fluctuation components, 0p/Ox can no longer be 
calculated as usual from the flow outside the boundary layer 
With the usual assumptions for boundary-layer study, the dy 
namic equation in the y direction gives (reference 1, p. 337) 

1 Op fe) a 


- = ol oo 
Ui 


u 
p Oy Ox oy 


Assuming v’v’ = Oat y = 0, y = 6, one obtains from Eq. (4) 


p= pete 


| Op 1 dps 


p Ox 


p dx 


"a (: a 
uv tt v 
Jo ox \ox oy 


Thus 0p/Ox is a variable in y through the last term. As under- 


stood from the momentum theorem applied to the previous con 
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trol volume, (5/p) (0p/Ox) in Eqs. (2) and (3) is the mean Value 
in the y direction; therefore, 


6 Op 6 dps o 
7 =o +8 
P OX} mean value in p dx 9 oO 


y direction 


6 7] a ( a 
u'y’ 
0 0 ox \ox 


Putting Eq. (5) in Eq. (3), dividing by U?, and rearranging, we 


u'v’ dy 
awa 
T ed y dy (5) 


get 


56 dps TT) ra) 2 a4 
-dy = 


+ + = 
pl? dx pl? ox f x l 
6 
U* Jo 


| 2 
? Jn Ox e 
: 2 (2 wo + 2 we") aya 
I uv v'v’ | dy dy 
l? Jo 9 ox \ox “ oy ydy @ 


The left-hand member of this equation is recognized to be the ong 
obtained with the mean value of the velocity by von Ka4rmég 
With the usual definition 


u'u' dy 


(reference 4, p. 132). 


Displacement thickness 6* 


Momentum thickness 6 


it yields the classical equation, E 


2) 6 dl 1 l [ re) 
o U dx U2 0 Ox 
1 


mediately as 


dé (= int d 
ye u'u' ay 

re) 

zs vay 


u'v ae 


° 6 
5 oz) o 1 "6 "y a Py 
: u'y’ dy + — 
U2 J 0 Ox U? 0 0 Ox \Ox 


Together with reference 3, this study gives a parallel result to 
those of Gruschwitz and von Karman mentioned previously. 
Experimental work is now directed toward the correlation of the 


time components of the velocities 
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